Fakultat fur Physik und Astronomie 

Ruprecht-Karls-Universitat Heidelberg 



Diplomarbeit 

im Studiengang Physik 

vorgelegt von 

Tobias Aurelius Knoch 
aus Mannheim 


4. Februar 1998 








Dreidimensionale Organisation 
von 

Chromosomen- Domanen 


in 

Simulation und Experiment 


Die Diplomarbeit wurde von Tobias A. Knoch ausgefiihrt 


am Deutschen Krebsforschungszentrum Heidelberg 
in der Heidelberg 3D Human Genome Study Group 
als Teil des Deutschen Human Genom Projekts 


unter der Betreuung von 
Herm Prof. Dr. Jorg Langowski 


sowie von 

Herm Prof. Dr. Dr. Christoph Cremer 
Institut fiir Angewandte Physik 




Zusammenfassung 

Um Aussagen iiber die dreidimensionale Organisation von Chromosomen in der Interphase zu machen, wurden 
Monte-Carlo - und Brownsche-Dynamik-Simulationen von Chromosomenterritorien mit Hilfe eines 
Polymeransatzes durchgefuhrt. Daraus wurden Observable (verschiedene Abstandsverteilungen zwischen 
genomischen Markem, Aufenthaltswahrscheinlichkeit und Dichte) berechnet, die fiir die Vorhersage und 
Interpretation von Experimenten dienen. Damit konnen Modelle der chromosomalen Organisation, insbesondere 
das Random-Walk-Giant-Loop- (MLS) und das Multi-Loop -Subcompartment- Modell, unterschieden werden. 

Fur die Bestimmung von Abstandsverteilungen zwischen genomischen Markem der Prader-Labhart- 
WillUAngelmann (PLWS/A) Region wurden Experimente mit einem konfokalen Laser-Scanning-Mikroskop 
(CLSM) durchgefuhrt. Zellkeme von Fibroblasten wurden mit schonender Fluoreszenz-in-situ-Hybridisierung 
(FISH) prapariert. Vber die Signal-Brightness, das Signalvolumen, die Intergrierte Fluoreszenz Intensitat (IFI) 
und den Abstand homologer Chromosomen wurden ebenfalls Aussagen gemacht. 

Vergleiche zwischen Simulation und Experiment ergaben starke Hinweise fur das MLS-Modell, sowie eine 
unterschiedliche Organisation zwischen mutterlichem und vaterlichem Chromosom 15. 

Im Hinblick auf zukunftige komplexere Fragestellungen der Gesamtorganisation des Zellkems stellte sich die 
Theorie der Fraktale fiir neue Analysemethoden als geeignet heraus. Die Ergebnisse der Simulationen zeigte 
dabei multifraktales Verhalten. Erstmals konnte damit die Fraktalitdt und Multifraktalitdt einer zelluldren 
Organisationsstruktur gezeigt worden sein. Damit besteht nun die Moglichkeit aus der fraktalen Struktur 
Aussagen iiber Diffusioneigenschaften und andere Organisationsmerkmale im Zellkem zu treffen. 


Three - Dimensional Organization of Chromosome Territories 

in 

Simulation and Experiment 

Abstract 

To draw conclusions about the three-dimensional organization of interphase-chromosomes, Monte-Carlo and 
Brownian-Dynamic simulations of chromosomal territorries were conducted, using a polymer model for the 
chromatin fiber. From the simulated configuration ensembles observables (various distance distributions between 
genomic markers, radial space distribution and density) were calculated and used for predictions and 
interpretation of experiments.lt was possible to distinguish between different models of chromosomal 
organization, especially between the Random-Walk-Giant-Loop- (RWGL) and the Multi-Loop-Subcompartment- 
model. 

For the determination of distance distributions between genomic markers of the Prader-Labhart-Willi/Angelmann 
region (PLWS/A) experiments were made with a confocal Laser-Scanning-Microscope (CLSM). Cell nuclei of 
fibroblast cells were prepared with soft Fluorescent-in-situ-Hybridization (FISH). Conclusions can be drawn 
from the signal-brightness, signal-volume, Integrated Fluorescent Intensity (IFI) and the distance between 
homologous chromosomes. 

Comparisons between simulation and experiment lead to strong hints for the MLS-model, as well as ab different 
organization of maternal and paternal chromsomel5. 

With respect to future aspects of more complexe questions on the total organization of cellnuclei the theory of 
fractals proved very worthwhile for new methods of analysis. The results of the simulations thereby showed multi- 
fractal behaviour. Thus for the first time the fractal and multi-fractal behaviour of a cellular-organization- 
structure was determined. Consequently possibilities exist to draw conclusions from the fractal structure for 
diffussion behaviour and other organizational patterns of the cell nucleus. 
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EINLEITUNG 


Sen mehr als dmhundert Jahren steht die Zelle im Blickpunkt der Forschung 1839 konnte 
Schwann m,t Hilfe verbesserter Mikroskopietechniken die Zelle als Grundeinhelt ^ er Pknren 

J ? \ VTJT Tul 5l l b ^ eutendste Struktur innerhalb der Zelle stellt der schoT^ 

17.Jahrhundert von Robert Hook entdeckte Zellkem dar. 

Der Zellkem enthalt ein langes fadenformiges Molekiil, die Desoxyribonukleinsaure (DNA) 

Snkemslrffi r T,?"' Ausgestreckt ist die gesamte DNA eines menschlichen 

ZeUkems ca. 2,6 m lang. Em Zellkem besitzt aber im allgemeinen nur einen Durchmesser von 
etwa 10 lira, so daB die DNA in drei Stufen kompaktiert wird, die Nukleosomen, Chromatin 
und Chromosomen genannt werden. Wie Untersuchungen der letzten Jahre gezeigt haben ist 
diese dreidimensionale Organisation des Genoms wesentlich fur die Funktion des Zellkems 
Inwieweit die Translation und die fur die Zellteilung benotigte Replikation von dieser 
dreidimensionalen Organisation beeinfluBt werden, liegt jedoch bis heute weitcehend im 
Dunkeln. 


Aufgrund verschiedener Beobachtungsmethoden wurden seit Beginn des Jahrhunderts 
verschiedene ModeUe fur die raumliche Struktur von Zellkemen wahrend der Interphase, in 
der Gene transkribiert werden, vorgelegtrRabl und Boveri postulierten 1885, bzw. 1909, 
zunachst eine Trennung der Chromosomen in Territorien [Ra 85, Bo 09], Aufgrund von 
elektronenmikroskopischen Aufnahmen stellte Comings 1968 ein Modell auf, bei dem sich 
Chromosomen zufallig und ungeordnet uber den ganzen Zellkem erstreckten [Co68], 

Durch Mikrobestrahlungsexperimente konnte jedoch in den 70er Jahren durch Cremer gezeigt 
werden, daB Chromosomen im Zellkem lokalisierte Strukturen bilden [Cr 74, Zo 79, Cr 82a, 
Cr 82b und Cr 84], 1984 stellten Pienta und Coffey ein Modell auf, das aus ca lOOkbp groBen 
Schleifen bestand. Diese Schleifen sollten sich wahrend der Zellteilung zu Rosetten 
zusammenlagem, aus denen sich wiederum ein Metaphasechromosom aufbaute [PiCo 84], 

Mit der Entwicklung der Fluoreszenz-in-situ-Hybridisierung (FISH) und der konfokalen Laser- 
Scanning-Mikroskopie (CLSM) konnte man 1993 die territoriale Organisation einzelner 
Chromosomen erstmals sichtbar machen [Cr 93, Zi 93]. 1992 benutzte van den Engh FISH, um 
raumliche Abstande von DNA-Sequenzen in Chromosomen zu messen [vdE 92]. Es zeigte 
sich, daB auf Skalen bis l,5Mbp die dreidimensionale Organisation durch ein Random-Walk- 
Modell beschrieben werden kann. Der geometrische Abstand ist bei diesem Modell 
proportional der Quadratwurzel des Genomabstands. Bei groBeren Skalen wachst der 
geometrische Abstand langsamer an. 

Sachs stellte daraufhin 1995 ein Random-Walk-Giant-Loop-(RWGL) Modell vor, das dieses 
Verhalten erklart [Yo 95]. Ein Chromosomenterritorium besteht dabei aus 3Mbp groBen 
flexiblen Chromatinschleifen, die endang eines Random- Walk-Backbones lokalisiert sind. Die 
Verteilung der Chromatinschleifen im Chromosomenterritorium ist bei diesem Modell rein 
zufallig. 

Das von Cremer schon 1993 vorgestellte Inter-Chromosom-Domanen- (ICD) Modell postulien 
dagegen einen weit hoheren Ordnungsgrad [Cr 93, Zi 93]. Hier wird ein hypothetisches 
dreidunensionales Netzwerk zwischen den Territorien angenommen. Genreiche DNA- 
Sequenzen sind an der Oberflache lokalisiert, genarme dagegen im Innem des Territoriums. 
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Einleitung 


Her Lase, die Gesamtheit der experimentellen Daten Zu 
Keines der Modelle ist jed«* ® de shalb das sogenannte Multi-Loop-Subkompartment- 
Schreiben. Munkel en^cke^ « Schleifen auJ chromatin sind dabei zu R osetten 
(MLS) Modell [Mtt 98]. ^ wie derum durch erne Verbindung aus DN A 

zusammengelagert. den ganzen Zellzyklus erhalten und bilden disUnkte 

miteinander verknupft, bleioen 

Einheiten. 

Organisation des Chromatins in Chromosomen ist bis heute 
Die genaue dreidimensionrnc Verfahren besitzen zwar hohe Auflosungen 

unklar geblieben. Etektr °"® Pr aparationstechnik aber die dreidimensionale Organisation’ 
f elsl6re r^ 1 b pi ‘t. V w J h M sind zwar schonend, besitzen dagegen nicht die notwcndige 

hohe Auflosung. 

Hier wurde die dreidimensionale Organisation von Chromosomen in der Interphase simuliert 
Da ^Organisation im ModeU bekannt ist, konnten Observable berechnet werden, d,e fiir die 
Interpretation und Vorhersage von schonenden Experimenten mil gennger raumlicher 
Auflosung dienten.Durch diese Vorgehensweise konnten die oben genannten Einschrankungen 

umgangen werden. 

Die Simulationen wurden mit Monte-Carlo- und Brownschen-Dynamik-Verfahren durch- 
gefiihrtHauptziel der Simulationen war dabei zunachst die Unterscheidung zwischen RWGL- 
und MLS-Modell auf kleinen genomischen Skalen. Neben der Simulation ganzer 
Chromosomentenitorien wurden innerhalb der jeweiligen Modelle verschiedene Parameter, 
wie z. B. die SchleifengroBe oder die Schleifenanzahl, variiert. 

Als Observable wurden aus den Simulationen verschiedene dreidimensionale Abstande fiir den 
Vergleich mit dreidimensionalen Abstanden zwischen genetischen Markem bestimmt. 
Dariiberhinaus wurden weitere Observable gesucht und gefunden (z.B. die radiale Dichte), die 
es erlaubten, genauer zwischen den Modellen zu unterscheiden. 

Um die Simulationen zu verifizieren, wurden Experimente mit hochauflosender konfokaler 
Laser-Scanning-Mikroskopie (CLSM) unter Verwendung der Fluoreszenz-in-situ-Hybridi- 
sierung (FISH) durchgefiihrt. Das Hauptziel bestand dabei in der Messung von Abstands- 
vertedungen zwischen genomischen Markem der Prader-Labhart-Willi/Angelmann-(PLWS/A) 
S T,° SOm 15 * DaS PLWS/A S y ndrora ™ Gegensatz zu anderen genetisch 
Chromosom is" wallrsch einlich auf die dreidimensionale Organisation von 

F “ r * E,perimeme W " d ' — d * d " ** 

vorgenommen Wtileriiin soUte dr 6 " WUrde e ‘ ne Unterscheidung zwischen den Modellen 
werden. Hier konnte gezeiet ** Faltun 8 sstr uktur der PLWS/A Region bestimmt 

wahrscheinlich eine unterschiedlirh^n ^ daB mUtterlic hes und vaterliches Chrorosom 
unterschiedhche dreidunensionale Organisation aufweisen. 

Ze like ms wurde versucht, weiter^if^f^^^ 1111 ^ 11 ^ ezd S Ec h der Gesamtorganisation des 

zugmndeliegenden Prinzipien sollten j melh<Kien und Observable zu finden. D ‘ e 
g ugen Gesetzen und Organisationsprinzini * 16 blolo S ische Fragestellung mit allg erae ‘ n 
Z^ammenhang bringen. P P n aus anderen naturwissenschaftlichen Bereicbe 
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Die Theorie der Fraktale stellte sich dabei als besonders geeignet heraus. Mit ihr lassen sich 
Phinome und Organisationsprinzipien verschiedenster Systeme vergleichen. So konnte erst 
kiirzlich gezeigt werden, daB Skalierungsgesetze zwischen verschiedenen Saugetieren rait nur 
drei Annahmen erklart werden kdnnen [WeBrE 97]. Eine dieser Annahraen besteht darin, daB 
das Transportsysteme (Blutkreislauf) in diesen Saugetieren fraktal ist 

Die Frage war nun, ob es auch im Zellkem ein solch allgemeines Prinzip gibt, das die 
Organisation erklart 

Die Ergbnisse der Simulation zeigen dabei nicht nur fraktales, sondem sogar multifraktales 
Verhalten. Erstmals konnte damit die Fraktalitat und MultifraktaiitSt einer zellularen 
Organisationsstruktur gezeigt worden sein. Erstaunlich ist hierbei die GroBenordnung des 
betrachteten Systems. 

In der Zukunft wird es damit moglich sein, weitere Organisationsmerkmale des Zellkems 
aufzudecken. Zum Beispiel lassen sich aus der fraktalen Struktur Aussagen liber das 
Diffusionsverhalten gewinnen. 


Die Arbeit ist in sechs Kapitel untergliedert: in Kapitel I werden kurz biologische Grundlagen 
und in Kapitel II die Entwicklung der Zellkemorganisation erlautert In Kapitel II werden die 
Simulationen und in Kapitel IV werden die Experimente in ihrer Methodik, Durchfuhrung und 
Ergebnissen beschrieben. In Kapitel V erfolgt ein Vergleich zwischen den Simulationen und 
Experimented SchlieBlich erfolgt in Kapitel VI eine Zusammenfassung der Ergebnisse. 


I Biologische Grundlagen 


Wie schon in der Einleitung erwahnt, ist die Erbsubstanz in vier hierarchische 
Organisationsniveaus gegliedert: DNA, Nukleosomen, Chromatin und Chromosomen.In den 
folgenden Abschnitten sollen diese Organisationsmveaus genauer erlautert werden. 


I.l.DNA 

Miescher isolierte 1869 als Hauptbestandteil des Ze like ms Desoxyribonucleinsaure (DNA) aus 
Eiterzellen und Fischspermien [Br 91]. Erst 1944 wurde aber durch O. T. Avery C M 
MacLeod und M. McCarty durch Transformationsexperimente an Pneumo'kokken 
experimentell nachgewiesen, daB die DNA der Trager der Erbinformation ist [Br 91] 

DNA ist ein fadenformiges Makromolekul, ein Polymer, das aus vier Grundeinheiten den 
sogenannten Nucleosiden, aufgebaut wild. Die Nucleosiden bestehen aus einem Zucker (B-D- 
r-Desoxyribose) und einer von vier heterozyklischen Basen Adenin (A), Guanin(G) Thvmin 
CT) und Cytosm(C). Die Monomere werden durch eine 3 '-5' Phosphodiester-Bindung zu 

Zuck™^ S>h2a£ 3,ymenS1 t rt : i Die Basen sequenz kodiert die genetische Information die 
zucker und Phosphatgruppen erfuilen pnmar struktureile Aufgaben 

*£?££££ SSSi S SSSIS" sclK T hr sul w 

durch James D. Watson und Francis H C Si? 3 VOtl Ront genstrukturanalysen 

Franklin aufgeklart werden f pTu 8 6 L M a „ '/Ta LmUS C PauUn g und Rosali " d * 
Doppelhelix zusaramenlagem. Zwischen Adenin’ ***! Einzelstran 8 e sich zu einer 
zwischen Guanin und Cvtosin drpi w ain . und Thymin werden dabei zwei und 

Einzelstrange in der Doppelhelix eine V ass ® rstoffbrucken bindungen ausgebildet, so daB die 
Die Doppelhelix betoa CAbbXl>. 

Einzelstrange sind dabei be^gHch^hrer Pokf gememSame *** gewunden. Die beiden 
Phosphat- und Desoxyribosereste befinrfon ierun g sdcbtu ng gegenlaufig orientierL Die 

Durch diese dreidimensionale Anordmmo » 1° augen ’ d * e Basenpaare befinden sich innen. 


Durch diese d rc idimeSoS 0 S! n e S rS Base, W befinden sich innen. 

groBere Stabilitat, die sich in einem hnhp Doppelhelix gegeniiber dem Einzelstrang 

AuBerdem kann sich zwischen den Basen ein^Tc aStizitatsmod ul und Torsionsmodul zeigt. 
Die DNA-Doppelhelix zeigt auf t ^ 3Usbilde " *3]. 

durch die Kippung der Basenn^™ che euie schmale und breite Furche, deren Tiefe 


de FCh ^ ^ PPUng der Ba senpaare gegen Und breite Furche - deren Tiefe 

fn dTS* rchC N ' Und O-Atome mit Protein^ ^ lurchmesser bestimmt wild. Wahrend in 
Proteine Furche weni ger basenspezifisch en !T )eSOnderS gUt wech selwirken konnen, sind 

Begl d , e DN ~ U "~ * - - 
meUUnder ^gewandelt werden kSn^^ui 8 ® 0 Einflussen, generieri. 
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Abb.1.1 


DNA-Doppelhelix: als Strukturformel links und als Kalottenmodel rechts. Gezeigt ist die Struktur des 
palindromischen Oktamers GGTATACC, die a us Rontgenstreudaten bestimmt wurde. [VoVo 95]. 


Die DNA-Kondensation im Zellkem zu hohen Packungsdichten erfordert daruber hinaus eine 
Flexibilitat des Doppelstranges. Bisweilen kann eine Verdrehung der Basen mnerhalb eines 
Basenpaares so weit gehen, da6 die Strange an dieser Stelle bereits teilsepariert sind. 
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hfrlikaler Drehsinn 

Durchmesser \ nml 

Basenpaare pro 

A -DNA 

rechts 

ca. 2.4 
11 

B-DNA 

rechts 

ca. 2,4 

10 

Z-DNA ~~ 

links 

ca. 1.8 

12 (6 Dimere) 

neiiKaier yv l 1 -j — — 

helikale Windung pro 

0,33 

0,36 

0,60 (pro Dimer) 

pasenpaar i j — — 

helikale Ganghohe 
A nstieff AVindune [nml 

2.8 

3.4 

4.5 

helikaler Anstieg pro 
Basenpaar [ nm] 

0,26 

0,34 

0,37 

Basenneisungswinkel [°] 

20 

6 

7 

groBe Furche 

eng und tief 

breit und tief 

flach 

kleine Furche 

breit und flach 

eng und tief 

eng und tief 

Ribose-Konformation 

C 3 -endo 

C 2 -endo 

Pyrimidin: C r -endo 
Purine: C^-endo 

glykosidische Bindung 

anti 

anti 

Pyrimidine: anti 
Purine: syn 

in Lsg. begiinstigt durch 

K + , Na + , Cs + 

geringe 

Salzkonzentration 

hohe 

Salzkonzentration 


Tab.1.1: 


Kenndaten der B-DNA sowie der A-und Z-DNA Konformationen [Eb 93]. 


1.2. Nukleosomen 

Bas^einheit^ fi^Chrom ^? er ^ chie der Kompaktierung die erste Stufe dar und bilden die 

1111(1 FelSenfeld zei S ten 1971, daB es m Chromatin DNA- 

Komberg fand daraufhin 1973 geschiitzt sind [KoKu 81]. 

,3eads on a String“ nannte (Abb tlUvZ Aufnahmen Stuktur, die er 

daB die Nukleosomen, die sie zunachst v Bod'l ' ^ ?° naW 0Uns zeigten dann 1974, 
[0101 74]. St V ® oc * ies nannten, Bestandteil dieser Struktur waren 

Bei Nukleosomen werden nun etwa nnn., 

oktameren Proteinkomplex gewunden nip “Pf 3 * 5 m 1 3/4 Windungen linkshandig um einen 
bezeichnet wird, verbindet die Nukleocnm hinaus § ehende DNA, die als Linker-DNA 

Basenpaare lang. Das Oktamer besteht ' en miteuiander und 1st je nach Organismus bis zu 80 
der Klassen H2A, H2B, H3 un d H4 au s zwe i Proteinen, den sogenannten Histonen, 

Uber die relative Lage der Linker-DNA w »• 

entsduadet ein weiteres Protein, Histon Hi rv^ Und Austritt a us der Wicklung (Abb.1.2) 
a er er einer Haamadel ahneln, uber die dpi- 7 ^ Un ^' :ion von Hi 1st zwar umstritten, konnte 
ornberg [KK 81] fand mit Hilfe einer „ Uga f g ^derer Proteine zur DNA steuerbar Lst. 
ontgen eugung, daB das Histonoktamer einp 0 ™ * nat * on von Elektronenmikroskopie und 
5-5 nm Hohe besitzt. Die Bindung der DNA ^d^* 16 Gestalt mit 1 1 nm Durchmesser 

8 der DNA erfolgt uber elektrostatische Wechsel- 
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Abb.1.2: 

Darstellung der 
Nukleosomenstruktur mit 
DNA, in der Aufsicht Oinks) 
und in der Seitenansicht 
(rechts), in einer Wire-Frame 
Darstellung. [Lug 97]. 


^iSL™^x e ?3r A ' phospaMckBr, ‘ m,d 120 p ° si,iven u ' )u " sen von l>s “" und 

Neueste rontgenWslallographische Aufnahmen von Lnge, [Log 97] mi. einer Anfliisung von 
.8 A zeigen, daB die DNA hier nur erne Ganghohe von 10,2bp besitzt. Dies konnte eine 
Erklarung fur bevorzugte Bmdungspositionen der DNA an das Oktamer sein. Sie fand 
weiterhin, daB die aminoterminalen Gruppen der Histonschwanze eine sehr geordnete Struktur 
bilden. Dies sollte entscheidenden EinfluB auf die Inter-Nukleosomen-Wechselwirkungen 
ausuben, und uber Acetylierung auch auf biologische Regulationsmechanismen der 
Erbsubstanz. 


13. Chromatin 

Chromatin stellt die zweite Stufe in der Kompaktierungshierarchie von DNA dar und setzt sich 
aus Nukleosomen zusammen. 

Isoliert wurde Chromatin, ebenfalls durch Miescher, bereits 1868. Er nannte es Nuklein. 



Links: EiektronenmikToskopische Aufnahme von Chromatin on a String" in SmM Salz fVoVo 95]. 

Rechts: Lie krone nmikroskoDiscbe Aufnahme der 30nm Chromatinfaser I \ o V o 95 !. 



Hen sich Nukleosomen wie Perlen zu einer Kette auf n e 
Wie oben erwahnt, reine Bereich kondensieren sie zu einer Faser mit ca Id* ! 

Saizkonzentralionen im phys ^ fllr Aufl>au dieser 30nm-Faser vorgeschlag^” 

“'^/ osomen einer Helbrwindung zusammengoiagert Die ? , 
^"ISii EStSSL. dad gleiche KrSfre zwischcn den Nukleosomen Z 

D'rSuptrte\d l Mo*r»o" ll 2»'8 raf ” nd Frante [ZeFr 841 Schlagl ' “SehMBige 

Ansammlungen von 6 bis 48 Nukleosomen vor. wobe. die m schonenden PrSpam,, orien 
ertahenen F^sem aus diesen sogennannten Superbeads bestehen sollen. 

EMemative dazu ware das von vielen Autoren [vHZ 95, WoHo 95; siehe dort wei terc ] aus 
experimentellen und theoretischen Daten vorgeschlagene 3D-Zick-Zack-Modell. Die 
Nukleosomen folgen hier beim Aufbau der 30nm-Faser einem zickzackformigen Weg und 
nicht einer Helix (Abb.1.3, links). 

Generell laBt sich Chromatin in zwei Klassen einteilen, die sich in der Packungsart 
unterscheiden: Heterochromatin ist inaktiv und besitzt keine oder nur wenige Gene 
Euchromatin ist dagegen aktiv oder in aktionsbereitem Zustand und genreich. 


1.4. Chromosomen 

Chromosomen (griech.: Farbkorper) als grofite strukturelle Einheit im Kern und die dritte 
Kompaktierungstufe von DNA wurden zunachst durch C. W. Nageli (1842) und W 
Hofmeister (1848) beobachteL Anfang der 70er Jahre des 19. Jahrhunderts wurden sie durch 
Farbemethoden als eigenstandige Struktur im Zellkem erkannt, erhielten jedoch erst von 
Waldeyer 1888 ihren Namen [Wa 88, Br 91]. Durch die damaligen FSrbemelhoden und 
“'f der Auflosung der Mikroskopieverfahren konnlen Chromosomen jedoch 

Werden ' 1411 ^“idhnet man das ZeUsladium. in 

(Abb.1.4) die in die heid T h rora ° somen stellen in der Mitose definierte Gebilde dar 
diskutierL ob Chromosnm ° ° C ^ trans P ort ^ ert werden. Schon damals wurde dariiber 
der Mitose neu eebildet werd^ it , egrenzte ^inheiten sind oder ob sie jedes Mai zu Begin 
™scben Zellteilungen) eine denmemSiXT^hito * ” ““ ‘ nt " PhaSe (d “ S'**™ 

Chrom^men^ai^auf^o^^^^^J^ 6 ? elle im ZdWtem den gleichen vollstiSndigen 
zwischen 50Mbp (<50^ 21^3^”“ 2 "■>* 22 Chromosomen. die 

homologen Paare bezeichnet man dabei a. (C hromosom *> enthalten (Tabm.l). Die 22 
Geschlechts-Chromosomen bezeichn^r iu • utosomen ’ wahrend man die iibrigen zwei als 
Chromosom, wahrend die weiblichen zwei Y.r^ SaugerzeUen enthalten ein X- und ein Y- 
Chromosomen, das auch oft als Barr-Kn u ^ ‘ lromosor nen besitzen. Eines der beiden X 
vollstandig inaktiv (mit Ausnahme der j” bezeichnet wird, ist in weiblichen Zellen fast 
letzel belegen, daB sich beide Chromosomen ° t Ut ° SOmalen Re 8 ion )- Untersuchungen von 

end der Metaphase, die ein Stadium der m Cf ^ altun 8 s struktur unterscheiden [Ku 96]. 
paarweise zusammen. Dabei stellt man <.i i.° Se re P r ^ sent iert, lagem sich die Autosomen 
rS C . hne1, Und Chrom °som in den kiira EmschnUrun g fest, die man als Zentromer 

“ d d “ la ”Sere„ ,-Ann te.h (Abb.14. 

Durch spezielle Farbetechniken TFr Qai , 

Querbanden-muster (Ideogram me; Abb.ll ^ !“*** C ^omosomen m der Metaphase 

m 5 ) sichtbar machen. Im allgemeineD 
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Abb.L.4: 

Links: Metaphaseplatte, Giemsa gefarbte Ideogrammbanden. [AI 94]. 

Rechts: Ekktronenmikroskopische Aufnahme eines Metaphasechromosoms. [BIFa]. 


wird dieses Muster auf den unterschiedlichen Adenin- und Thymin-Gehalt und den Zeitpunkt, 
in denen Gene in diesen Bereichen repliziert werden, zuriickgefiihrt. Als G-Banden werden 
dabei Banden bezeichnet, die hohen Adenin- und Thymin bzw. Gen-Gehalt aufweisen und friih 
replizieren. Mit R-Banden werden die restlichen Banden bezeichnet [Co 78], Kurzlich konnten 
diese Muster auch fiir die Interphase nachgewiesen werden [Zink unpublished], so daU also 
Teile der Chromosomenstruktur der Mitose auch in der Interphase erhalten sind. 

Bis heute ist jedoch die dreidimensionale Organisation von Chromosomen nicht bekannt. Die 
verschiedenen Modelle, die diese Faltung des Chromatins beschreiben, werden in Kapitel n 
naher erlautert. 


0 CUnmosoneW™*™ 



D Chromosomenorganisation im Interphasezellkern 


Die erste TTieorie der Zellkemstruktur postuUerte Rabl 1885 unter der Annahme, dafi in r ler 
und Pflanzenzellen eine struktureUe und genet, sche Identitat wahrend des Zell zyklu ' 
. , nm rum»rincnmpn Hip in der Metanhase distinct** 


und Pflanzenzellen eine struktureUe und geneuscne laentnat wahrend des Zell zyk T 
beibehalten wild [Bo 09]. Chromosomen, die in der Metaphase distinkte Objekte darstellen 
soUten auch wahrend der Interphase solche defmierten Temtonen einnehmen. 

Die mit dem Aulkommen der Elektronenmikroskopie verbundene Hoffnung, die einzel ne 
Chromosomen im Interphasezellkern auflosen zu konnen, erfiillte sich jedoch auch nich 
wahrend der zahlreichen Experimente in den 1960er und 1970er Jahren. Lediglich d' 
beobachtete Feinstruktur schien auf lange, sich im Zellkem bewegende Chromatinfasern 
hinzudeuten. 

Hieraus entstanden u. a. die ModeUe von Comings [Co 68], sowie Vogel und Schroeder IV 
74], die im Gegensatz zu Rabl von einer nichtterritorialen Organisation ausgingen Zu di 
Zeitpunkt war bereits allgemein akzeptiert, daB sich das Euchromatin vollstandig entpackt 
Chromosomen ihre gesamte dreidimcnsionale Organisation verlieren [St 77, Wi 73] Un 
Hauptsachlich Mikrobestrahlungsexperimente an lebenden chinesischen Hamster^,, 
sprachen jedoch in der zweiten Halfte der 70er und den friihen 80er Jahren fiir eine territoriale 
Orgamsauon der Chromosomen in der Interphase [Cr 74, Zo 79, Cr 82a, Cr 82b und Cr 841 
Mittels Fluoreszenz-in-situ-Hybridisierung (FISH) mit 7 iKamm?n 01 »c» t . , 

spezi&chen DNA Sonden gelang schUeBUch 1988 de S K^TcT”'"' 
somentemtorien im Interphasezellkern TCr 88 I i ss p; so f 8 INa ^P weis von Chromo- 

Als Model, wurde hieraTf do “h C rTmer [& 93) Id Lrto m ^ 

Domanen (ICD) ModeU vorgeschlaeen hw a ■ Z 93 J das Inter-Chromosora- 

DieCpneg^'T Eadun S der Chromosomenoberflachen^^osnihert^ V °” Zwb5chen “> die 

unbeklmnt^obirohl ma^ih^^ r°l C £ r0 ™ OSOmen blieb allerdings weiterhin 
Schon 1984 batten deshatt Pi en ? a T/rn be Fu ^n des Genoms zuschrieb. 
phasechromosomen vorgeschlaa^n \a . . Cy ein de taxlreicheres Modell fiir Meta- 
tinschleifen bestehen, die zu Rosetten P asecbrornosorn en sollten dabei aus Chroma- 
diese Rosetten aber wieder in einzelne meng f lagert sind. In der Interphase sollten sich 

Mu der Weiterentwickluno ' ScWeifen aufld sen. 


^ 

R^dn^w 1 !? Zeigen ’ da6 man GTn^Ah? tU H HybridiSierUng k0nnte van den Engh [vdE 
1997 "hi a ^’ Giant ' Loo P (RWGL) Modells ® de von DNA Sequenzen mit Hilfe eines 
1991 chlug Miinkel das Multi-L 00 n £f kann. 

d^Mod "i P aten ^sererklane ££ 1 J^**"™** ModeI1 (MLS) vor, das die experi- 
C VOn Pienta “"d Coffey, sowie deslfwri*^ 6116 ' Es stellt eine Weiterentwicklung 

“e ModeUe weeden „ ^ ^ 

en ^°*8enden Abschniden genauer erlaulerl. 
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Abb.D.1: 

A: Rabls Modell: 

Der Ort und die Orientierung jedes Interphase- 
chromosoms ist durch den Verlauf eines primaren 
Kernfadens (rechts) gegeben. Von diesem gehen 
sekundare und tertiare Faden als late rale Projek- 
tionen aus (links), die das dreidimensionale Chro- 
matinnetzwerk im Interphasekern bilden. 

B: Comings Modell: 

Die Befestigungspunkte der Chromatinfaser sind im 
Euchromatin weit, im Heterochromatin (H) eng 
zusammen. Die Chromosomen sind polarisiert, d. h. 
daB die Zentromere (C) sich auf einer Seite des 
Zellkerns und die Telomere sich auf der anderen 
befinden (B, A). 

[Cr 93]. 


II.l Modell von Rabl (1885) und Boveri (1909): 

Aufgrund von lichtmikroskopischen Untersuchungen an Zellen von Salaraandra maculata und 
Proteus Zellen veroffentlichte Carl Rabl 1885 seine Theorie iiber die interne Struktur des 
Interphase-Zellkems. Chromosomen behalten in dieser Theorie ihre Anaphase-Telophase- 
Orientierung wahrend der folgenden Interphase bei (Abb.n.l, A). Nach diesem Modell gehen 
von einem primaren Kemfaden dabei sekundare und tertiare Faden als laterale Projektionen 
aus, wobei sie das dreidimensionale Chromatinnetzwerk des Interphasezellkems bilden. Die 
primaren Kemfaden gehen von einem sogenannten Polfeld aus, das auf der einen Seite des 
Zellkerns lokalisiert ist und in dem sich die Zentromere bevorzugt aufhalten. Die Kemfaden 
laufen bevorzugt entlang der Kemhiille und enden auf der Gegenpolseite. Die Telomere der 
Chromosomen enden an den Kemhiillen. 

Theodori Boveri erweiterte dieses Modell dahingegen, indem er spekulierte, daB die 
Chromosomen bestimmte Territorien einnehmen [Ra 85, Bo 09]. 

IL2 Modell von Coinings (1968): 

Durch die Entwicklung der Elektronenmikroskopie und die Moglichkeit, nun mit hoher 
Auflosung die Topologie des Zellkerns zu erforschen, gelangten Comings [Co 68] und andere 
allerdings zu da von abweichenden Modellen (Abb.ILl, B): 

Die Chromosomen sind so an der Kemhiille, den Nucleoli und an einer intemen Kemmatrix 
befestigt, Haft sie relativ genau definierte Positionen im Zellkem wahrend der Interphase 
einnehmen. Entzieht man durch mehrmaliges Waschen in NaCl (2M) dem Zellkem emen 
GroBteil der DNA und Histone, so bleibt eine sogenannte Kemmatrix iibrig, an der bei 
unvollstandiger Preparation der groBte Teil (abgesehen von der Kemmembran) des Chromatins 
befestigt ist Die Lokalisation von Euchromatin im zentralen Bereich nahe der Kemmatrix und 
von Heterochromatin im Bereich der Kemhiille unterstutzt diese These. Ein weiterer Hinweis 
isl daB man fruhreoliziere nde DNA nahe der Kemmatrix. SDatreDlizierende daeeeen un 
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0-84 ^ 


Bp Dichte 
lObp 1 

80bp 6-7 
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60kbp 680 


IMbp 12000 


18 12000 
Schlei- 
fen 



1 M. IVllUl 


»r^<«^, A p ° r |3‘'Zl,!!!L S n b r > "‘ Ch ”“" l! ' ““'"P*"* P™ ' V1 » d ""8 

mit 1908 Schleifen in 106 Minibandern a \*\ m h “ man€n Meta Ph^echron,osoms 4, ohne Zentromer, 
Chroma tin-Faser. Das Modell kfi/c c e en zu 60kbp. Die weiBe Faser entspricht der 30nm 

Modell 1 st 26,5cm lang und aquivalent zu einem Chromosom von 3,19pm Lange. 

peripheren Bereich fand Comine h#* 

Polarisation der ^^rphase-Chromosomeli^ 16 ^ V ° n aufgestellte Annahme liber die 


VOnP, enta und Coffey ( 1984 ): 

Ausgehend von dem Modell vn r* 

SSg" StUdle " 1984 tPiCo PiCnta Und C ° ffey aufgrund e ‘ 

VogeU^ visualisierte 1980 DNa .< . ““ 

bromid ( y 0 " ^“P^elikalen die 80 die Kernrnatrix gebunden 

Halo urn di e p eri . ° 8 °^' Die Vis ualisieruno h”" 1 Cmer niedr igen Konzentration Ethi< 
einer DNA-Schleif ^ der Kernm atrix Aus d P 1 ^ er DNA erzeugt dabei einen sogena 
&hle,fe 01 ca - 90 000 Basenp^e Halos wurde der DNA-^ 
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Sollten diese Schleifen an einera Basispunkt an der Kemmatrix befestigt sein, so sollte sich bei 
hohen Konzmrationon v °n Ethidiumbromid (lOOug/ml) die urspni^gliche Struktur wiedS 
ausbilden [BeWo 76], Vogelstem konnte dies ebenfalls 1980 zeigen [Vo 80]. BesaB die DNA 

SST 3 " 8 ’ trat der Effekt nicht ein ’ was die Hypothese end ^ - 

Wahrend Paulson und Laemmli [PaLa 77] bereits ahnliche Schleifen fanden, wurde die 
Bedeutung der Kemmatnx durch Experimente von Pardoll [Pa 80] gesteigert, der zeigen 
konnte, daB in diesen Halos Replikation stattfindet. Es wardabei auch moglich, zeitaufgeloste 
Studien mit Hilfe der Autoradiographie gewinnen. 6 

Die erste einfache Uberlegung von Pienta und Coffey bestand nun darin, die DNA-Schleifen 
nur wahrend der Replikation als freie DNA-Schleifen anzusehen und sie fur den ubrigen 
Zeitraum, sicher jedoch fiir die Metaphase, als kondensiertes Chromatin zu betrachten Durch 
Mittelung liber aUe verfugbaren veroffentlichten Daten (genauere Ubersicht in [PiCo 84] 
versuchten sie nun, allgemeine Daten fiir die Schleifenstruktur im Kern zu erhalten Diese 
wurden in die Metaphase hinein extrapoliert. Ein Vergleich mit bereits existierenden 
Vorstellungen ergab nur eine sinnvolle Ubereinstimmung mit dem Radial-Loop-Modell von 
Marsden und Laemmli [MaLa 79] bzw. Adolph und Kreisman [AdKr 83] fiir die Metaphase. 
Chromosom 4 besteht demzufolge aus ca. 106 Minibandem, die wiederum jeweils ca. 18 
Schleifen aus Chromatin enthalten, was die gemessene Metaphasenausdehnung sehr gut 
wiedergab. Mit elektronenmikroskopischen Aufnahmen, in denen man die Schleifenanzahl pro 
Chromatidwindung zahlte [Ut 81, Ad 80, La 79 und Du 70], konnte dies in etwa bestadgt 
werden. 

Zusammen mit weiteren Ergebnissen von Cook [Cook 76] schlossen Pienta und Coffey, daB, 
obwohl man in der Interphase keine Strukturen direkt sichtbar machen konnte, DNA-Schleifen 
erhalten bleiben. Die Minibander liegen hier wahrscheinlich in aufgeloster Form vor. Die 
Kondensation bzw. die Dekondensation zu Minibandem sollte auf chemomechanischen 
Wechselwirkungen beruhen. 


II.4 Inter-Chromosom-Domanen (ICD) Modell (1990-1996): 

Die territoriale Organisation der Chromosomen im Interphasezellkem konnte 1988 mit Hilfe 
neuer leistungsfahigerer konfokaler Laser-Scanning-Mikroskopie und Fluoreszenz-in-situ- 
Hybridisierung erstmals sichtbar gemacht werden [Cr 93, Zi 93]. Dies bestatigte die Hypothese 
von Rabl und Boveri. Es konnte jedoch nicht ausgeschlossen werden, daB aus den relativ 
kompakten Territorien weit ausgedehnte, mit aktiven Genen besetzte Chromatinschleifen 
herausragen [KoCo 88]. 

Neuere Daten von R. Eils [Ei 95] unter Verwendung von dreidimensionalen Voronoi- 
Berechnungen von FISH-markierten Chromosomen zeigen, daB die Chromosomenterritorien 
sich nur in einem recht kleinen Randbereich iiberlappen. Dies bestatigen auch 
Untersuchungen, bei der die Lage FISH-markierter Gene relativ zum Territorium [Cr 93, Ku 
96, Di 97], bzw. die Lage der Chromosomenarme und Chromosomenbanden bestimmt wurde. 
Die Entstehung dieser Territorien kann jedoch bis heute noch nicht erklart werden. 
Erklarungsmodelle, die eine groBe Kemmatrix verwenden (siehe z. B. Comings), besitzen den 
Nachteil, daB die Komponenten (Strukturproteine) der Matrix bisher nicht gefunden werden 
konnten [Dr 95, Mat 97]. Maniotis [Ma 97] beschrieb jedoch kiirzlich Experimente mit 
Mikromanipulatoren und Videomikroskopie, die auf ein kraftiibertragendes Netzwerk im 
Zellkem und im Cytoplasma schlieBen lassen, das, wie die Autoren glauben, aus DNA besteht. 


i interphasezellkern 




Fibrillen 


ICD-Space 


groBe, mikroskopisch 
sichtbare Aggregate 


Abb.II.3: 
Schematische 
Darstellung des ICD- 
Modells: 

Zwischen den Terri - 
torien befindet sich 
der ICS-Raum, in dem 
der Transport von 
Proteinen, Replikation 
und Transkription 
stattfinden soli. 

[Cr 93, Zi 93] 


Das ICD-Modell (Abb.n.3) beruht nun auf der Vorstellung, daB durch kurzreichweitige 
elektrostatische Krafte, sogenannte Donnan-Potentiale, zwischen den Chromoso- 
menoberflachen ein dreidimensionaler Zwischenraum entsteht, das den Kem als Netzwerk 
durchzieht. 

In diesem Netzwerk sollen neben dem m-RNA-Transport zu den Kemporen solch wichtige 
Prozesse wie DNA-Replikation, Transkription, Splicing und Reparaturprozesse stattfinden. 
Die Oberflachen der Chromosomenterritorien erlangen damit groBe funktionelle Bedeutung. 


II.5 Random- W alk-Giant-Loop-(RW GL) Modell (1992-1997) 

Cremcr und Zirbel konnten 1988 zeigen, daB Chromosomen der visuellen Beobachtung 
zufolge offenbar voneinander getrennte Interphaseterritorien einnehmen. Uber die 
sne^rr^ T 1 emt0nen * uf g en omischen Langenskalen von 5Mbp bis 300Mbp und 
Pienta und Coffe h onnten J ed °ch nur Vermutungen angestellt werden. Das Modell von 
jetzt nicht nachweisbar smT SC Clfenstrukturen S ef ordert, welche aber in der Interphase bis 

Sequenzenln Meroh^Chromoso" 1 m latenUen Pro J ekti onen Genom-Abstande von DNA- 

menschlichen Fibroblasten-Zellen iT^CO/cf Ph^ h ^' ^ MeSSUngCn WUrde " h “ 
Klonabstande variierten zwischen 0 ful? PhaSC dur chgefuhrt. Die genomischen 

Praparationsbedingungen ergab sich dahc .' 1 ^ ,p und 4Mb P- Unter den gewahlten 

zufalligen Weg folgte. Der eenm^t i! ei ^ ^ rgan ^ sation » bei der die Chromatinfaser einem 
einem Chromosom war deshalb nrn 6 • ^ Stand zwischen zwei markierten Punkten auf 
(Abb.n.4). p oportl °nal zur Quadratwurzel des Genomabstandes 

Fur Genomabstande, die groBer als 1 smk 

Anwachsen des gemessenen geometnsrh^ *if waren > er gab sich jedoch ein langsameres 

unabhangig (Tr 91], Wahrend Hahnfeldi rw oder er wurde schlieBlich sogar davon 

l a yjj diese Abweichungen durch die Einbettung 
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Daten von Yokota [Yo 95]. Offene Kreise: Hypertonlsche Preparation. Voile Kreise: isotonische 
Preparation. 

A: Quadrat des mittleren raumlichen Abstandes fur kleine genomische Abstande. 

B: Quadrat des mittleren raumlichen Abstandes fur groBe genomische Abstande. 

der Interphasechromosomen in kugelformige Doraanen erklarte, schlug Ostashevsky [Os 94] 
vor, daB die Chromatinfaser in mehrere Mbp groBe Schleifen gefaltet ist. Abstandsmessungen 
von Yokota 1995 [Yo 95] fur Genomabstande von 0.15Mbp bis 190Mbp auf den 
Chromosomen 4, 5 und 19 ftihiten zu der Folgerung von Sachs [Sa 95], daB es wahrscheinlich 
nur zwei Ebenen der Chromatinorganisation auf Skalen bis 200Mbp gibt. Als einfachste 
Interpretation ergibt sich eine Organisation von flexiblen Chromatinschleifen von einer GroBe 
zwischen IMbp und 5Mbp, deren Befestigungspunkte entlang eines Random-Walk-Backbone 
plaziert sind, der aus Proteinen, ahnlich denen der Kemmatrix bestehen soli,. 

Diese Vorstellung ist aufgrund der mittleren Ausdehnung des Backbones (G1.IL13, siehe 
Unten) mit der Ausbildung von Territorien vereinbar. Dies laBt aber eine systematische 
Kompartimentalisierung von genreichen und genarmen Unterregionen im Territorium nicht zu. 
Korrespondierende Schleifen konnen daher mit gleicher Wahrscheinlichkeit im Inneren wie an 
der Peripherie des Territoriums gefunden werden. Dies steht zu den Ergebnissen des ICD- 
M odells in klarem Widerspruch. Wahrscheinlich ist die verwendete Auflosung der bisher 
verwendeten FISH und Laser-Scanning-Mikroskopie-Verfahren nicht groB genug, oder die 
Feinstruktur der Chromosomenterritorien wird durch die Praparation beeinfluBt Die beiden 
postulierten Ebenen der Chromatinstruktur sollen im weiteren nun naher beschrieben werden. 


II.5.1 Random- Walk-Modell vod 0.1 bis 1.5Mbp 

Das von van den Engh beobachtete Verhalten entspricht dem Abstandsverhalten eines 
PolymermodeUs, das zur Beschreibung von groBen, flexiblen. Unear aufgebauten Molekiilen 
verwendet wird. Spezifische lokale Strukturen auf kleine Skalen werden durch ein gemitteltes 
Verhalten auf groBem Skalenbereich angenahert. Die Chromatinfaser wird als Folge von 
Segmenten angenahert, die in gleichen Genomabstanden lokalisiert sind. 

Die GrbBc des statistischen Segmentes ffeier DNA betragt ca. 1000 Basenpaare, wahrend dies 
bei der 30nm Chromatinfaser ca. lOkbp bis lOOkbp ist [Ha 93, Os 94], Deshalb wird fur die 
Segmente ein Wert von lOOkbp als sinnvoU angenommen (Abb.II.5, Unks). 

Fiir die Ableitung des Zusammenhangs der Projektion des geometrischen Abstandes mit dem 
genomischen Abstand werden die Segmentpositionen in kartesischen Koordmaten X,, l - 0, 1, 
... . N beschrieben. Die Z-Koordinate stimmt mit der Strahlrichtung im Mikroskop uberein. 


II Chromosomenorganisation im Lntcrphasez 



Abb.IL5: Mathematisches Modell des RWGL-Modells 

Links: Bereich von 0,1 Mbp bis l,5Mbp. Die Segmente (0 bis N) dieses Modelbsmd entbng der 
Chroraatinfaser in gleichen Genomabstiinden lokalisiert. Die Segmente 0 und N beze.chnene d.e Telomere. 
Die Abbildung zeigt die Projektion dieser Segmente auf die XY-Ebene. 

Rechts- Bereich von 2Mbp bb 200Mbp. Schleifen mit J Segmenten sind an den Punkten 0, J, 2J, etc. nut 
dem Backbone verbunden. Die Segmente i, j bezekhnen mit FISH markierte DNA-Sequenzen. 


Die folgenden Rechnungen betreffen somit nur die (X,Y)-Ebene. Die Komponenten konnen 
wie folgt beschrieben werden: 


x i =(x,,Y i ,z i )=(x„z 1 .). (n.1) 

Man kann die statistische Segmentposition von Polymeren nun als Zustande eines 
mechanischen Systems betrachten. Fiir eine bestimmte Segmentanordung im thermischen 
Gleichgewicht ergibt sich folgende Wahrscheinlichkeitsdichte in Form einer 
B oltzmann verteilung : 


-u 

G(X 0 ,...,X w ) = Ce*' r . 


(n. 2 ) 


C: Normierungskons Xante 
u(x 0 ,... t X N ): potentielle Systemenergie 
k B : Boltzmannfaktor 
T: Systemtemperatur 


GemaB dem Genomabstand fiir zwei markierte Genomsequenzen wird nun ein entsprechendes 
Segmentpaar i,j gewahlt und die Projektion des geometrischen Abstandes r zu 

*,-i/(x,-x J )’*(r,-r j ) ! =| i( - i .| w.3) 


besummt. 
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D.e W^rschemUchkeitsdichte P(r) fiir den Abstand r erhalt man, indem man zunachst die 

t ^ C f 61 * C te 2 f*** zwe ' kartesische Intervalle u,v durch Integration fiber alle 
Intervalle, auBer u und v, von G aus GLIL2 berechnet: 

P *(*» ~*>) s P j(“ ~ v ) = p A X i ~ Xj.Y, - Yj), (IL4) 

p t(*' -Xj) = Cj dsdw dZ t dZ j | dZ ( ^ 


Hierbei ist: 




(n.6) 


Die potentielle Energie U = U 0 mit 


U o = 




(H.7) 


beschreibt die Wechselwirkung zwischen benachbarten Segmenten, wobei die nominate 
Federkonstante Kq liber alle Beitrage der Chromatinstruktur in kleineren Skalenbereichen 
summiert. | i-j | ist hierbei der projizierte Genomabstand der Segmente i und j, die durch n 
Segmentlangen voneinander getrennt sind. 

Somit ergibt sich aus GLEL5 mit H.7 die GauBform: 



Mit Hilfe der fiir zwei Dimensionen angenaherten Jakobigleichung 

lx 

P(r) = rdr $ ddP 2 {r cos d,r sind), 

0 


(n.8) 


(n.9) 


erbalt man nun die Rayleigh-WahrscheinUchkeitsdichte P(r) fur Genomabstande von O.lMbp 

bis 1.5 Mbp: 



Der mittlere Abstand <r> betragt dann: 
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<r> = J rP(r¥ r = J(*[) 


(II- 1 1) 


bzw. 


(r 2 ) = J r 2 P(r )dr = 2I\ 


(H.12) 


so daB also das Quadrat des projizierten, mittleren Abstandes <r > zweier Segmente 
Die P ^lgeradne Beriehung fur einen uneingeschrankten Random Walk lautet fiir n-»~ 


(r 2 ) = b 2 n, 


(H-i3) 


wobei n die Anzahl der Segmente der Lange b bezeichnet [Doi 88]. 


n.5.2 Random- Walk-Giant-Loop-(RWGL) Modell von 2 bis 200Mbp: 

Wiirde die Chromatinstruktur eines Chromosoms iiber den gesamte^Bereich von O.lMbp bis 
200 Mbp hinweg einem Random Walk folgen, wiirde der Abstand <r > nach GLEL13 auch im 
gesamten Bereich linear vom Genomabstand n abhangen. Yokotas MeBergebnisse lassen sich 
zwar mit einer Geraden der Steigung 2.15 (im /Mbp im Bereich von O.lMbp bis 0.15Mbp 
annahem, im Bereich groBer als 2Mbp betragt die Steigung jedoch nur noch 0.08 1 pm 2 /Mbp 
(Abb.II.4). Sachs folgerte daraus, daB es noch zusatzliche Wechselwirkungen zwischen 
Segmenten, die weiter auseinanderhegen, geben muB (Abb.H.6). 

Der Wechsel der Geradensteigung entspricht dem Ubergang von einem Zufallsweg (Random 
Walk) bei kleineren Skalen zu einem dichteren bei groBeren Skalen, der sich am besten durch 
ein Modell mit groBen Schleifen (Giant Loop) beschreiben laBt Das Modell besteht nun der 
achheit halber aus Schleifen gleicher Lange (3,13Mbp), die entlang eines Random Walk 
lokalisiert sind (Abb.n.5, rechts). Die Segmente des Random-Walk-Weges, welche durch 0, J, 
R ’ f .' l n • * -» N - J, N beschrieben werden, sind gleichzeitig die 
0 mThSII ^" e - ,hr Absl ” d sich mit der experimentellen Regression 

m“i« ™a P <^ 01D l 3 “ CA “»"">• “e Segmente i und j beaeichne. mi. FISH 
beiden Segmente rni^R^rin formale Herleitung des physikalischen Abstandes dieser 
SeX^ “ der POtentieUen Ener ^ ie U au6er U o noch ein weiterer 




(11.14) 


hinzugefugt werden. 

Wenn N/J die Gesamtzahl der Schi^if 

Wechselwirkung benachbarter Befectionrf 11 e ^ c ^ re *t> 1 . dann bezeichnet dieser Term die 
Segment 0 und J, bzw. J und 2J bis (N ^ punJcle ’ al5 ° die Wechselwirkung zwischen dem 
zwischen benachbarten Segmenten und k J N ' U ° bt wieder der Wechselwirkungstei® 
Ertnh des Abstandes aw de„T Ko ““'“ e«Be, 0 fit, den Random Walt F * 

hen awe, ntadtietten Segment ij rand wiede, & 
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Wahrscheinlichkeitsdichte P(r) mit G1.II.14 und einer etwas anderen Konstante V berechnet 
werden, die die in GLII.12 gezeigte Rayleighform besitzt. 


II.5.3 Uneindeutigkeit des Random-Walk-Giant-Loop-Modells : 

1997 zeigte Bing Liu [LiSa 97] die Nichteindeutigkeit des Sachsschen Modells. Die 
MeBergebnisse von Yokota [Yo 95] lassen sich auch durch ein Modell mit zumindest zwei 
Random- Walk-Backbones beschreiben. Hierfiir muB nun ein zweiter Erganzungsterm in der 
Gleichung fur die potentielle Systemenergie eingefiihrt werden, wobei nun die 
Befestigungspunkte J/2, 3J/2, ... , 2N/J-1 fiir Schleifen am Backbone hinzugefugt werden: 


Mit denselben Argumenten, die schon oben verwendet wurden, erhalt man nun fiir die 
Wahrscheinlichkeitsdichte P(r) wiederum eine Rayleigh Verteilung. Diese fiihrt auf einen 
mehrere Seiten langen Ausdruck fiir die <r>-Abhangikeit in Bezug auf den genomischen 
Abstand d, weshalb hier auf die Orginalpublikation verwiesen sei [LiSa 97]. 

Eine nunmehr nur noch numerische Analyse ergibt folgendes Bild fiir die drei Parameter k 1? k 2 
und die LoopgroBe: Fiir genomische Abstande unterhalb von IMbp ergibt sich mit k x = 0.46 
<r > = 0.025 +1.88d in guter Ubereinstimmung mit den experimentellen^Daten. Fiir groBe 
genomische Abstande ergibt sich mit einer SchleifengroBe von 6Mbp ein <r > = 2.37 + 0.08d, 
wobei k 2 = 1.0, also der schon oben erwahnte Wert. Weitere Verallgemeinerungen unter 
emeuter Einflihrung eines Oder mehrerer weiterer Backbones sind somit plausibel. 


II.6. Multi-Loop-Subcompartment- (MLS) Modell 

Miinkel entwickelte 1996 das Multi-Loop-Subkompartment-ModeU, da die bisherigen Modelle 
die Gesamtheit der experimentellen Daten nur unzureichend beschreibt [Mii 98]. 
Chromatinschleifen mit einer GroBe von ca. lOOkbp werden zu Rosetten, den sogenannten 
Subkompartimenten zusammengeschlossen. Die Rosetten werden mit DNA, der sogenannten 
Linker-DNA, verbunden. Somit entfallt die Notwendigkeit eines Backbones (abb.n.6). 

Indem man den Linker zu einer Schleife zusammenzieht oder offnet, ergibt sich ein natiirlicher 
Ubergang von der Meta- in die Interphase. Hierzu konnte ein Protein zum Beispiel eine 
Verbindung zwischen den Nukleosomen herstellen. Damit ergibt sich eine einfache Kontrolle 
der Funktion des Genoms durch die Veranderung der dreidimensionalen Organisation. 


U = U 0 + U X + U 2 = 
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Sell rofie: 5Mbp 


Backbone 
(Nicht DNA) 



R-Bande 


Befestigungs 

Punkte 


G-Bande 



30nm dick 


RosettengroBe l-2Mbp 
(gemaB der GroBe der 
Interphase-Ideogramme) 

Linker besteht aus DNA 
(Im Gegenssatz zum Backbone) 


Abb.n.6: 

Links: Das RWGL-Modeli. Rechts: Das Multi-Loop*Subcompartment-ModeU. 

Die GroBe der Rosetten entspricht einer Ideogrammbande in der Interphase. Die Ideogramme 
der Interphase unterscheiden sich von denen der Metaphase durch ihre Bandenzahl. Beim 
Ubergang von der Metaphase in die Interphase wachst sie von 850 auf ca. 3000. Allerdings 
verschwindet das Bandenmuster bei ublichen Farbemethoden normalerweise in der Interphase. 
Kurzhch konnte jedoch [Zink, unpublished] nachweisen, daB es auch dort existiert. 

Fakten d * e Erzeu S un g von Bandenmustem auf organisatorische 

vermutet wird [Co 78] ^ asenzusammensetzun g zuruckgefiihrt werden, wie schon lange 

Pr^Lfr: ^ok^as^ Experim^t g h emeSSenen Daten Yokotas teilweise mit 
hypertonisch behandelte Zellkeme LLnKetKp -T* Pra P aratlonsmethode " ; 
Zellkeme mit Paraformaldehyd fixiert wuZ a ^m S6SSlg flXiert ’ w °g e S en isotonische 
sich dabei fur die isotonische Pr o nar af 1 , Cn ; Auf klemen genomischen Abstanden ergeben 
Die hypertonische Praparahon Ken z “™ RWGL-ModeU(Abb.U.4). 

Methanol Eisessig zerstort die Chroma. . , Cm Stark anschwe llen, und die Fixierung mit 
Histone [Co 78, ReRu 77] Es bestpht^^ 1 ^ erheblich durch die radikale Beseitigung der 
stiitzenden Ergebnisse als Preparation dlC Verrautun g> daB die das RWGL-Modell 

spricht bei dieser Preparation von nicht h ^^ ang6Sehen werden konnen. Yokota selbst 
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UI Simulationen 


Im simulatorischen Tei] dieser Arbeit wurden Chromosotnenterritorien als Polymerkette 
versehiedenarlig modelherl. Aus diesen Simulationen warden expenmentell enLSS? 

°r c rr ^ iisr ,,efen “ •- * ***— 

Neben raumlichen Abstanden fllr Vergleiche mit intragenomischen Abstanden wurden neue 
Observablen gefunden und getestet. Mit ihnen laBt sich in den Simulat.onen zwtschen den 

Modellen unterscheiden. Daruberhinaus konnen Experiments mit den erhaltenen Vorhersagen 
mterpretiert werden. 5 

to H _ m ^ ck auf zukunfuge komplexere FragesteUungen der dreidimensionalen Organisation 
des Zellkemes envies sich dabei erne Analyse mit Hilfe von Fraktalen als besonders 
ohnenswert So konnte erstmals gezeigt werden, daB die Chromatinfaser multifraktales 
Verhalten aufweist. 


to folgenden wird zunachst auf die verwendeten Methoden der Modellierung eingegangen 
werden. Danach werden einige theoretische Aspekte der Analyse komplexer fraktaler Systeme 
erlautert. An die VorsteUung der Observablen, die zur Simulationsanalyse verwendet werden 
schlieBt eine Darlegung der Simulationsdurchfuhrung an. SchlieBlich werden die 
Simulationsergebnisse prasentiert und es erfolgt eine Zusammenfassung des simulatorischen 
Teils dieser Arbeit. 


III.l. Simulationsmethoden 

In dieser Arbeit wurden Monte-Carlo- und Brownsche Dynamik-Verfahren zur Simulation 
benutzt Brownsche-Dynamik-Verfahren berechnen die Trajektorie eines Systems in Raum 
und Zeit unter EinfluB thermischer Fluktuationen. Daraus konnen dynamische GroBen, wie 
z.B. Newtonsche Krafte, Beschleunigungen und Geschwindigkeiten berechnet werden. 
Resultierende Observablen wie freie Energie Oder Entropie konnen ebenfalls erhalten werden. 
Monte-Carlo- Verfahren dagegen dienen dazu, den Phasenraum reprasentativ abzutasten, ohne 
zeitliche Informationen zu liefem. GroBe Systeme wie Chromatin oder ganze Chromosomen 
miissen dabei aus groBeren geometrischen Untereinheiten zusammengesetzt werden, da eine 
Simulation auf atomarer Ebene iiber groBe Zeitskalen die Kapazitaten heutiger Rechner bei 
weitem ubersteigt. 


MI-1.1: Chromatin als Polymer - AUgemeines 

Die Chromatinfaser wird in dieser Arbeit vereinfacht als fadenformiges Polymer dargestellt 
(Kapitel 11.5). Die atomaren Details der Chromatinfaser gehen dabei in Form elastischer 
Eigenschaften ein. 

Ein Polymer kann naherungsweise durch eine Kette aus steifen Segmenten beschrieben 
werden, die durch isotrop bewegliche Gelenke miteinander verbunden sind. Fur die 
Wahrscheinlichkeit des End-zu-End-Abstandes einer solchen „GauBschen Kette“ gilt: 



(m.l) 
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„ ... „ Hip pin Ma 8 fur die Flexibilitat der Kette darstellt: je steifer das 

b ist dabei die Segmentlange, die ein Mau 

Polymer, umso gr°B er b - . p ers istenzlange L P , die die Steifigkeit einer sogenannten 

Fine verwandte GroBe list die Pers ^ ^ die mitdere Projektion des En d-zu-End- 

IZSf to Antogsrictong «Jet Keue Oder to Marion eweie, Rich.ungsvetoren 

der Kette wie in [GrKh 94]: 


(«£>*), =( cos (*(*))) 


(III.2) 


oc e 


III. 1.2: Chromatin als Polymer - die modellierten Eigenschaften 

Die Chromatinfaser wurde in dieser Arbeit als lineare Kette von Chromatinsegmenten 
(Abb.III.l) bei einer absoluten Temperatur von 310 K modelliert [Mu 98]. Als Segmentlange 
1 q wurde die Kuhnsche-Lange gewahlt, die sich fiir die GauBsche Kette wie folgt ergibt. 


<«=)= 



(III.3) 


<R 2 > ist dabei der mittlere quadratische End-zu-End-Abstand einer Chromatinkette. Mit Hilfe 
der von van den Engh [vdE 92] bestimmten Abhangigkeit des raumlichen Abstandes vom 
genomischen Abstand ergibt sich mit einer Chromatindichte von d * 105bp/nm, eine Kuhnsche 
Lange von 260 nm (Abs. II. 5. 2). Diese zeigt gute Ubereinstimmung mit friiheren Ab- 
schatzungen, die sich zwischen 200 nm und 359 nm bewegen [Os 94]. Fiir die Simulationen 
wurde ein mittlerer Wert von 300 nm gewahlt, der ca. 31000bp entspricht. 


R(i+1) 



Abb.III.l: 

Segmentanordnung, ihre 
Freiheitsgrade und 
Potentiale. 

Zwei Segmente sind iiber 
ein Gelenk miteinander 
verbunden. Die Koordi- 
naten der Gelenke 
werden mit R(i) be- 
schrieben. Die Segment- 
lange L ist variabel. Die 
rote Linie ist die 
Resultierende eines 
spharischen Excluded 
Volume Potentials um 
die Gelenke. 
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Die Segmente konnen in ihrer Lange variieren, so daB ihnen die folgende Streckenergie 
zugewiesen wurde: 


U sV) = *f{l-lo)\ (m4) 

wobei 1 die aktuelle und 1 0 die mittlere Lange darstellt. Da die Strecksteifigkeit fur Chromatin 
nicht bekannt ist, wurde aus numerischen Stabilitatsgriinden 5 = 0.1 gewahlt 
Fiir die Biegesteifigkeit zwischen zwei benachbarten Segmenten wurde ebenfalls ein einfaches 
harmonisches Potential gewahlt. 


U B (P) = 



(ffl.5) 


wobei P der Winkel zwischen diesen Segmenten ist, und \|f mit der Kuhnschen Lange gemaB 

, _A_ an.6) 

Lk ~ 2 y/ 2 


verkniipft ist [Mii 98]. 

Eine Torsionssteifigkeit wurde in der Modellierung nicht beriicksichtigt, da dies eine 
erhebliche Reduktion der Rechenzeit erlaubt. AuBerdem laBt sich diese Naherung durch die 
Prasenz torsionsstreB-reduzierender Proteine, wie Topoisomerase n, rechtfertigen. 
Chromosomen oder Chromatinfasem durchdringen sich nicht. Ura dies zu beriicksichtigen, 
wurde eine sogenannte Excluded-Volume-Wechselwirkung (EV-WW) eingefuhrt: 


U E y=U° Ev k l 


J , r A -2r 2 c r 2 ] 


(m.7) 


r ist dabei der Abstand zu einem Segment. In den Simulationen stellte es sich als sinnvoll 
heraus, r c auf das I.5fache der Dicke einer Chromatinfaser zu setzen. Die Hohe des zu 
iiberwindenden Potentials und damit die Wahrscheinlichkeit fiir das gegenseitige 
Durchdringen der Chromatinfasem kann durch Auswahl von U gy variiert werden. 

Die 46 Chromosomen verteilen sich auf den verfugbaren Raum im Zellkem. Fiir die 
Simulation eines einzelnen Chromosoms muB daher eine spharische Energiebamere Inbound 
eingefiihrt werden. Diese hat eine ahnliche Form wie Ugy und ist genauso hoch wie bei der 
EV-WW. Das Volumen ist proportional zum Verhaltnis aus DNA-Gehalt des Chromosomes 
c Chrom ^ dem DNA-Gehalt des gesamten humanen Genomes C genom , so daB man mit dem 
Kemradius R^n, in Ubereinstimmung mit der Literatur [HH 89] den Territorienradius 




'err ^Kem 


{ 


■'Ckrom 


(ffl-8) 


setzt Eine einfache Abschatzung fiir den End-zu-End-Abstand einer freien Chromatinkette 
ergibt fiir Chromosom 15 mit einer Lange von 106Mbp (Tab.ni.l) sofort 


\l\un. 


(m.9) 


U _ l2°2™x 1,06 10 s bp « 

d y 105-^- 


nm 



einen Wert von 27 \xm. Dies iiberschreitet bei weitem die GroBe 


Schleitenbiiaung iad s. u.j uu ~ r 

ahnlich dem Streckpotential aus Gl.m.4 zusammengehalten: 



(DUO) 


III.13: Brownsche-Dynamik-Algorithmen 

Die Theorie der Brownschen Bewegung wurde fiir Systeme von Teilchen entwickelt, deren 
Masse und raumliche Ausdehnung um ein Vielf aches groBer sind als die der Molekiile des - 
normalerweise fluiden - Einbettungsmediums. Diese Bedingung ist auch fiir Chromosomen im 
Zellkem erfiillt. Die Simulationen wurden gemaB dem Ermak-McCammon-Algorithmus 
durchgefiihrt [ErMc 78]. 

Die Langevin-Gleichungen beschreiben ein System von N Teilchen in einem solchen viskosen 
Einbettungs medium. Thermische Fluktuationen des Mediums iiben eine zufallige Kraft f { auf 
das System aus. Dadurch tritt in der klassischen „newtonschen“ Bewegungsgleichung neben 
der Summe Fj iiber die inneren und auBeren Krafte und einem Reibungstensor fiir 
geschwindigkeitsabhangige Krafte noch ein Term fiir die zufalligen Krafte f { auf: 


wobei rMie Masse, v.die Geschwindigkeit und i,j den Index der Koordinaten (1 < i, j <3N) der 
1 euchen bezeichnen. J 

Fur klassische Brownsche Systeme ist gauBverteilt mit Kovarianz = 0 und Mittelwert = 0: 
(/.(*)) = 0 und (fi(t)fi(t')) = 20 ^ 3(1 (m.12) 


w ^=-Z^ + ^ + Za,/., 


(m.ii) 


D* ist der Diffusionskoeffizient von i Fiir die ^ , 

dem Boltzmannfaktor k, der Zusammenhang K ffl2ienten “ij 8 Ut mit der Temperatur T und 


^ IJ k ji, ( 111 . 13 ) 


(m. 13) 



(m.i4) 
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In enter Nahemng kann die Hydrodynamik einer Polymerkette durch Kugeln gleicher GroBe 
(RacLus c 0 heschrieben werden, die an den Verbindungspunk ten der Segment 

lokahsiert srnd [ChLa 94] Der Diffusionstensor laBt sich dann als Rotne-Prager-Tensor mit 
der Viskositat T| darstellen [RoPr69]: 5 mK 


D = 5 , — (5 -l) 
6nrjo ' ? ' 


k„T 


rl 


V 'N 

1+^2- 

2a 2 

H 

I Vvl 

Ll 4 J 

4 

l 3 ij. 


(m.15) 


Nach [ErMc 78] relaxiert die Impulsverteilung viel schneUer zur Gleichgewichtsverteilung als 
die Positionsverteilung. Fur groBe Zeiten erhalt man dann fUr die Zeitentwicklung der 
Koordinaten den Ausdmck 


r f (/ + St) = r,(r) + £ D ij F j ~ + R,(t) (DL16) 

;=l 

und fur die wiederum gauBverteilte Verschiebung Rj(t) den Mittelwert und die Korrelation: 

(3(0) = 0 und (R l (t)R i {t')) = 2D y St. OH. 17) 

Mit einer gegebenen Anfangskonfiguration laBt sich mit G1.HI.15 und Gl.m.16 eine neue 
Konfiguration im nachsten Zeitschritt berechnen. Die dazugehorigen Krafte erhalt man fur 
unser Problem durch Differenzierung der entsprechenden Potentiale aus Abschnitt HI. 1.2. 


III. 1. 4: Monte-Carl o-Algorithmen 

Bei Monte-Carlo- Verfahren wird ein mathematisches Problem naherungsweise dadurch gelost, 
daB man einzelne reprasentative Werte unter bestimmten Randbedingungen berechnet. 
Urspriinglich wurden Monte-Carlo-Verfahren am 2 Ende des Zweiten Weltkrieges durch von 
Neumann, Ulam und Metropolis entwickelt [MR 2 T 2 53]. 

Die reprasentative Auswahl wird hierbei durch eine Wahrscheinlichkeitsverteilung p erzeugt. 
Da nicht bekannt ist, wie p fiir einen gegebenen Phasenraum aussieht, muB der Phasenraum 
durchsucht werden. 

Hierzu erzeugt man eine sogenannte Markov-Kette von Zustanden, die einer kiinstlichen 
Phasenraumtrajektorie entspricht. Um die Markov-Kette zu erhalten, muB das System von 
seinem Istzustand in einen neuen Zustand iiberfiihrt werden. 

Zwei Zustande m,n sind durch eine Ubergangswahrscheinlichkeit 71^ verbunden, so daB aus 
einem Anfangszustand pj ein neuer Zustand mit P 2 = Pi 71 ergibt. hn Limes fur unendlich viele 
Ubergange erhalt man die Eigenwertgleichung p7t = p. 

Man zeigen kann zeigen, daB 


5>--‘ <ml8) 


isl bezeichnet man als stochastische Matrix. 

Um erne Markov-Kette und damit eine Phasenraumtrajektorie zu erzeugen, wahlt man nach 

(MRV 53] 
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n = cc m 

'*mn m 


Pn ~ Pm 


(H-19) 





Pn<Pm 


fur m^n und a als eine symmetrische stochastische Matrix. 

Um von einem Zustand in den nachsten zu gelangen, fiihrt man das System mit einer 
zufalligen Verschiebung, einem sogenannten Monte-Carlo-Schritt, in einen anderen Zustand 
fiber. Daraus resultiert eine Energieanderung 8V mn . Um das Verhaltnis p„/p m und damit die 
Ubergangswahrscheinlichkeit zu bestimmen, bildet man den Boltzmannfaktor der 
Energiedifferenz [AlTi 87]: 



e (-py.) e (-fw.) 




(IL20) 


1st < 0, so wird der Monte-Carlo-Schritt immer akzeptiert, ist dagegen > 0, so 
wird die Veranderung jedoch nur mit einer bestimmten Wahrscheinlichkeit akzeptiert. 

Die Kunst der Metropolis-Monte-Carlo-Methode besteht darin, die Schritte fur ein System so 

zu wahlen, daB die Akzeptanzrate moglichst hoch ist, um den Phasenraum moglichst schnell 
abzutasten. 

Um die Chromatinfaser zu modellieren, wurden nun folgende Schritte verwendet: 


1 ) 

2 ) 


Zwei Segmente werden um einen zufalligen Winkel um die Verbindungslinie des 
Anfangs des emen und des Endes des anderen Segmentes gedreht. 

B !! k vp°rH a h? S n le i fen c (AbS ' 11 ' 5 md a6) werden um einen zufalligen Winkel um ihre 

(siehe Basispo,e "“ 

bezefcbtrtm^ ak^S^hritL Um ^ Verbindun g slinie seiner Enden gedreht. Dies 

besteht aus einem lsl em sogenannter Relaxierungsschritt notwendig. Er 

‘ aus emem Brownschen-Dynamik-Schrit. ohne ZufaUskrafi. 

Nachdem ein Schritt durcheefiihn miM. ■ , . 

die Energie, die Energiedifferenz 5V mU der Potentiale aus Abschnitt HI. 1-2 

Die SchrittgroBe wurde so ecwahU°Tn Afczeptanz die s er neuen Trajektorie bestimmL 

Akzeptanzrate unterschritten so werden h- -“f Akzeptan z ra te bei 50 ± 5 % liegt. Wird die 

’ Crden dlC Schr *ttgro6en verkleinert, und umgekehrL 


3) 

4) 
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III.2: Fraktale Analyseverfahren komplexer Systeme 


III. 2.1: Einfiihrung 

r 1C I he t nC d t F Fra . ktale ™ urde von Benoit Mandelbrot in den 70er und 80er Jahren zur 
Eesehreibung komplexer Systeme entwickelt [Man83]. Fraktale sind selbstahnl.che' odcr 
se bstaffm^ skaleninvarmnte Objekte, die sowohl als regulare Fraktale TexaSe 

au““ aUCh 3,5 ZU,allS ' ° dCr <«*«**» sSSnSS 

Em Beispiel fur ein Fraktal stellt eine Kustenlinie dar, die immer prinzipiell ahnlich aussieht, 
wenn man ste aus verschiedener Hohe, also mit verschiedener Auflosung, betrachtet. Fur die 
Lange der Kiistenlinie als Funktion der Auflosung ergibt sich dabei ein Skalengesetz 
Fiir die Unterscheidung verschiedener Strukturen erhiilt dabei man sofort ein Mali in Form dcs 
Exponenten des Skalengesetzes. Weiterhin erhiilt man die Moglichkeit, verschiedenc Klassen 
von Skalengesetzen miteinander zu vergleichen. 

Insbesonders gelten solche Skalengesetze auch fiir die raumliche Anordung von Objektcn und 
die in einer solchen Anordnung ablaufenden dynamischen Prozesse. Die Theorie der Fraktale 
stellt hierdurch ein Mali fur die Irregularitat und die Raumerfullung eines Systems dar 
(Abb.III.2). Die Irregularitat der Kiistenlinie laBt sich mit der fraktalen Dimension D (dem 
Exponenten eines Skalengesetzes) charakterisieren. 

Fiir die Lange der Kustenlinie mit verschieden groBen MaBstaben (Auflosung) erhalt man 


N(h)och~ D bzw. L(A) = A^(A)- A oc A 1 ’ 


(III.21) 


N(A) ist dabei die Zahl der MaBstabe der Lange A, mit der die Lange L der Kurve bestimmt 


r'% 







Topologische 

Fraktale 


Dimension 

Dimension 


1 

1,00 


1 

1.02 

Abb.DI.2: 

1 

1,25 

Hier wird gezeigt, wie sich die 
fraktale Dimension einer Kurve 
bei immer starkerer Faltung 
verandert, bzw. den hier 
zweidimensionalen Raum der 

1 

1,45 

Flache immer mehr ausfuilt. 
[Ka 89J 


1 Selbstahnliche Objekte = Vereinigungsmenge isotrop-reskalierter Kopien ihrer selbst 
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nachsten in gerader Linie gemessen. 


euklidischer Dimension E: 


D up <D Hml =D<E. (m.22) 


(m.22) 


Mit der so erzeugten fraktalen Dimension konnen nun Aussagen iiber die Entstehung des 
untersuchten Systems oder die Prozessdynamik, insbesondere aber Diffusionseigenschaften 
gemacht werden. 


m.2.2: Das Stechzirkel und Box-Verfahren 

Die fraktale Dimension einer Kurve wird aus den verschiedenen Kurvenlangen L in 
Abhangigkeit von der Lange der MaBstabe bestimmt. Hierzu wird im allgemeinen die 
Stechzirkelmethode verwendet (Abb.n.3): 

Mit einer Zirkelspanne einer bestimmten GroBe beginnt man an einem beliebigen 
Einstichspunkt und sucht den Schnittpunkt des Zirkelkreises mit der Kurve. Dieser gilt als 
nachster Einstichspunkt. Der Zirkel wird entweder von auBen oder von innen an die Kurve 
herangeffihrt, da in realen Datensatzen die genauen Kurvenkoordinaten meist unbekannt sind. 
Fehler bei der Vemachlassigung von Einbuchtungen konnen durch Mittelwertbildung fiber 
verschiedene Anfangspunkte verkleinert, aber nicht beliebig minimiert werden. Dieses 
Verfahren wiederholt man fur verschiedene Zirkelspannen, um Gl.m.21 anwenden zu konnen. 
Mit dem Box-Verfahren bestimmt man die fraktale Dimension einer Massenverteilung. Der 
Raum wird dabei mit Wfirfeln einer Kantenlange fiberdeckt und das vom Objekt belegte 
Volumen bestimmt. Variation der Kantenlange ffihrt ahnlich zu Gl.m.21 auf 


N boM) x ^ D bzw. V' fca (A) = N(A) • A 3 oc A 3 ' d (HI. 26) 



Einschwingverfahren 



Ausschwungverfahren 


Abb.HI.3: 

Das Stechzirkelverfahren: 
Gezeigt wird sowohl die 
Methode des Ein- 
schwingens, als auch des 
Ausschwingens. 
Einbuchtungen der Kurve 
konnen leicht vernachlassigt 

werden. 
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IIL3. Observablen der Simulationsanalyse 


Fur die Analyse der Simulation wurden GroBen berechnet, die es zulieBen, zwischen den 
simulierten Modellen der Chromosomenorganisation zu unterscheiden. Die Observablen 
dienen zur Interpretation und zur Vorhersage von experimentellen Daten Dies gUt 
insbesondere fur intragenomische Abstande. Neben diesen wurden weitere Observable 

“ ber ^ » “on 

In den folgenden Abschnitten werden die verwendeten Observablen kurz erlautert. Zuvor soli 
aber noch auf das geschnebene Analyseprogramm Geometry eingegangen und allgemeine 
Aussagen uber die Observablen gemacht werden: 


Analyseprogramm Geometry: 

Fur die Analyse der Simulationen wurde ein Program m namens „Geometry“, das die 
Observablen berechnet, in der objektorientierten Programmiersprache C++ entwickelt. 
Objektorientiert bedeutet, daB das Programm in separierten Einheiten, den sogenannten 
Objekten, organisiert ist. Die einzelnen Objekte konnen dabei mit anderen Objekten nur uber 
genau definierte Schnittstellen in Verbindung treten. Die detaillierten Vorgange innerhalb 
eines Objektes sind somit fur andere Objekte irrelevant. Objektorientierte Programmicrung bot 
die Moglichkeit, Anderungen schnell und effizient zu implementieren, da diese nur in den 
einzelnen Objekten vorgenommen werden miissen. Anderungen im gesamten Programmablauf 
wurden dabei selten notig. 


Einordnung der Observablen in einzelne Klassen: 

Die verwendeten Observablen lassen sich generell in solche mit globaler, halb-globaler und 
lokaler Aussagekraft trennen: 

• Globale Observablen machen eine Aussage iiber das Gesamtsystem, wobei iiber die 
gesamte lokale Struktur des Systems gemittelt wird. 

- Halb-globale Observable mitteln iiber Teilbereiche des gesamten Systems. Sie machen 
somit Aussagen iiber lokale Strukturen, betrachten aber deren Verhalten im gesamten 
System. 

- Lokale Observable machen dagegen eine Aussage iiber eine genau spezifizierte lokale 
Struktur eines Systems. 

Lokale und halb-globale Observable konnen dabei durch globale Eigenschaften des Systems 
beeinfluBt werden. 


1113. 1 Energies, Geometriedaten und Segmentlange (global) 

Fur die Energien der einzelnen Potentiale, den Gyrationsradius, die Eigenwerte des 
Tragheits tensors und die Scgmentlange wurden die Haufigkeitsverteilungen besummt 
Alle GroBen machen dabei eine Aussage uber das Relaxierungsverhaiten e'ner S'mulaUon von 
der Anfangskonfiguraiion ins Gleichgewicht. Ist die Standar. ia b w e .chung der Seg^ 
menilangenvene.lung zum Beisp.el klein, so ist die Simulation gut relax jert und 
umgekehn.Der Gyrauonsradius und die Eigenwerte steUen dariiberhinaus noch ein gules MaB 
lur die Aubdehnung eines Chromosoms dar. 
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•i „„„ crhnell zu berechnen sind, dienten sie zur Online 

Da die Haufigkeitsverteilungen schneu zu 
Uberwachung dcr Simulationen. 


II, JU Radiate Aufenthaltswahrschein.ichkeit und Dichte (global) 

Die radiale Aufenthaltswahrscheinlichkeit wurde als die Haufigkeitsverteilung der Masse in 
Abhangigkeit des radialen Abstandes zum Massenschwerpunk berechnet. Die radiale Dichte 
isi die auf das Volumen einer Kugelschale normierte, Aufendialts-wahrscheinlichkeU in einem 
bestimmten radialen Abstand zum Massenschwerpunkt. Die Breite der Kugelschale wurde 

durch die Auflosung der Haufigkeitsverteilung vorgegeben. 

Fur die dreidimensionale Organisation von Chromosomen spielt die Massenverteilung eine 
wesentliche RoUe. Sowohl bevorzugte Regionen, als auch ihr Randverhalten konnen bestimmt 
werden. Eine schnell abfallende Dichte spricht bei Chromosomen fur eine geschlossene 
Struktur mit geringer Wechselwirkungsflache. Diese Angabe ist jedoch relativ zur maximalen 
Dichte und macht eine Aussage iiber die ..Rauhigkcit der Chromosomenoberflache. 


III33 Raumliche Abstande 

FromBegin - Abstand zu einem festen Punkt (global): 

Der FromBegin-Abstand bestimmt die raumliche Distanz als Funktion des genomischen 
Abstands, bezogen auf ein Ende der Chromatinfaser des Chromosomes. Dabei lauft man durch 
die Chromatinkette und bestimmt den Abstand zum Ende des p-Armes des Chromosoms.Der 
FromBegin-Abstand macht damit eine Aussage liber die globale Faltung der Chromatinfaser 
beziiglich eines festen Referenzpunktes. 


All Pairs und Random Pairs (halbglobal): 

Im Gegensatz zu FromBegin wurden auch die raumlichen Abstande aller Chromatinsegmente 
(All Pairs) zueinander als Funktion des genomischen Abstands, bestimmt. In einem System 
endlicher GroBe, wie es hier vorlag, ist dabei die Anzahl der paarweisen genomischen 
Abstande antiproportional zu den genomischen Abstanden selbst. Da die Abstande nicht mehr 
bezuglich ernes festen Referenzpunktes bestimmt werden, ergibt sich eine Mittelung iiber 
verschiedene raumliche Strukturen mit gleichem DNA-Gehalt. 

Damn ergibt sich erne AbstandsverteUung, die einem variablen Schleifenbasispunkt im MLS- 
zu dfesem R^«n w ark ^ r P a ^. mit einem gegebenen genomischen Abstand wurde relativ 
Verteilune der ^Srhl f S ° Ig * okalisiert se *n.Dies fiihrte zwar zu Aussagen uber die 


UU.4 Raumliche Abaind, twidchc g.nomische. Mmk.ro (total| 

Raumliche Abstande zwischen 

eines festen Referenzpunktes bestimmt werH^ geilomischen Abstands konnen bezuglich 
wahlt den Referenzpunkt auf der Chrnm r r Wahlt dabei eine genomische Distanz und 
lokalisien ist. Die resulUerende Absiandsver.l’n^J* 115 ’ a " der einer der beiden Marker 
iiber die An der Basispunkte eine Aussaee m » -h 1 "^ Ver 8 ,e * c h m it dem Experiment 

Schleifen-basispunkie fest sind oder nichi ° ^ ns ^ >eson dere laBt sich feststellen, ob die 
Da nicht bekannt ist, wo ein reales M 

zwischen verschobenen Anfangsp^terseTr^^ 115 ^" ,St ’ mUssen raumliche Abstande 

S wer den. Es wurde festgestellt, daB 
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1113^ Aufenthaltswahrscheinlichkeit von genomischen Markern (global) 

Analog zu Abschnitt III.3.2 kann auch die radiale Markeraufenthaltswahrscheinlichkeit und die 
radiale Markerdichte bestimmt werden. Dies fiihrt zu Aussagen iiber dynamische Prozesse wie 
zum Beispiel Umlagerungen bestimmter Substrukturen an die Oberflache eines 
Chromosomenterritoriums. Eine nichtlokale Markeraufenthaltswahrscheinlichkeit deutet dabei 
auf eine variable und dynamische Organisation hin und umgekehrt. 

Im Vergleich mit Experimenten konnte somit eine Entscheidung iiber die Validitat eines 
Modells beziiglich zusatzlicher topologischer Einschrankungen (Kernmatrix) getroffen 
werden. Ebenfalls wird eine Abschatzung moglich, inwieweit Abstandsmessungen von 
Markem, relativ zur Chromosomenoberflache, mit realistischen Modellen beschreibar sind. 

III-3.6 Fraktale Dimension (komplex [global, halb-global, lokal]) 

In der Strukturbiologie wurden bisher fraktale Analysen nur bei Proteinen eingesetzt, um deren 
Peptid-Riickgrat und Oberflachenrauhigkeiten zu bestimmen. Als erster Ansatz wurde 
zunachst ein exaktes Stechzirkelverfahren entwickelt (Abb.m.4). Dieses Verfahren berechnet 
die fraktale Dimension der Chromatinfaser. 

Von einem zufallig ausgewahlten Segment wird die Kurvenlange der Chromatinfaser mit 
verschiedenen MaBlangen zu beiden Enden der Chromatinfaser hin bestimmt. Der Schnitt- 
punkt mit der Chromatinfaser fur das exakte Stechzirkelverfahren ergibt sich in zwei Schritten: 
Zunachst wird der Abstand zwischen jetzigem Einstichspunkt und nachstem Segmentanfang 
bestimmt. 1st dieser kleiner als das gewahlte MaB, wird dies fiir den folgenden Segmentanfang 
wiederholt, bis der Abstand groBer wird als die MaBlange. Im dem Segment, in welchem die 
MaBlange iiberschritten wird, muB also der nachste Einstichpunkt liegen. 


R 



LI < M < L2 


Exaktes Stechzirkelverfahren. 
Zunachst wird untersucht, in 
welchem Segment der nachste 
Einstichpunkt des Zirkels der 
Lange M liegt (LI < M < L2). 
Dann wird aus den Vektoren 
von R nach R+l und R+l nach 
R+2 und der Lange M der neue 
Einstichpunkt exakt berechnet. 


Abb.m.4: 
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. cppmentanfangs und des Zirkeleinstichpunktes laBt sich dann 
Mit den Koordmaten des Segm > ^ t berechnen (AbbJL4). 

algebraisch der neue Zirke ren die Kurvenlange sehr gut. Fur groBe MaBe mu Bte 

S.lung iiber mehret. Anfangspunkte erfolgen, da hie, die spezifische Stntktu, 
jedocti eine mmci e def M i tte lungen wurde bei verschieden groBen MaBen so 

ge e wSu,1S S dTe n ltandardabweichung iiber den ganzen MaBbereich vernachlassigbar wurde, 

Folgende Mittelungen wurden vorgenommen: 


Onm < MaB < 
50nm < MaB < 
lOOnm < MaB < 
500nm < MaB < 
lOOOnm < MaB < 


50nm => 
lOOnm => 
500nm => 
lOOOnm => 
L Cut-Off =* 


1 Anfangspunkt 

2 Anfangspunkte 
4 Anfangspunkte 
6 Anfangspunkte 
8 Anfangspunkte 


Als AnfangsmaB wahlt man normalerweise ein Lange von ungefahr 10% der kleinsten 
strukturaufbauenden Einheit. Im allgemeinen miBt man dabei also die fraktale Dimension eines 
euklidischen Objektes aus. Uberschreitet das MaB eine gewisse Lange, so gelangt man zu 
einem Endlichkeits-Cut-Off, da reale Objekte raumlich limitiert sind. MaBe, die mehr als 70% 
bis 80% der Gesamtausdehnung eines Objektes ausmachen, sind fur die Bestimmung der 
fraktalen Dimension deshalb nicht mehr sinnvoll [Ka 89]. 


IH.4. Simulationsdurchfiihrung 

In dieser Arbeit wurden unterschiedliche Modelle der dreidimensionalen Organisation von 
Chromosomenterritorien simuliert, um zwischen den Modellen zu unterscheiden und um 
Vorher^gen fur experimentell meBbare GroBen zu machen. Vor allem wurden dabei Kriterien 
fur die Unterscheidung von RWGL- und MLS-Modell entwickelt. 

exnJmp^'if n * n '^ rCr p 653011 * 16 * 1 zu untersuchen und realistische Aussagen fiber 
variiert werHp 6 aten zu a * ten ’ mu ^ ten bei den Simulationen verschiedene Parameter 
Organisation im 7 H 1 W ^ ^ ^ ussa S en fur zukiinftige komplexere Fragestellungen der 
2SZT.2SS; n~ beiei " flUS “ dabei Simulationsansatz, bzw^hAe 

der Linker-DNA N h,CI,CnSrdGe ’ d ' e Schleife narizahl und die Lange 

groBen Schleifen, besondere Srhlp f r ° 6 mul ^ ten ^ a bei, zur schnelleren Relaxierung der 

Dk S sichd^°" d T g innCrhalb ^"^sTtte SIT ™ ldaL Im MLS ' Mode11 wurde noch 

undeiner hohen (l.OkT) Excluded if J Vurde, i ^^^binaus noch mit einer niedrigen (O.IkT) 
Die Simulationen (EV-WW) simuliert. 

gewahlten Anfangszustand abhangen oder nicht ^ 1611 Werden ^ onnte > ob die Ergebnisse vom 
t-riekte, die bei der Chro 

W«Xturfcfglt Ch h'"'* 115 ™' Br0WnSCh “ Dynamik e ' ne 
“ Chromosomen in ihrem nativen Voiumen simuliert. um 
* den EinfluB von globalen auf lokale u 

ermeiden und abzuschatzen. ™ tUren zu unte rsuchen, und um Randeffekte zu 




El Simulationen 



Schleifenwechselwirkung in, Hinblick ™, 

dX^te^Sr 1 - - eszeggen, 

gewinnen ist. Dabci zeigt sich auch ob es'^ " C » 8 elChaare statistische Awssagekraft zu 
Kemverbandes zu shSm ' ,n “ lne Chro ">‘«»"»" autohalb das 

d M ^: Mode " aussssm 

sszr^s. 


rhrnmntnm 7^ T ^ S |™ ulationsdur chfuhrung fur Chromosom 15 genauer besprochen 
Chromosom 1 5 wurde gewahlt, um den Vergleich mit dem experimentellen Teil dieser Arbeit 
so genau wie moglich durchfiihren zu konnen. Arbeit 


m.4.1 Ideogramme 

Da eme Rosette im MLS-Modell einer Ideogrammbande entspricht und deren DNA-Gehalt 
esitzen sollte, muBte eine zuverlassige Liste fur die genomische GroBe der 850 ISCN(1985) 
Mctaphasebanden [ISCN 85] erstellt werden (Abb.ID.5). Hierzu wurden die hochauflbsenden 
Daten uber Chromosomemdeogramme von Uta Francke [Fr 94] verwendet. Die mit 
verschiedenen Farbetechniken erstellten Human-Ideogramme stellen bisher die genaueste 
Bandenlokalisation dar (Fehler < 4%). Zunachst wurden die gesamten Daten digitalisiert und 
die nichtenthaltenen ZentromergroBen mit einem Lineal aus den maBstabsgetreuen 
Abbildungen bestimmt. Die von U. Francke und anderen vermerkte unterschiedliche 
Schattierung der Banden wurde nicht berucksichtigt, da bis heute nicht bekannt ist, was hierfiir 
den eigentlichen Grund darstellt (Konformationsanderung, andere Basenkonzentration, 
Proteinbesatz oder eine Mischung aus allem [Co 78]). Unter der Hypothese, daB jede 
Interphasenbande einer Rosette entspricht, wurde es notwendig, die heteromorphen Regionen 
lql2 und 9ql2 zu halbieren, sowie Yql2 zu dritteln, da ihr DNA-Gehalt beziiglich dem 
durchschnittlicher Banden zu hoch war. In den Testsimulationen hatte dies zu irreal hohen 
Schleifenzahlen und Chromatindichten gefiihrt. Mit Hilfe von Tabelle II. 1 wurde es moglich, 
den DNA-Gehalt einer bestimmten Bande mit einer fur die Modellierung ausreichenden 
Genauigkeit von 5% zu bestimmen. 


IIL4.2 Erzeugen der Anfangskonfiguration 

Fiir das Erstellen der Anfangskonfiguration einer Simulation wurde ein weiteres 
objektorientiertes Programm namens „ChromCreate“ entwickelt. Im folgenden sollen die 
Moglichkeiten des Programms und die einstellbaren Parameter kurz erlautert werden. 

ChromCreate ermoglicht es, die Chromatinsegmente im Raum gemaB den Modell- 
vorstellungen zu positionieren. Es erzeugt dabei eine Konfigurationsdatei, die neben den 
Koordinaten auch eine Liste mit Attributen eines jeden Segments enthalt. Die Attribute geben 
die Lange des Segmentes und die Zugehorigkeit zu einer Schleife oder einem Linker an 
( Abb.III.6). Aufgrund dieser Angaben kann das von Munkel entwickelte 
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Simulationsprogramm entscheiden, welche Monte-Carlo oder Brownsche-Dynamik Schritte 
durchgefuhrt werden. 

Mit ChromCreate kann jedes gewunschte Chromosom oder Subkompartiment simuliert 
werden. Die Kenndaten des jeweiligen Chromosoms werden aus den digitalisierten 
Ideogrammdaten bestimnit. Somit ist es auch moglich, nur einige Ideogrammbanden zu 


ChromCreate wurde weiterhin so gestaltet, daB mit den unten angegebenen Parametern vollig 
beliebige Anordnungen von Schleifen und Linkem erzeugt werden konnen. 


Lange der Segmente: 

Die Segmentlange bedingt den Grad der Diskretisierung des Modells und damit seine 
Auflosung. Die Lange der Chromatinsegmente (in Nanometem) kann frei gewahlt werden. Ihr 
Basenpaargehalt wird dabei automatisch bestimmt. Fur Segmente, die sich im Linker befinden, 
und fur Segmente, die sich in einer Schleife befinden, konnen unterschiedliche Werte benutzt 
werden, so daB eine Feinabstimmung der Modelle erfolgen kann. 


Dicke der Chromatinfaser: 

Die Dicke der Chromatinfaser und damit die Wechselwirkung bei einer Annaherung zweier 
Segmente kann eingestellt werden. Fur Monte-Carlo-Simulationen hat sich ein Wert von 
60nm, fur Brownsche-Dynamik-Simulationen von 30nm (30nm-Faser) bewahrt. 

SchleifengroBe: 


Die SchleifengroBe wird normalerweise in Basenpaaren angeeeben. Im MLS-Modell wird mit 


simulieren. 



Schleifenanordnung: 



Schleifenform: 



-w'fu f ° lgen einem Rand om-Walk in der Ebene. Bei 
h _ ™ ^ re *d‘ mens * ona len Raum. Ab einem DNA- 
chleifen schneller ins Gleichgewicht, als 


zwei weitere: 
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Abb.III.5: 

Ideogramm- 
Darstellung des 
humanen Genomes 
wie es verwendet 
wurde um das 
MLS-Modell zu 
eichen [Vo Vo 95]. 
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Linker: 

Band: 

Lange: 


1111 
0 0 0 0 

50 50 50 50 


0 111 
N 0 0 0 

50 50 50 -1 


Abb.H.6: 

Attribute eines Segmentes j e 
nach seiner Position: 
Generell kennt j e d es 
Segment seine Lange. 
Linker besitzen den Linker- 
wert 1 und den Bandenwert 
0. Segmente einer Schleife 
besitzen fur beides eine 0, 
das erste Segment einer 
Schleife bestitzt IN (N = 
Schleifennummer), das 
letzte ON. Das allerletzte 
Segment eines Chromsomes 
besitzt eine Lange von -1. 


Lange des Linkers: 

Der Linker verbindet im MLS-Modell die Rosetten und im RWGL-Modell einzelne Schleifen 
raiteinander. Die Linkerlange kann in Segmenten oder prozentual zur IdeogrammbandengroBe 
gewahlt werden. Beim MLS-Modell wird der benotigte DNA-Gehalt aus der Ideogrammbande 
abgezweigt. Im RWGL-Modell dagegen verringert sich die Schleifenzahl eines ganzen 
Chromosoms. 


Form des Linkers: 

Die Form des Linkers entscheidet iiber die Ausdehnung und das Erscheinungsbild der 
Anfangskonfiguration. Im MLS-Modell bildet der Linker eine ebene Schleife zwischen den 
Rosetten. Diese Struktur entspricht in Form und Ausdehnung einem Metaphasechromosom. 

Die 2D- und 3D-Schleifenform dagegen macht es bei hoher EV-WW notwendig, den Linker 
durchSgen ^ 3D ' Random - Walk zu wahlen, damit die Schleifen sich nicht 


Implementierte Modelle: 

weiteren ak TAK Mndln ^ R y5’ L " Mode11 sind hier noch zwei weitere Modelle, die im 
bei denen die Schleifenirtft 6261 ^ 61 Werden ’ ““Plementiert: sie stellen RWGL-Modelle dar, 
Man nimmt dabei an, dffi dfe^n^n^cWdf 1 ^ ^ 255 °- Bandenmuster angepaBt wird. 
eine groBe Schleife ubergehen. f euier Rosette dur ch Praparationsartefakte m 


01.4 J Simulationsansatz des MLS-Modells 

MoctefZ Ctomram^™ u d “' n C e h r "“ S ,' m ' Jlat ' one " ™ n C. Miinkel bekaniu, daB ein MLS- 
20 Schleifen p„ Subhomp^,,, dle DaJ^Tytl^^^ U " d ‘° ^ 

Fiir die Simulation des MLS-Modell 

Schleifenanordnung, die DNA-Gehalt d * C ^ nker ^nge, die SchleifengroBe, die 

Wechselwirkung variiert. 8 samten Chromosoms und die Excluded Volume 
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Wire Frame BUd und Ray Traced Bild der Anfangskonfiguration von eines MLS-Modells. In rot und grtin 
sind abwechselnd die Rosetten angefarbt. SchleifengroBe 126kbp, Linkerlange 1200nm. 


Fiir eine sinnvolle Simulation wurden zwei Felder mit 16 beziehungsweise 6 Parametem 
angesetzt. 


Feld 1: 

Die Linkerlange wurde mit 600nm, 1200nm, 1800nm und 2400nm variiert. Fur jede 
Linkerlange wurden sodann die Schleifen in cis- und trans-Stellung angeordnet. Jede der acht 
Kombinationen wurde mit niedriger (0,lkT) und hoher (1,0 kT) EV-WW simuliert. Die 
Schleifenlange betrug 126kbp, mit der in Testsimulationen gute Ergebnisse erzielt worden war. 
Fur die Form der Schleifen und des Linkers wurden nur normale Schleifen verwendet Die 16 
sich ergebenden Simulationen sind im Anhang im Detail zusammengestellt (Tab.A.l). 


Feld 2: 

Fur die SchleifengroBe wurde 80kbp, lOOkbp, 158kbp und 251kbp gewahlt. Die 
SchleifengroBe von 126kbp kann aus Feld 1 entnommen werden. SchlieBlich wurde noch die 
ChromosomengroBe von 82Mbp, 90Mbp, iiber 106Mbp nach 95Mbp variiert, wobei der 
relative Volumenanteil des Chromosoms am Kem beibehalten wurde. Es wurde eine trans- 
Anordnung der Schleifen, eine Linkerlange von 1200nm, normale Schleifenform und 
Linkerform, sowie eine niedrige EV-WW angenommen. Die 6 sich ergebenden Simulationen 
sind im Anhang im Detail zusammengestellt (Tab.A.2). 


m.4.4 Simulationsansatz des RWGL-Modells 

Fiir die Simulation des RWGL-Modells wurde die SchleifengroBe, die Linkerlange, die 
Schleife- und Linkerform und die Excluded Volume Wechselwirkung variiert 
Als Schleifenlange wurden die folgenden Werte verwendet: 5Mbp, 4Mbp, 3Mbp, 2Mbp, 
IMbp, 0.5Mbp, 0.25Mbp und 125Mbp. Dadurch wurde fiir gro Be Schleifen eine ausreichende 
Zahl von Moglichkeiten fiir den Vergleich mit Experimenten und fiir kleine Schleifen ein 
AnschluB an das MLS-Modell gewahrleistet. 

Die verschiedenen SchleifengroBen wurden mit niedriger und hoher EV-WW simuliert. Da das 
RWGL-Modell eine andere dreidimensionale Organisation zeigt, konnten Ergebnisse aus 
Simulationen des MLS-Modells nicht ohne weiteres iibemommen werden. Wahrend nun fiir 
die SchleifengroBen 5Mbp, 4Mbp, 3Mbp und 2Mbp bei hoher EV-WW stets eine 2D-Schleife 
verwendet wurde, wurde bei niedriger Wechselwirkung stets eine 3D-Schleifenform 
verwendet. Die Schleifen wurden durch einen Randomlinker verkniipft. 
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Schlei- 

fen- 

lange 

[Mbp] 

Schleifen 
lange 
in ganzen 
Segmente 

[Mbp] 

Schlei- 

fen- 
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Schlei- 
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300nml 
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12 

19 
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26.36 
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10 

25 

960.00 

10.24 
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3.0 

3.016145 

35.14 

809.73 

7.29 

8 

33 

835.57 

7.76 

8 => 8 

2.0 

2.010763 

52.72 

661.08 

4.86 

(6) _ 

48 

669.82 

5.33 

- 5 =»(6T~ 

1.0 

1.005381 

105.43 

467.48 

2.43 

2 

94 

495.08 

2.72 

3 =» 4 

0.5 

0.502690 

210.87 

330.54 

1.21 

(2) 

187 

351.01 

1.37 

1 =>2 

0.25 

0.251345 

421.73 

233.74 

0.61 

(2) 

333 

263.04 

0.77 

-1 => (2) 

0.125 

0.125672 

843.46 

165.28 

0.30 

(2) 

561 

202.66 

0.46 

-L=>(2)_ 


Tab.III.2: 


Bestimmung der Linkerlangen fiir das RWGL-Modell. 


Um Abhangigkeiten von der Linkerform, bzw. von der Ausgangskonfiguration, zu bestimmen, 
wurden fiir 5Mbp- Schleifen 3 weitere Simulationen angesetzt: Fiir 2D-Schleifen wurde eine 
normale Linkerform mit niedriger und hoher EV-WW und fiir eine 3D-Schleife eine normale 
Linkerform mit niedriger EV-WW (damit sich die Schleifen entwirren konnen) gewahlt. 
Dariiberhinaus wurde eine Simulation doppelt angesetzt, um die Entwicklung von gleichen 
Anfangskonfigurationen zu untersuchen. Fur SchleifengroBen ab IMbp wurden nur noch 
normale Schleifenformen und normale Linkerformen erzeugt. Zum Nachweis sollte aber eine 
Simulation mit 3D-Schleifenform bei normaler Linkerform durchgefiihrt werden. Sie sollte 
auch als Validitatstest fiir die dritte zusatzliche Simulation mit der SchleifengroBe von 5Mbp 
dienen. 

Bei verschiedenen SchleifengroBen besitzt ein RWGL-Chromosom eine unterschiedliche Zahl 
c°u, S ^ hleifen ’ so daB Chromosom 15, mit einem DNA-Gehalt von insgesamt 106Mbp, 20 
Schleifen mit 5Mbp und 561 Schleifen mit 126kbp besitzt. Das globale Verhalten wird aber 

i U n? i dlC Link „ erIa J n ^ bestimmt und ist nur in engen Grenzen variabel (Abs.II.5). Die 
erlange muBte deshalb spezifisch angepaBt werden (siehe Tab.A.3): 

Snalte' Ti^nH rTn f ° Igt ’ k ° nnte mit Hilfe der Anzahl der Linker (Tab.ffl.2, 

quadratische AbstanH m ^ rlan S e bestimmt werden. Der dazu benotigte mittlere 

< r 2 > _ 4800+2nn n- q > WUrde aUS Testsimulationen bestimmt und besitzt einen Wert von 
Walk. Die Zahl der Seom eBm ^ nte ’ d ' e d * esen Linker aufbauen, folgen wieder einem Random- 
und wiederum GUI. 13 ** Segmentlange VOn 3 °° nm 

aufgespalten. Eine vTrgroBemi^XrTinkr^ 6 ™ RWGL ' Mode11 in Schleifen und Linker 
Schleifen. Da jedoch der DNA Gehlh H * somit a utomatisch den DNA-Anteil der 

Rechnung wiederholt werden da sich Hi P "° ch nicht einb ezogen ist, muB die 

Linker. c eifenzahl anderte und damit auch die Zahl der 
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Abb.ffl.8: 

Ray-Traced Bilder von Anfangskonflgurationen des 
RWGL-Modells. 

I. 0.2 5Mbp-2DSchleifen, Random-Linker, 20 
Schleifen 

r.o.: 5Mbp-3DSchleifen, Random- Linker, 20 
Schleifen 

l.u.: 5Mbp-3DSchleifen, normale Linkerschleife, 20 
Schleifen 

r.u.:126kbp-Schleifen, normale Linker- und 
Schleifenform, 561 Schleifen 


Durch die Diskretisierung auf 300nm-Segmente, erfolgte fur 126kbp-Schleifen eine zu groBe 
Linkerlange mit einer Abweichung von ca. 30%. Diese konnte aber zur Bestimmung einer 
Observablenabhangigkeit von der Linkerlange im RWGL-Modell benutzt werden (analog zur 
Simulation des MLS-Modells), obwohl ChromCreate differenziertere Linker zulaBt. 

Die so erhaltenen 21 Simuiationen werden im Uberblick in Tab.A.3 des Anhangs gezeigt. 
Abb.n.8 zeigt Visualisierungen der Anfangskonflgurationen des RWGL-Modells. 


1II.4.5 Simulationsansatz anderer Modelle 

AbschlieBend wurden noch die TAK-Modelle und drei lineare Ketten simuliert. Die linearen 
Ketten, die bei niedriger und hoher EV-WW simuliert wurden, dienten dabei als Vergleich und 
zu Fehlersuche in der Simulationssoftware. Sie unterscheiden sich in den Attributen der 
Segmente und damit in den ausgefiihrten Monte-Carlo-Moves. In Tab.A.4 sind wiederum die 
Simulationsansatze niedergelegt. 
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Abb.m.9: 




Oben: Darstellung des Gyrationsradius (griin) und der Eigenwerte des Tragheitstensors (nach absteigen- 
der GroBe: blau, rot , rosa) fur das MLS-Modell bei links niedriger und rechts hoher EV-WW. 


m.4.6 Simulationsdurchftihrung: 

Insgesamt wurden 49 Simulationsansatze erzeugt. Da eine Simulation mit hoher Auflosung, 
d.h. kleiner Segmentlange, in Anbetracht der Rechenzeit nicht mbglich war, wurden die 
Simulationen in zwei Schritte aufgespalten. 

Im ersten Schritt wurden Monte-Carlo-Simulationen bei einer Segmentlange von 300nm 
durchgefuhrt. In einem zweiten Schritt wurde die Segmentlange dann auf 50nm gesenkt und 
mit Brownscher Dynamik nachrelaxiert, so daB man zu einer hohen Auflosung gelangte. 


Fur die Simdationsdurchfiihrung wurde das von Miinkel entwickelte Programm „VirtNucSim“ 
verwendet. Die Simulationen wurden dabei auf der IBM-SP2 mit 80 Knoten im Deutschen 
^ebsforschungszentrum und auf der IBM-SP2 mit 512 Knoten der Universitat Karlsruhe 
Brownth. n 236 MilUonen Monte-Carlo Schritte und 4,9 Millionen 

Z SS tei '" S8eS “ , ‘ b '" 8 “ gtt " d " Sto ““°" 5 - 5 CPU - 


Monte-Carlo-Simulationen: 

statistisch aus^^el^^ge^^K^^on^Monte Carlo 1 V ° n 3 °° Hm WUrde ^ 

Konfiguration, jede lOOste Energie und jeder lOOste Gyratio^r^ 
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dazugehorigen Tragheitstensors wurden abgespeichert. Aus dem Gyrationsradius und den 
Eigenwerten wurde mit Hilfe der jeweiligen Haufigkeitsverteilungen entschieden wann 
genugend Konfigurationen vorhanden waren. Dies beruhte auf Untersuchungen in denend.e 
Werte des Gyrationsradius als Funktion der Monte-Carlo-Schritte gegeneinander verschoben 
und mil -dem Statistikprogramm JMP [SAS96] eine Korrelationsanalyse durchgefuhrt wurde 
D.eses Verfahren wurde st.chprobenhaft vorgenommen. Die daraus gewonnenen Ergebnisse 
ergaben mit zusatzhch erzeugten Visualisierungen der Konfigurationen, daBjede 2000ste bis 
5000ste Configuration statistisch unabhangig war. Zwischen 200.000 und 500.000 Schritte 
wurden benotigt, um genugend Konfigurationen fur eine Auswertung zu besitzen (Tab A 1 
Tab.A.2, Tab.A.3 und Tab.A.6). 


Brownsche-Dynamik-Nachrelaxation: 

Nachdem die Monte-Carlo-Simulationen genugend viele unabhangige Konfigurationen erzeugt 
hatten, wurde die Aufldsung der Modelle erhoht. Die Segmentlange wurde dabei durch die 
Teilung der Segmente von 300nm auf 50nm (- 5300bp) abgesenkt. AnschlieBend wurden 1000 
Brownsche-Dynamik-Schritte durchgefuhrt, um eventuelle nichtphysikalische Artefakte der 
Monte-Carlo-Dynamik zu beseitigen. Pro Simulation wurden hierzu 100 Konfigurationen 
ausgewahlt, d. h. jede 2000ste bis 5000ste Monte-Carlo-Konfiguration. Die statistische 
Unabhangigkeit konnte durch eine Analyse der Energie, des Gyrationsradius und durch 
Visualisierungen wie oben beschrieben gewonnen werden. Fur das Herauskopieren und die 
Segmentlangenteilung wurden hierfur die Programme RelaxCopyConfig und 
RelaxKuhnToBending entwickelt. 

RelaxCopyConfig ergzeugt dabei auch die fur die SP2 notwendigen Process of Execution Files 
(POE). RelaxKuhnToBending verlangt neben der Eingabe der neuen Segmentzahl zusatzlich 
auch die Segmentlange und des DNA-Gehalts fur Segmente in Schleifen und Linker. Durch 
geschickte Wahl konnen dabei die Segmente gedehnt oder gestaucht werden, d. h. die 
Schleifen relaxieren zu etwas groBeren oder kleineren Langen. Hiermit wird eine Verfeinerung 
der Modellparameter erreicht, die gegeniiber einer volligen Neuberechnung mit Monte-Carlo- 
Schritten ca. 40% CPU-Zeit einspart. 

Die 4900 Nachrelaxationen bilden eine breite und liber verschiedene Modelle gestreute 
Grundlage fur die Untersuchung der Chromosomendynamik im Interphasezellkern, allerdings 
ohne die wiinschenswerte Dekondensierung aus der Metaphase in die Interphase. 


III.5 Ergebnisse der Simulationsanalyse 

Eine Gesamtanalyse mit Geometry liefert 51 Kurven pro Simulation. Observable, die lokale 
Aussagen machen, sind nur fur die nachrelaxierten Simulationen (hohe Auflosung) sinnvoll. 
Bis zur Abgabe dieser Arbeit konnten zwar die Nachrelaxationen vollstandig durchgefuhrt 
werden, eine Analyse war jedoch nur fiir die Monte-Carlo-Simulationen moglich. 

Somit wurden hier nur die Observablen FromBegin, RandomPairs, radiale Aufenthalts- 
wahrscheinlichkeit und Dichte, die Segmentlangenverteilung und die fraktale Dimension fiir 
die Analyse verwendet. 

In den folgenden Kapiteln werden die Ergebnisse der einzelnen Observablen detailliert 
dargestellt. Eine Ubersicht liber die gefundenen Parameterabhangigkeiten wird dabei in der 
Zusammenfassung der Simulationen (Abs.VI.l) gegeben werden. Die Ergebnisse der 
Simulation „anderer“ Modelle wird in Abschnitt III. 5. 4 erfolgen. 


Ill Simuiationen 


42 



Abb.ni.10: 

rr^r ener SimUU * ti0nen nach 200000 Carlo Schritten. Banden sind abwechsebd 

O.L.: MLS-Modell mit niedriger EV-WW .. , 

Subkompartimente zu sehen. ’ satz bch nachrelaxiert. Deutlich sind die distinkten 

O.R.: MLS-Modell bei hoher EV-WW in 2f h v 

StE"-* S “ P ' Stk '* if “ ^2S2. - bilden k „« d . s „., k „ 

ll-R.: RWGL-Modell mit 126kbp-Schleifen. Kirin, ciu., 

Subkompartimente. “ cnleifen mischen sich nicht, bilden aber auch keine 

Fehler der Observablen: 

Fur die Analyse warden die Observablen hk a- 

Simulation gemittelt. Die Fehler der ■ , er ^ versc ^ lc denen Konfigurationen einer 

Standardabweichung dieser Mitteluno S en Observablen ergeben sich nun aus der 

Mittelung, dividiert durch die Wurzel der Anzahl der unabhangi- 
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FromBegin in Abhangigkeit von der Linkerlange und EV-WW links fur die Cis-Anordung und rechts fur 
die Trans- Anordnung (die vier unteren, wellenformigen Kurven . 

gen Konfigurationen, iiber die der Mittelwert gebildet wurde. Untersuchungen ergaben, daB 
der Fehler in der GroBenordnung von weniger als 3% des Wertes der Observablen liegt. Da 
dies sowohl durch die Testsimulationen, als auch durch Stichproben in diesen und 
Simulationen von Miinkel bestatigt wurde, wird im weiteren nur noch bei besonderem AnlaB 
auf Observablenfehler eingegangen. 

Allgemeine Aussage iiber die Monte-Carlo -„Dynamik“: 

Schon wahrend der Durchfiihrung unterschieden sich Simulationen mit hoher und niedriger 
EV-WW stark in ihrer Monte-Carlo-„Dynamik“ Bei geringer EV-WW (Abb.III.10, oben) 
konnte die Chromatinfaser relativ leicht andere Chromatinfasem durchdringen oder in andere 
Chromosomen-Territorien eindringen. Weiterhin zeigte sie hohe raumliche Flexibility und 
eine relativ groBe Standardabweichung der Verteilung des Gyrationsradius. Hohe EV-WW 
fiihrte dagegen zu wesentlich kompakteren Strukturen mit statischerer Dynamik (geringer 
Standardabweichung, Abb.in.9). Beim MLS-Modell bildeten sich kugelformige Strukturen mit 
wenig stmkturierter bzw. zuganglicher Oberflache (Abb. III. 10, o.r.). Falls es Teilen des 
Subkompartimentes gelungen war, sich weit genug aus dem Potentialtopf zu entfemen, so 
blieben sie dort quasi gefangen. 

Die unterschiedliche Form und GroBe der Anfangsstruktur fiihrte hierbei zu unterschiedlicher 
Enwicklung des mittleren Gyrationsradius: groBe Anfangskonfigurationen zeigten abneh- 
mende Gyrationsradien und umgekehrt. Im Mittel wurden ca. 50.000 Schritte gebraucht, um 
den Mittelwert zu erreichen. Fiir die mitunter groBen Schwankungen des Gyrationsradius 
(Abb.in.9, o. r., und zum Vergleich idem, o.l., u. 1., u. r.) konnte bisher noch keine Erklarung 
gefunden werden. 

IIL5.1 Analyse des MLS-Modells 

Da die Simulationen in zwei Felder von Parametem aufgespalten wurde, soli diese Trennung 
auch hier durchgefiihrt werden: 
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Abb.m.12: 


RandomPairs von oben nach unten mit verschiedener Auflosung, links fiir niedrige und rechts fiir hohe 
EV-WW, bzw. die Cis-Anordnung (sie ist jeweils in dem helieren Farbton dargestellt). 

Feld 1 


FromBegin: 

CV «ni 7 gm 1St m Abhangigkeit vom genormschen Abstand, fiir verschiedene Linkerlangen und 

„ in HI- 10 links fiir die cis-Anordnung und rechts fiir die trans-Anordnung der 

Schleifen aufgetragen. 

EV-WW = 0,1 kT: 

dehnt TicifanT El ^ e !^ n § e nimrnt Her mittlere raumliche Abstand zu, d. h. das Chromosom 
eerineen auRerrn "7 ro ™ osom ve rhalt sich damit wie eine Polymerkette, die man bei 
Potentialtopfes gegebln " 85 C m8Ungen ver 8 r66ert - Dies ist durch die niedrige Hohe des 
EV-WW = 1,0 kT: 

£££ Abstand nich. zu, da die Hbbe des 

S“*ES - ^ Pedcde. Das Cheotnosot. 

Nachbarschaft zwischen 7 ™,-; ru 0rus zusammen > wobei das Zustandekommen einer 

das Chromosom sich von seinem °^ SOmenteilen statistisch ist. Visualisierungen zeigen, da 6 
dann in die Torusform faltet Anfor ^ n ^ sz ^ stan< ^ ausdehnt, an den Topfrand stoBt und sich 
Teile des Chromosoms kaum mehr start c«h r* 16 " EV ' WW fl nden nun Umlagerungen groBer 
Generell liegen die Abstande niedno^r 1 u a _ eme Relaxation nicht mehr erfolgen kann. 

Die Ausschlage zu hoh^We^nT EV 'WW =0,1 kT. 

Chromosoms sich iiber den Rand . “Hterlange v on 2400 nm zeigen, da B ein Teil des 

dCS Potentia ltopfes hinausbegeben hat und dort verblieben 
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Radial er Abstand zum Massenschwerpunkt [p.m] 

Abb.m.13: 



Radiale Aufenthaltswahrscheinlichkeit links und Dichte rechts, oben fiir niedrige, unten fiir hohe EV-WW, 
fiir verschiedene Linkerlangen, wobei die helleren Fartone die Cis-Andordnung darstellen. 


ist Dies ist dadurch erklarbar, daB die benotigte Energieanderung viel zu hoch ist, so daB die 
Akzeplanzrate des Monte-Carlo-Verfahrens zu stark absinkt. 

Ob dieses Verhalten und die zufallige periodische Schwankung ein intrinsischer Effekt der 
hohen EV-WW oder zu geringer Simulationsdauer ist, ist bisher unentscheidbar. Der Effekt ist 
jedoch noch bei 750.000 Monte-Carlo-Schritten festzustellen, wogegen bei niedriger EV-WW 
400.000 Schritte ausreichen, um das Chromosom zu relaxieren. 

RandomPairs: 

RandomPairs ist fiir verschiedene Linkerlangen und Schleifenanordnungen fiir niedrige EV- 
WW (0,1 kT) in Abb.ni.12, linke Seite, und fiir hohe EV-WW in Abb.m.12, rechte Seite, 
darlegt. Durch die enge Verwandschaft von FromBegin und RandomPairs lassen sich die 
Begriindungen und Folgerungen von FromBegin iibemehmen: 

EV-WW =0,lkT: 

Die raumlichen Abstande nehmen mit wachsender Linkerlange zu . 

Es zeigt sich jedoch eine Abhangigkeit von der Schleifenanordnung: Die cis-Anordnung hefert 
groBere Abstande fiir Linkerlangen ab 1200nm (Abb. HI. 12, links unten). Eine Erklarung gibt 
es hierfiir bisher nicht. 

EV-WW = 1,0 kT: 

Die raumlichen Abstande nehmen mit wachsender Linkerlange nicht zu. 

Aufgrund der Mittelwertbildung ist das Uberschreiten des Potentialtopfes nicht mehr 
lesizustellen. Eine Abhangikeit von der Schleifenanordnung ist nicht zu erkennen (Abb. HI. 12, 
unten). 
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Abb.m.14: 

Segmentlangenverteilung fur links niedrige und rechts hohes EV-WW, fiir verschiedene Linkerlangen. 


Generell sieht man unterhalb von IMbp bereits die Schleifenstruktur des Modells. Durch die 
Erhohung der Auflosung in den Nachrelaxationen wird dieser Effekt dann starker ausgepragt 
sein, wie Testsimulationen zeigten. Fiir kleine genomische Abstande und hohe EV-WW 
resultieren wiederum niedrigere Abstande. Unterhalb von IMbp liegen die Abstande von 
Konfigurationen mit Linkerlangen oberhalb von 2400 nm iiber denjenigen fiir niedriges EV- 
WW und einer Linkerlange von 1200 nm. 

Radiale Aufenthaltswahrscheinlichkeit und Dichte: 

Aus unterschiedlichen raumliche Abstanden wie FromBegin und RandomPairs resultiert auch 
eine Veranderung der radialen Aufenthaltswahrscheinlichkeit und Dichte. 


EV-WW = 0,1 kT: 

Das Maximum der Aufenthaltswahrscheinlichkeit nimmt mit zunehmender Linkerlange ab 
(Abb.m.13, oben links). Es verschiebt sich dabei zu groBeren radialen Werten, und die Breite 
der Verteilung nimmt zu, da das Chromosom ausgedehnter ist (FromBegin, RandomPairs) 

Das Plateau konstanter Dichte (Abb.III.13, rechts oben) nimmt bei einer Erhohung der 

Linkerlange ebenfalls ab. Kompensiert wird dies durch einen flacheren und langeren Abfall 
der Dichte zu groBeren Radien hin. 

Der Durchmesser (soweit man bei einem so dynamischen Gebilde hiervon iiberhaupt sprechen 
0,::Z^°rT..l 5b ^ somit ^so zwischen 4 nm und 5 (am. Ein Vergleich beider 
des Chromosf ZU1 ^ h" U d* 1 " aus ’ vor a ^ em bei groBer Linkerlange der groBte Teil 

MhMe ten von H r fufT befindet und frei zu ^g |ich ist- Die bereits beschriebenen 
V ° n dCr ^etfenanordnung lassen sich in unausgepragter Form ebenfalls 

EV-WW = 1,0 kT: 

“ isen (Abb.in.i3. 

D.ch B is, hie, „o«o hcmh. befu > den ** W % des Chromosoms, die 

etwas mehr als die mittlere Ausdehnnn h d " 8 1St untcr 1,5 ^ beschrankt. Diese betragt 

des Chromosom genau auf sem Volimi SehSIlS z“ be ' ho1 '" EV ' WW M S lich “! 
RandomPairs wird nun klar, daB das f'h n ' C m Zusammenhang mit FromBegin und 
sondem auch gegen die Wand es PotenriuhT'f ° S ° m SlCh mcht nur zu einem Torus einfaltet, 
im Randbereich ist im Innem die Dichte oe ^ CS ^ e P resst w * rc * Durch die starke Konzentration 
Durch den scharfen Peak ist die RanH"^ 

,43berflachenschicht“ entsteht. ZOne sehr dilnn, so daB eine unzugangliche 
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MLS-Feld 2: FromBegin oben links, RandomPairs oben rechts, radiale Aufenthaltswahrscheinbchkeit 
Mitte links, radiale Dichte Mitte rechts, Segmentlangenverteilungen fur 200nm-Segmente unten links und 
300nm-Segmente unten rechts. 


Segmentlange: 

EV-WW = 0,1 kT: 

Die Segmente sind fur niedriges EV-WW (Abb. III. 14, links) sehr gut relaxiert, da die 
Verteilung scharf und symmetrisch ist. 

EV-WW = 1,0 kT: 

Bei hoher EV-WW (Abb.m.14, rechts) zeigt sich jedoch eine schlechte Relaxierung, da die 
Verteilung zu kleinen Segmentlangen unsymmetrisch verschoben sind. Der Chromatinfaden 
steht also unter erheblichem DruckstreB, da er gegen die Wand des Potentialtopfes gedriic t 
wird, wie die Dichte anschaulich zeigte. Die schlechte Relaxierung deutet dabei auf erne 
Simulierung mit zu wenigen Monte-Carlo-Schritten hin. 
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Feld 2: 

Abhangiekeit von der SchleifengroBe: 

Aus FromBegin (Abb.ID15, oben links) folgt, daB mit wachsender SchleifengroBe und damit 
abnehmender Schleifenzahl die Abstande zunehraen. RandomPairs zeigt diesen Effekt noch 
deutlicher (Abb.III.15, oben rechts). Da jedoch in Feld 2 die Linkerlange konstant gehalten 
wurde, laBt sich dies nur auf eine Wechselwirkung zwischen den Rosetten zuriickfiihren 
Zunachst sollte rnan jedoch annehmen, daB Rosetten mit vielen kleinen Schleifen 
abgegrenztere Einheiten darstellen und wenige groBe Schleifen sich mehr inneinander 
verhaken konnen. Das Umgekehrte scheint aber der Fall zu sein. Eine Erklarung gibt es hierfiir 
bisher nicht. 

Die radiale Aufenthaltswahrscheinlichkeit (Abb.III.15, Mitte links) und die Dichte (Abb.ni 15 
Mitte rechts) nehmen dementsprechend ab. Dieser Effekt wird teilweise durch die Erhohuns 
der ChromosomengroBe kompensiert. 

Abhangigkeit von der ChromosomengroBe: 

Mit wachsender ChromosomengroBe nehmen die Abstande und die Dichte zu (Abb III 15 
oben und mitte rechts), da der Potentialtopf sich nicht verandert hatte. Die Wechselwirkung 

der Rosetten, laBt sich also weder fur lokales, noch fiir globales Oder halbelobales Verhaken 
vemachlassigen. 

Abhangigkeit von der Segmentlange: 

Die Segmentlangenverteilung zeigt fur 200nm und 300nm gleich gute Relaxation. 


m.5.2 Analyse des RWGL-Modells 
FromBegin: 

EV-WW =0,1 kT: 

Fur verschiedene SchleifengroBen in tAhK it i* l 

ahnhches globales Verhalten wie h,p c- , '. 6 ’ oben und links unten) zeigt FromBegin 

Linker fur RGWL- und MLS- Model l des MLS-Modells, da der DNA-Gehalt im 

Linker, aber nicht die Schleifen ProRe • ~ Cb warden. Dies bedeutet, daB nur der 

Fiir kleine genomische Abstande erlST* EmfluB a “ f globale s Verhalten haben. 
chleifen handelt. Die scheinbare Ve ui 30 6 06 starIce Modulation, da es sich um groBe 
unscharfere Mittelung fiir groBe gennmi Jn Tl. 6 ™ 8 der Schleifen laBt sich durch immer 
VergrbBert man die Linker W e " f AbstSnde ^en. 

Verh A |! SUndeS festzuste Ben. Fiir hohe^y 1 ww 56 * 01 MLS 'M°dell ein generelles Anwachsen 
Verhalten wie in Feld 1 der Simulationen f lgt Abbmi 6, rechts unten, ein ahnliches 

-Modells, allerdings ohne ausgepragte 
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Abb.ni.16: 

FromBegin fur verschiedene SchleifengroBen, oben und links unten mit niedrigem und rechts unten mit 
hohem EV-WW. Die Farbskala rechts geht von gelb iiber blau nach grtin und folgt so Tab.A.3 von oben 
nach unten in ihrer Bezeichnungsweise. 

Periodizitat.. Wieviele Monte-Carlo-Schritte fur eine Relaxation benotigt werden, bleibt 
wiederum unklar. 

RandomPairs: 

EV-WW = 0,1 kT: 

Fiir niedrige EV-WW (Abb.IH.17, linke Spalte) und verschiedene SchleifengroBen verhalt sich 
RandomPairs fiir groBe genomische Abstande ebenfalls ahnlich zu Feld 1 aus den 
Simulationen des MLS-Modells. Auf kleineren Skalen ist jedoch klar die starke Abhangigkeit 
von der SchleifengroBe und die Schleifenstruktur zu erkennen. Ein deutlicher Effekt der Form 
der Schleife und des Linkers und somit der Ausgangskonfiguration ist zu erkennen. Die 
Steilheit des Anstieges fiir kleine genomische Abstande ist mit dem Random- Walk einer freien 
Kette erklarbar. Die Zwangsbedingungen fiir die Segmente in der Schleife nimmt dabei mit 
steigender SchleifengroBe ab. Die VergroBerung des Abstandes mit steigender Linkerlange 
verhalt sich analog zu Feld 1 der MLS-Simulationen. 

EV-WW = 1,0 kT: 

Fur hohe EV-WW ergeben sich grundsatzlich groBere Abstande als bei niedriger EV-WW. 
Wie in FromBegin schon erlautert, falten sich die Anfangskonfigurationen von auBerhalb in 


HI Simulationen 


50 



ndomPairs mit verschiedenen SchleifengroBen, links fiir niedriges, rechts fur hohe EV-WW. 


VduroeTS P? ^ SiCh “ ^ 

verstarkt. g * ur S ro ^ ere Linkerlangen, was diesen Effekt noch 


Rj-M. **---*«***». „„„ Djditg. 

niedrige EV-WW fUr verschiedcne SchleifengroBen und 

SchleifengrSBe wachst der Maximumw!^ vAbb.III 18 .oben links). Mit abnehmender 

Linker gegenuber dem Anfaneswert ok, 4 "’ 'I™ fUr ^ eine Schleifen wieder mit groBerero 
(Abb.ni.18, oben rechts). Der Effekt ist dar^f mCn ' Fiir d ' e Fichte gilt entsprechendes 
ichter in den Zwischenraumen anderer Srhi . zurticlczu fuhren, daB kleinere Schleifen sich 
^Weifen auch ein ausgepragteres PlatM e ” au ^ ia ^ ten konnen. Hierdurch entsteht bei 
“ ***** »"Sdehnen die sich such einf.cher 

■ deshalb da, Plareau konslan, ab. FUr den 
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Abb.ffl.18: 

Radiate Aufenthaltswahrscheinlichkeit, links, und Dichte, rechts, flir oben niedrige und unten hohe EV- 
WW fiir verschiedene SchleifengroBen. 


Durchmesser eines Chroraosoraenterritoriums ergeben sich somit auch etwas groBere Werte 
als in Feld 1 der MLS-Simulationen. 

EV-WW = 10 kT: 

Fiir hohe EV-WW tritt eine Starke Konzentrierung des Chromosoms auf (Abb.in.18, links 
unten). Die Aufenthaltswahrscheinlichkeit besitzt ein Maximum, das, wie bei den 
Simulationen des MLS-Modells, bei 1,24 |xm liegt. Daran schlieBt ein zweites Maximum oder 
Plateau an. 

Die Dichte (Abb.in.18, rechts unten) zeigt eine starkere Plateaubildung mit einer Kante, bzw. 
bei scharfem und hohem ersten Maximum einen leichten Peak bei 1,24 pm. Innerhalb der 1,24 
|im-Zone tritt dabei eine gleichmaBige Dichte auf, die etwas starker als bei niedriger EV-WW 
abfallt. Die hohe Aufenthaltswahrscheinlichkeit auBerhalb dieser Zone ist auf die lange Zeit 
zuruckzufuhren, die benotigt wird, um sich aus der Anfangskonfiguration in durch andere 
Schleifen freigegebene Bereiche hineinzufalten. Anhand von Visualisierungen kann gezeigt 
werden, daB die Schleifen dabei sehr dynamisches Verhalten zeigen, d. h. sie bewgen sich 
stark. Im Vergleich zu niedrigem EV-WW gibt es keine systematischen Abhangigkeiten. 


Segmendange: 

EV-WW = 0,1 kT: 

Fur hohes und niedriges EV-WW und verschiedene Schleifenlange ergeben sich gegenUber 
dem MLS-Modell scharfere Segmentlangenverteilungen (Abb.D3.19). Die Segmente sind also 
sehr gut relaxiert und stehen unter wesentlich geringerem StreB. Dieser nimmt mit sinkender 
SchleifengroBe und zunehmenden topologischen Randbedingungen jedoch leicht zu. Eine klare 
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Abb.m.19: 

Segmentlangenverteiliuigen fur links niedrige und rechts hohe EV-WW fiir verschiedene SchleifengroBen. 

Systematik ist aber nicht zu erkennen. Ein aussagekraftiger Unterschied zwischen hoher und 
niedriger EV-WW ist ebenfaUs nicht festzustellen, da groBe Schleifen grundsatzlich besser 
SegmentstreB abbauen konnen. 


HI.5J Multifraktaiitat von Chromosomen 

Aus Abb.m.9 ergibt sich, daB die geometrische Anordnung der Chromatinfaser ein auBerst 
komplexes Verhalten aufweist, wie es typisch fur die in Abschnitt n.2 beschriebenen fraktalen 
Systeme ist. Hierzu wurde die fraktale Dimension aus den Log-Log Kurven bestimmt. Der 
Werzu benotigte Fit (KaleidaGraph, Version 3.08 fiir Macintosh) wurde so optimiert, daB der 

faktote SeSif ent m5gliChSt gr ° B Wurde ' Die ne 8 ative Steigung des linearen Fits ist die 
Allgemein wurden folgende Phanomene fiir alle Simulationszyklen festgestellt: 

Euklidischer Cut-Off: 

Dimei^L ^ Se S mentla nge, dannwird die fraktale 

Da die Segmentlanee 300 nm e “ er Luiie mit fraktaler Dimension 1) bestimmt. 

fraktale Dimension von 1 05 bis 1 08 ^ T" 60 Ma6bereich von 5 nm bis 300 nm eine 

ein kleiner Anstieg in der n f fu " den - Kurz vor dem MaBwert 300 nm befindet sich 

messen sie dort die Segmente besonH 16 MaBstaB e selbst ganze euklidische MaBe darstellen, 

der fraktalen Dimension von 1 ist auf d™ v** DlC Breite des Anstiegs und die Abweichung 

auf die Verteilung der Segmentlangen zuriickzufuhren. 

Korrelationskoeffizient war nicht zu^rr» meiSten FSllen 0,999±0,001. Ein solch hoher 

mogliche Anzahl von Punkten verw^nH t n ^ ewesen * FUr die Anpassung wurde die maximal 
^*nbreite einer Kurvendarstellune M^ot \ * der Fetl * er eines MeBpunktes unterhalb der 
Besummung der fraktalen Dimension nicht’ ^ f schlossen werden, daB der Fehler bei der 
w ° n mcht mehr als ± 0,01 betragt 

Multifraktaiitat: ’ 

Objetami^e" 6 " Uberdecken Dunension zwischen 300 nm und 5000 nm 

remen mathematischen Sinne uLdhctf 1 " 6 , Gr5Ben o^nng. Obwohl ein fraktales 

viele GroBenordnungen iiberdeckt, sollte 
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Bestimmung der fraktalen Dimension fiir links niedrige und rechts hohe EV- WW, fur das MLS-Feld 1. 


dies fiir reale Objekte zumindest flir mehrere GroGenordnungen gelten. Avenir konnte zeigen, 
daB im aligemeinen schon weniger als eine GroGenordnung ausreicht, um ein Objekt als fraktal 
zu bezeichnen und die entsprechenden Folgerungen zu ziehen [Av 89, Av 98]. 

Als Vergleichswerte seien hier schon die fraktalen Dimensionen von vemetzten Polymeren, 
Gelen und Kettenpolymeren in gutem Losungsraittel genannt. Die typische fraktale 
Dimensionen betragt bei Polymeren und Gelen 2,50 und bei Kettenpolymeren 1,67. Diese sind 
jedoch auf einem weit groGeren Skalenbereich giiltig [BrBa 86]. 

Der Nachweis des euklidischen Cut-Offs und die Moglichkeit, die Fits mit solch hohen 
Korrelationskoeffizienten zu bestimmen, stellt die Leistungsfahigkeit einer fraktalen Analyse 
dar und zeigt, daG das Verfahren korrekt implementiert wurde. 


Fraktale Analyse des MLS-Modells - Feld 2: 

Fiir niedrige EV-WW (Abb.ffl.20, links, und Tab.A.5) findet man drei Oder vier lineare Zonen 
mit verschiedener fraktaler Dimension. Die erste Zone reicht von 300 nm bis 560 nm, die 
zweite von 560 nm bis ca. 1000 nm, die dritte von ca. 1000 nm bis 2500 nm und die vierte von 
ca. 2500 bis 5000 nm. Im ersten Bereich ist das MaG dabei sensitiv auf eine Schleife, im 
zweiten auf eine Rosette, im dritten mehrere Rosetten und im vierten schlieGlich fast ein 
ganzes Subcompartment. Die unterschiedlichen fraktalen Dimensionen entsprechen also der 
fraktalen Dimension der jeweiligen strukturellen Einheit. Mit wachsender Linkerlange nimmt 
die fraktale Dimension ab, da sie ein MaG fiir die Raumerfiillung darstellt. Die oben 
beschriebenen Observablenabhangigkeiten sind also mit der fraktalen Analyse koharent. 
Allerdings ist keine Abhangigkeit von der Schleifenanordnung festzustellen. Betrachtet man 
die Werte genauer, so ergeben sich fiir Zone I Dimensionen von 2,67 bis 2,20, fur Zone II 
Dimensionen von 3,13 bis 2,42, fiir Zone HI Dimensionen von 2.07 bis 1,98 und fiir Zone IV, 
sofem sie auftritt, Dimensionen von 2,43 bis 2,15. Rosetten stellen im aligemeinen also au erst 
dichte Gebilde dar, wahrend die Rosetten im Subkompartment recht locker verteilt sind. 

Fiir hohe EV-WW (Abb.m.20, rechts) ergeben sich jedoch neben den zwei GTSt ™Zonen, tot 
grundsatzlich hohere Werte annehmen (Zone I von 2,85 bis 2,62, Zone II von , is ’ 
und der dritten Zone mit grundsatzlich niedrigeren Werten im Bereich von 1,70 bi&W, ,emt 
vierte Zone mit sehr hohen Werten fiir die fraktale Dimension im Bereich von 4,32 1 bis WJ. 
Bei diesen hohen Werten handelt es sich um den Maximum- ut- , er * e 
Ausdehnung beschreibt. Hier verliert die fraktale Dimension ihre ii tig ™ ’ 

Umstanden konnen hier jedoch Aussagen iiber den Rand ernes Objektes getro 
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Abb.ni.21: 

Bestimmung der fraktalen Dimension 


fur das MLS-Feld 


2 . 




2.3 2.4 2 J5 2.6 2.7 2.8 2.9 3 3.1 3 2 3.3 3.4 35 3.6 3.7 

Logarithm us des MaBcs [am] 


Abb.EU.22: 

Bestimmung der fraktalen Dimension fiir links niedrige und rechts bohe EV- WW, fiir das RWGL-Modell. 


Die dritte Zone mit ihren niedrigen Werten bestatigt die oben gemachten Aussagen iiber die 
Aufenthaltswarscheinlichkeit und die Dichte (Abb.IH.16 und S.48). 

Die hohen fraktalen Dimensionen in Zone H ermangeln jedoch bisher einer Erklarung, es 
konnte sein, daB hier ein Maximum-Cut-Off der Rosetten vorliegt. 


Wie nach Feld 1 nicht anders zu erwarten, zeigt sich bei Feld H ebenfalls eine Zonierung der 
fraktalen Dimension in drei Bereiche (Abb. III. 21 und Tab. A. 6). In Zone II ist das 
Zusammenpressen der Schleifen aufgrund der groBen Schleifenanzahl in einer Rosette mit 
einer SchleifengroBe besonders gut zu sehen. Weiterhin ist aufgrund der unterschiedlichen 
Se^entlangen yon 200nm und 300nm der Anstieg auch bei 200nm zu finden. Durch die 
nnH e £ men Jdigenverteilung ergeben sich nun noch zwei kleinere Anstiege bei 100 nm 

der ^ ^ es P ea ks in den Kurven bestatigt. Die Abhangigkeiten von 

to SchleifengroBe und dem DNA-Gehal, sind wiederum eindeutig nachzuweisen 

Irak tale Analyse des RWGL-Modells: 

Zonen, die sich^^ef kleinYr Schl ^ 22 T* Tab A 7) er S eben sich genereU zwei lineare 
Dimension reichen in Zone I von 1 7/^^ “ drd aufs P dten - Die Werte der fraktalen 
SchleifengroBe von 126kb D ereeh™’ 76 u u 2 ’ 23 ’ U " d m Zone n von 3 - 05 bis 2,43. FUr eine 
Linkerlange bedingt sind. SlC ierbe * jedoch andere Werte, die durch die groBere 
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di“e“™o^2 r 0?‘ n ST ! d ' r | J°“" 1 u " d 11 « ""=**«' EV-WW. «, sLellt ™. rest, 

, • k t \ 4 1-1 ^ ei ^ ro ^ en Schleifen auf ca. 1000 nm bei kleinen Schleifen 

verschiebt Man lost also w.ederum verschieden groBe Strukturen auf, wobei Xrdings ^ 
L inkers truktur und die Schleifenstruktur nicht mehr so eindeutig trennbar sind 8 

Bei hober EV-WW (Abb.m.22 reehts) bleibt die Zonengren^doch s"^ erhalten (unterer 
ert nur bei LOOnm) obwohl die fraktalen Dimensionen sich kaura andem. Bei kleinen 
Sehleifen sind jedoch wie in Feld 1 die drei Zonen gemaB den Schleifen, mehreren Sch S 
und dem ganzen Subkompartment hin organisiert. 


111^5.4 Vergleich der Modelle 

AbschlieBend soil nun noch ein Vergleich zwischen den Simulationen im Hinblick auf ihre 
Verwertbarkeit beziiglich der experimentellen Daten erfolgen. Es steilt sich die Frage, ob die 
Modelle mit den vorhandenen Daten gegeneinander abgrenzbar sind. Da ihre statistische 
Relevanz eindeutiger ist, wurden hierfiir nur Simulationen mit niedriger EV-WW verwendet. 
Die in den Abbildungen gezeigten Kurven stellen einen reprasentativen Querschnitt dar. 
Sowohl FromBegin-Abstande als auch die Segmentlangenverteilung werden im folgenden 
nicht betrachtet, weil sie keine neuen Erkenntnisse iiber die aus obiger Besprechung 
resultierenden Ergebnisse erbringen. 

Die Analyse der TAK-Modelle ergab, daB die Variation des DNA-Gehalts der einzelnen 
Banden keine anderen Ergebnisse in den Observablen zeigte, wie beim RWGL-Modell 
unterscheidet. 

Die Simulation der linearen Kette entsprach den theoretischen Vorhersagen. Auf kleinen 
Skalen folgte sie einem exakten, auf groBen Skalen einem, durch den Potentialtopf 
eingeschrankten, Random- Walk, so daB Simulationsartefakte ausgeschlossen werden konnen. 
Die fraktale Dimension der linearen Kette besaB dabei, wie das RWGL-Modell, zwei Zonen. 
In Zone II verhalt sie sich exakt nach der theoretischen Erwartung [Man 83]. Die Zunahme der 
Dichte auf groBeren Skalen spiegelt sich in der groBeren fraktalen Dimension wider. 


RandomPairs: 

Die globale raumliche Abstandsabhangigkeit vom genomischen Abstand laBt keine 
Unterscheidung zwischen den Modellen zu (Abb.II.23, links oben). Fur kleine genomische 
Abstande konnen die Modelle jedoch separiert werden (Abb. n. 23, reehts oben, und Abb.II.23, 
links Mitte). Dieser Unterschied ist allerdings von der SchleifengroBe und der Linkerlange 
abhangig. Fur groBe Schleifen im RWGL-Modell ergibt sich ein klarer Unterschied zum MLS- 
Modell mit alien verschiedenen Linkerlangen fiir genomische Abstande unter 3Mbp. 
Unterschreiten die Schleifen im RWGL-Modell 500kbp, so ist dies jedoch nur noch 
eingeschrankt moglich. Benicksichtigt man insbesondere die zu groBe Linkerlange in der 
RWGL-Modell-Simulation mit Schleifen von 126kbp, so ist dies nur noch mit lokalen 
Observablen moglich, nicht mehr jedoch mit der Observablen RandomPairs. 

Dies gilt auch fiir raumliche Abstande der linearen Kette, die fiir kleine genomische Abstande 
wie bei groBen Schleifen im RWGL-Modell anwachsen. 

Radiale Aufenthaltswahrscheinlichkeit und Dichte: 

Weder die radiale Aufenthaltwahrscheinlichkeit, noch die radiale Dichte, wie sie in AbbJI.-3 
links unten und Abb. M. 23 reehts unten gezeigt sind, erlauben es grundsatzlich, zwischen den 
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Abb.III.23: 

Vergleich der verschiedenen Modelle. RandomPairs, oben rechts und links, sowie mitte links; Fraktale 
Dimension, mitte rechts; Aufenthaltswahrscheinlichkeit, unten links; Dichte, unten rechts. 

Fraktale Dimension: 

^ ra ^ ta * en D ime nsion ist eine Unterscheidung zwischen den Modellen iiber den 
^ ka ^ en ^ ereic ^ m ogli c h. Dabei ist zu beriicksichtigen, dafi nicht nurder 
Grundsatzlich sttm lmens * on von Bedeutung ist, sondem auch ihr GiiltigkeitsintervaU. 
dem RWGL Model 16 fl ? ktale die beste Differenzierungsmoglichkeit zwischen 

Au” al uber f n MLS - Moddl d ar. Sie bietet weiterhin die Moglichkeit, direkt 
Territorium zu machen^A 11 ^ 110 ^ Und die ^echselwirkungen in einem Chromosomen- 
mehrerer Observablen zu 616 ^ serv able lassen dies nicht oder nur im Zusammenspie 
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IV Experimente 


Z.el des expenmentellen Teils war es, dreidimensionale raumliche Abstande zwischen 
fluoreszent markierten DNA-Sequenzen zu bestimmen. Mit diesen Werten wurde im VeSch 
imt den Simulationen versucht, zwischen dem RWGL- und dem MLS-Modell zu 
u^rscheiden. Die genomischen Abstande dieser Sequenzen sollten dabei weit unteriuSb von 

euiem neuen konfokalen Laser-Scanning Mikroskop 
(FISHWe^ndef h ^ ^ Markierung wurde die Fluoreszenz-m-Situ-Hybndisierung 

Mehrere hundert Zellkeme wurden aufgenommen, da bisher existierende Experimente nur 
unzurmchendestatistische Aussagekraft besitzen. Als genomische Region wurde die Prader- 
Labhart-WiUi/Angelmann-Region (PLWS/AS) auf Chromosom 15 gewahlt. Im Vergleich mit 
Simulationen sollte daraus die Faltungsstruktur der PLWS/AS-Region bestimmt werden 


IV. 1 Experimentelle Methoden 

Wie schon in Kapitel II erwahnt wurde, scheiterte die Aufklarung der Organisation von 
Chromosomenterritorien bisher hauptsachlich daran, da6 bildgebende Verfahren mit hoher 
Auflosung strukturzerstorende Praparationsschritte benotigen. Die groBten Fortschritte der 
letzten Jahre auf dem Gebiet der Chromosomenorganisation wurden durch die konfokale 
Laser-Scanning-Mikroskopie in Zusammenspiel mit der Fluoreszenz-in-situ-Hybridisierung 
gewonnen. 


IV.1.1 FIuoreszenz-in-situ-Hybridisierung 

Die in-situ-Hybridisierung ist eine Methode, mit der man spezifisch DNA- 
Einzelstrangsequenzen markieren kann. Hierzu werden kurze Einzelstrangsequenzen beniitzt, 
die ihr komplementares Gegenstiick finden und sich daran anlagem. Die in-situ-Hybridisierung 
eroffnet dadurch die Moglichkeit, mit Nukleinsauresonden spezifische Sequenzen in Geweben, 
Zellen, Zellkernen oder Chromosomen sichtbar zu machen. Die ersten in-situ- 
Hybridisierungen wurden von J. G. Gall und M. L. Pardue [GaPa 69] sowie H. A. John [Jo 69] 
zum Nachweis von rDNA-Sequenzen in cytologischen Praparaten mit radioaktiv markierten 
RNA-Sonden durchgefiihrt. Die erste erfolgreiche in situ-Hybridisierung einer Einzel- 
kopiesequenz gelang Harper und Saunders [HaSa 81], die mittels radioaktiv markierter Sonden 
das menschliche Insulingen auf Chromosom 1 nachwiesen. 

Erne Hybridisierung umfaBt fiinf verschiedene Arbeitskomplexe: 

Herstellung des biologischen Praparats, Gewinnung sowie Markierung der Sonde, 
Hybridisierung, sowie Detektion und Darstellung durch mikroskopische Verfahren. Wahrend 
in diesem Abschnitt auf die ersten drei Punkte naher eingegangen wird, werden die letzten 
beiden in Abs. IV.2 und IV.4 beschrieben. 
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DNA-Doppelstrang 


Markierte DNA-Probe 


ACG AACGG AGCA 
A gcaa_gcca_cg_aa 


AACGG 

GCC 


I 


Hybridisierung und Detektion 
der Probe mit fluoreszenz- 
markierten Antikorpern 


ACG AACGG AGCA 
T1GCC 



Abb.IV.l: 


Schema der Hybridisierung: Ein DNA-Doppelstrang besitzt eine bestimmte Nukleotid-Sequenz. Eine 
markierte DNA-Sonde besitzt eine komplementare Sequenz. Nach Denaturieren des Doppelstranges und 
der Sonde kann diese ihr komplementares Gegenstiick finden. Beide bilden dann ein Hybrid. 


Der Erfolg einer in-situ-Hybridisierung hangt von einer moglichst guten Konservierung der 
raumlichen Lagebeziehung der Zellkomponenten ab. Die Fixierung begrenzt dabei nicht nur 
den Verlust von Struktur, sondem muB das konservierte Gewebe auch fur die Sonde und die 
Nachweisreagentien genugend permeabel lassen. Je hoher der Grad der Fixierung, desto 
undurchlassiger ist das Gewebe. Die Wahl des Fixiermittels hangt von den Zellinien, von der 
Sonde, von dem fur den Nachweis der Hybridisierung bildgebenden Verfahren und der 
gewiinschten Sensitivitat ab. Im allgemeinen unterscheidet man zwei Systeme von 
Fixiermitteln, wobei die erstere die Praparate quervemetzt (Glutaralaldehyd, Formaldehyd) 
und die zweite Proteine ausfallt (Ethanol oder Methanol). Die erste Fixierungsmethode fiihrt 
o t zu genngerer Permeabilitat, die dann durch zusatzhche Schritte erhoht werden muB.Die in 
leser Arbeit verwendeten FixierungsprotokoUe sind in Abb.IV.2 niedergelegt. 

Gewinnung und Herstellung der Sonden: 


Fixierung: 
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Fixierung der Zellen - 2D Praparation: 

1) Losen der Zellen in der Kulturflasche mit 0,05% Tripsin in Phosphat Buffered Saline (PBS) 
(5min). Zellen in 50ml Rohrchen sammeln. 

2) Zellen 10 min abzentrifugieren bei 1 100 bis 1300 U/min. 

3) Hypertoniebehandlung mit 10ml 0.55% KC1 und 1% NaCitrat fur 10 min bei 37°C. 

4) Zellen 10 min abzentrifugieren bei 1100 bis 1300 U/min. 

5) 3*Fixierung mit Methanoleisessig (3:1 in PBS), waschen und zentrifugieren. 

6) Pepsinverdau: 50ftg/ml in 0.01 M HC1 (pH 2,3) (99ml H 2 0+lml 1M HC1 + 50 pi Pepsin) 
vorgewarmt auf 37°C flir lOmin bei 37°C. 

7) Auftropfen der Zellen auf die Objekttrager, wobei das Zytoplasma verlorengeht. 

Fixierung der Zellen - 3D Praparation: 

1) 3* Waschen in PBS. 

2) Fixierung in 4% Paraformaldehyd (frisch, 4g in 100ml PBS,pH 7,0) fur lOmin. 

3) Waschen in PBS 3*3 min. 

4) Membranpermeabilitat: 0.5% Saponin + 0.5% TritonX-100 fur 20min. 

5) Einlegen in 20% Glycerin (20ml auf 80ml PBS) fur 30min. 

6) Langsam in flussigen Stickstoff tauchen, auftauen lassen und zuriick zum Schritt 5). 

5) Wiederholung 5 mal. 

7) Waschen in PBS 2 * 5min. 

Abb.IV.2: 

Fixierungsprotokolle fur die 2D Praparation und die 3D Preparation 

Die Sonden in dieser Arbeit wurden mit Nick-Translation markiert. Die Enzyme DNasel und 
die DNA-Polymerase aus Escherichia coli wurden zum Einbau der markierten Nucleotide in 
beide Strange der DNA Doppelhelix verwendet. 

Die DNasel fiihrt dabei Einzelstrangbriiche in die doppelstrangige DNA ein. Durch die 
Exonukleaseaktivitat der DNA-Polymerase werden Nukleotide in 3 Richtung vom 
Einzelstrangbruch entfemt und am anderen Ende die markierten Nukleotide eingebaut. Enthalt 
die DNA zu wenig Einzelstrangbriiche, werden die Sondenfragmente zu lang und mcht 
ausreichend markiert; fiihrt die DNasel dagegen zu viele Briiche ein, geraten die 
Sondenfragmente zu kurz. 

Denaturierung der DNA und Hybridisierung: . 

Bei der DNA:DNA-in-situ-Hybridisierung werden die Sonde und die Zielsequenz denatunert. 
Danach fmdet die Sonde ihr komplementares Gegenstiick, die Zielsequenz, bildet mit ihr einen 
Doppelstrang und hybridisiert. Die Stabilitat der Hybride hangt vom GC-Gehalt, der Lange es 
Nukleinsaurehybrids, seiner Umgebung und von Basenfehlpaarungen ab. Die Sonden e in en 
sich in einer Hybridisierungslosung, die Formamid, geloste Salze, Dextransulfat, oc a 
DNA, Natriumdodecylsufat (SDS) oder Rinderserumalbumin (BSA) enthalt. Formamid sorg 
dafiir, daB die Reaktion bei einer Temperatur stattfmden kann, die die Gewetestmktur nic 
schadigtDas in dieser Arbeit verwendete Hybridisierungsprotokoll befindet sich in Abb.I . 
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Hybridisierung: 

1) RNase(A)-Behandlung: /ccrw oaa„i , • 

-Arbeitslosung (lOOpg/ml in 2 * Salt Sodium Citrate (SSC)) 200pl pro Objektrager fur lh bei 
37°C, anschlieBend waschen: 3 * 5rnin 

2) Pepsin-Verdauung: 

— 50(ig/ml in 0.0 1M HC1 (pH 2,3) (99ml H 2 0+lml 1M HC1 + 50 pi Pepsin) vorgewamit auf 
37°C fur lOmin bei 37°C. 

— Waschen in PBS fur 2 * 5min und waschen in PBS+MgCl 2 fur 1 * 5min. 

Postflxierung in 1% Formaldehyd (100ml PBS + MgCl 2 + 2,7ml 37% Formaldehyd) fur 
10m in. 

— Waschen in PBS fur 1 * 5min und anschlieBend in 70%, 90% und 100% Ethanol einlegen. 

Chromatin-Denaturierung: 

1) 2min in 70% Formamid/SSC (auf 70°C vorgewarmt), dann jeweils 3min in eiskalte 
Ethanollosung (zunachst 70%, dann 90% und 100%) und schlieBlich trocknen lassen. 

DNA-Proben-Denaturierung: 

Zunachst bei 37°C fiir 20-30min vorwarmen, dann 5min bei 75°C denaturieren. 

Detektion: 

1) Herstellen von 11 4 * SSC mit 0.2% Tween (SSCT): 200ml 20*SSC + 800ml H,0 + 2mi 

1 ween-20. L 

2) Objektrager in SSCT + 3% BSA fur 20min bei 37° C abblocken 

3) Dererste Wechsel: 

m sscn 

4) Der zweite Wechsel: ' 

lOOpl biotinyliertes anti-Avidin D (goat 1 -900 in ccm , _ . . 

lOOpl Fluorescein- Avidin DCS (T200 in ssrr'i t- a* ■ . 

* 5min (Schuttler). ‘ SSCT) fur 45n »n bei 37°C. Waschen in SSCT fiir 3 

6) DN A-Farbu ng: 50ml of SSCT +2 s..i r>Ar»j t . lx 

all .1 I . . r 1 IJAp l-Lsp. (lmp/mh fiir Hn,;- 

rtOD.ivj. Das verwendete Hybridisieru ^ 

Beispiel fur die Detektion. ^ P r °tokoII, DNA- und Chromatin Denaturierung, sowie ein 


Detektion der Hybridisierungsstelle: 

Im AnschluB an die Hybridisie 

nachgewiesen werden. Fur de^N^hLf ^. S . 0nd f. mUssen die Hybridisierungsstellen 
rKjmalerwe |se die Immunogenitat der Markerm oft • V markierter Sonden nutzt man 

Bin "i ° r P ern ( z - B. Anti-Digoxigenin) naoh ° ° u c aus - Die Markermolektile lassen sich 
Affinitt ntlki ? er ° der mit Avidin (ein GlykT' 56 "' Bit>tin la6t sich en tweder durch Anti- 
£S£L ZU B,0Un ' K * = lO^M- 1 ) TdeTl J?*™ aus HUhnereiweiB mit sehr hoher 
Antikom f rnU . weileren AnUkorpern fuhrt 7 .. P tav,dln nachweisen. Eine Detektion der 
tSZLSL ^° ,ekU,e ^ebunden.die detekT VerStark -8 des MeOsignals. An die 

sich in Tab. IV. 1 ' ° VCWendet werden Eine Liste'gSi^^n" k6nnen ' Hierzu k6nne " 

8 g!ger Fluoreszenzfarbstoffe befmdct 
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FluoreszenzfarbstofT 

Anregungs- 
maximum Tiunl 

Emissions- 
maximum fnm] 

Far be der 
FI uoreszenz 

a) Fluorophor 




Coumarin-AMCA 

350 

450 

hi on 

Fluorescein Isothiocyanat 
(FITC) 

495 

515 

griin 

R-Ph ycoerythrin 

525 

575 

rot 

Rhodamin 

550 

575 

rot 

Tetramethy Irhodamin 
Isothiocvnat (TRIC) 

575 

600 

rot 

Texas Red 

595 

615 

rot 

Ultralights 

rot 

680 

dunkelrot 

CY3 



rot 

CY5 

648 

650 

dunkelrot 





b) DNA-Gegenfarbuns 




Chromomycin A3 

430 

570 

gelb 

DAPI 

355 

450 

blau 

Hoechst 33258 

356 

465 

blau 

Propidiumiodid 

340. 530 

615 

rot 


Zur Detektion eingesetzte Fluoreszenzfarbstoffe und ihre spektralen Eigenschaften. 


IV.1.2 Konfokale Laser-Scanning-Mikroskopie (CLSM) 

Die konventionelle Auflicht-Fluoreszenzmikroskopie wird oft durch verschwommene 
Bildbereiche beeintrachtigt, die durch unscharfe Abbildungen von Strukturen, die sich nicht in 
der Brennebene befinden, entstehen. Marvin Minsky entwickelte 1957 [Mi 57] das konfokale 
(Laser-) Scanning-Mikroskop (CLSM), um die Auflosung zu verbessem. 

Konfokal bedeutet, daB das Bild des Beleuchtungspunktes mit dera Detektionspunkt im 
Objektraum so zusammenfallt, daB sich beide in der Fokalebene iiberlagem (Abb. IV. 4). Eine 
groBere Auflosung in z-Richtung wird durch die Punktformigkeit der Beleuchtung und der 
Detektion mit einer Blende erreicht. Im Gegensatz zum Epifluoreszenzmikroskop muB der 
Fokuspunkt aber das Objekt abtasten (scannen). Die Anregungsstrahlung fallt im 
Fluoreszenzbetrieb durch die Gesamtprobe. Dadurch werden stets auch Fluoreszenzmolekiile 
angeregt, die ober- oder unterhalb des Focus liegen. Das Mikroskopobjektiv erfaBt auch diese 
Bereiche. Ihr Beitrag am detektierten Signal wird proportional zum Abstand von der 
Fokalebene durch die Blende unterdriickt. Als Detektoren werden hochempfindliche 
Photomultiplier verwendet, die noch geringe Photonenintensitaten messen konnen. Durch die 
Verringerung der Anregungsintensitat laBt sich das Ausbleichen der Fluoreszenzfarbstoffe 
stark reduzieren. 

Fiir die Intensitat einer Punktlichtquelle eines konventionellen Fluoreszenzmikroskops gilt 
[Wil89, WiCa 88]: 


si(v) = 



(IV. 1) 
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Laser 


Konfokale Blenden 


Detektor 


Objekt 
, in Fokalebene 
nicht in Fokalebene 


Abb.IV.4: 

Prinzip des konfokalen Laser-Scanning. 
Mikroskops: 

Ein Laserstrahl wird durch einen 
dichroitischen Spiegel und ein Objektiv 
auf die Probe fokussiert. Das von der 
Probe ausgehende Signal wird durch 
dasselbe Objektiv und den Spiegel auf 
den Detektor geleitet. 

Im Fokus wird dabei nur ein kleines 
Volume n beleuchtet und detektiert, das 
Anregungs- und Detektionssignal wird 
durch eine konfokale Blende gefuhrt. 


wobei J,(v) die Besselfunktion erster Ordnung der normierten optischen Einheit v ist, die sich 
wie folgt darstellt: 


v = — — rsin(a), (IV-2) 

A 

mit der Wellenlange X, dem Radius r des Unscharfekreises im Objektraum und dem halben 
Offnungswmkel a. Fur ein CLSM ergibt sich fur die Intensity 


W v)=(H(r>y, 


(IV.3) 


25% besser im Ver^ei^hz^einemk 1 P ° tenz des Abstandes ab. Die Auflosung ist damit urn ca. 
gleiCh zu euiem konvenuoneUen Mikroskop [Ste 86]. 


IV.1J Das konfokale Mikroskop Leica TCS NT 

Scanning- Mikroskop TCS* N'^ < def I p* C ° n ^°* Ca ^ en Atdnabrnen wurde das konfokale Laser- 
gesamte Mikroskopaufbau zu sehen- ™ a Leica verwendet. In Abb.IV.5 ist zunachst der 
mit dem Objekttragertisch im Zentmm n C1 ? em S^ampften Granittisch steht das Mikroskop 
seitlich angebrachten Quecksilberdam befuidel sich der Objektivrevolver. Neben den 

dem MikroskopfuB die Scannereinn^^^" fUr den D urchlichtbetrieb befindet sich auf 
Monitoren fur die Aufnahmekontr^ufe™ 8 ' dem Tisch rechts befmden sich neben den 
K steuert die Aufnahmen. Fiir die Anr^o manuelles Steuerpult fiir Mikroskopparameter. Em 

UI *g wird ein Argon-Ionenlaser benutzL 
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Abb.IV.5: 

Aufbau des Leica TCS NT (mit freundlicher Genehmigung der Firma Leica). 



Abb.IV. 6: 

Aufb.u dcs Scannerkopfes des Leica TCS NT (mit freundlicher Genehmigung der Finna Leica). 
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. , Hpr scannerkopf erlautert: Das LaserUcht wird von einem auf dem 
In AbbUdung rV.6 wud akustooptischen Filtem (A,) durch ein 

FUterrad (A) 1 ™_ Licht) in das System eingekoppelt und auf* ! 

Glasfaserop^ (Bp >■ 2 fokussiert Der von d er Blende kommende Lichtstrahl wird m 

Anregungsbler^e danach von einem dichroitischen Spiegel auf dem Filterrad ( E ) 

JSS Ein einzelner Spiegel (G), der durch zwei Galvanometer betrieben wird, lafit den 
EchSl nun in X-Y-Richtung das Praparat abtasten. Erne Lmse (H) und em Objektiv (J) 
erzeugen die punktformige Beleuchtung des Praparates_ Das von dem Praparat ausgehende 
Fluoilszenzlicht wird vom Anregungslicht durch den dichroitischen Spiegel getrennt und iiber 
dieselbe Optik wieder auf Filterrad (E) gelenkt, von wo es auf die variable konfokale 
Detektionsblende (K) fokussiert wird. Danach lenken dichroitische Spiegel (L) den Lichtstrahl 
separat in die ausgewahlten Emissionsbander auf die Photomultiplier (O). Die Filterrader (N) 
nehmen eine weitere Filterung des Signals vor. Einerseits soli laut Hersteller durch das 
kompakte Design die Lebensdauer, durch die Trennung der einzelnen Komponenten jedoch die 
Flexibility des Scannerkopfes erhoht werden. 


IV.1.4 Bildverarbeitung 

Fur die Prazisionsdistanzmessung und die volumenerhaltende Segmentation von Objekten 
nahe oder unterhalb der Auflosungsgrenze werden bei der CLSM Bildverarbeitungs- 
algorithmen benotigt.. Die in dieser Arbeit verwendeten Bildverarbeitungsmethoden wurden 
von Bomfleth [Bor 98] speziell fur diese Anwendung entwickelt. Sie sollen im folgenden kurz 
vorgestellt werden. 

Fast alien Bildverarbeitungsmethoden ist eine Rauschfilterung (hier eine Median-Filterung) 

des dreidimensionalen Bildstapels gemeinsam. Um Strukturanalysen zu betreiben, also GriiBen 

wie Objektvolumen und Distanzen zu messen, werden Segmentierungsverfahren eingesetzt. 

Fiir sehr kleine Strukturen, die in ihrem Volumen stark variieren und oft mit dem 

D. e VL t !T 8SV0 , lumen “berlappen, ist die GroBe dieser Strukturen stark von der 

. fi .. _ un ^ s . don (PSF = point spread function) des Mikroskops abhangig. Diese PSF 

durch MAM UOreSa ^ e Mikropartikel, deren Eigenschaften genugend genau bekannt sind, 

fiir das Leic/^rr c we ^ en - Oie PSF der im einzelnen verwendeten Farbstoffe wurde 

Mikrokuseln tRpaH i dlC < ^ u ^ na ^ lnie von fluoreszierenden (FITC und TRUC) 

Mikrokugeb (Beads) von Bomfleth gemessen. 

zunachst £e MaJdl^ Algonthmus von Bomfleth sucht in einem Bildstapel 

getrennt sein In einem , ^ M a mussen durch mindestestens zwei Voxel (Bildstapel) 
Objektvolumen in de^Probe^fmd^ wird T versucht, das „wahre“ 

der Intensitat (Schwellwert) zuv edeni Objekt wird dabei ein individueller Mindestwert 

Die Funktion, die beide Groften^ 16 ^.! der VOn d * esem wahren Volumen abhangt. 

Objekte ist das Verhaltnis aus Sch m » et ’ Wurde durch eine Simulation bestimmt. Fiir groBe 
Objekte ergibt sich folgende Abhangt "T* 1111(1 maximaler Intensitat konstant. Fiir kleine 

gig eit der relativen Intensitat vom Volumen 


IJV\ = ?Y 2 +pV+y 
100 


(TV.4) 


Parameter cxdc * 

Objekt ein Schwellwert vonSt^reL^K™^ bestimmt wurden. Zu Beginn wird jedem 

segenen, der iteriert wird. 
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Das Objekt wachst nach folgenden Regein: Ein Voxel gehort zura Objekt, wenn seine 

, L 6 l! ,Ch dCm SchweUwert ^ nachste Voxel auf der Verbindungslinie 
zum Objektmittelpunkt bereits zum Objekt gehort, die Intensity dieses nachsten Voxels grofer 

Oder gleichdem SchweUwert ist und schUeBUch das Voxel nicht bereits einem anderen Objekt 
im selten Farbkanal angehort Aus dem sich daraus ergebenden neuen Volumen ergibt sich ein 
neuer SchweUwert und das Verfahren wird wiederholt, bis der „wahre“ SchweUwert gefunde“ 

Daruberhinaus kann nach der Bestimmung der Schrittweite des Mikroskops (x y = 78nm und z 
= 200 nm) der Massenschwerpunkt eines Objektes bestimmt und die Distanz zwischen zwei 
Objekten berechnet werden. Fur Objekte, die in unterschiedUcher spektraler Signatur markiert 
sind, lassen sich Distanzen bestimmen, die kleiner sind als die abtesche Auflosungsgrenze 
Unterschiedliche Wellenlangen werden dabei nicht auf denselben Bildpunkt abgebUdet 
chromatischer Shift). Die Farbkanale mussen damit korrigiert werden. Zwischen FITC und 
TRITC wurde ein axialer Shift Ac z = 214nm gemessen, wahrend der Shift zwischen TRITC 
und CY5 zu vemachlassigen ist 


Der Fehler bei Distanzmessungen in einem CLSM hangt nicht nur von optischen Aberrationen 
Photonenanzahl, Detektorrauschen und ungenauer Positionierung des Scanners ab, sondem 
auch vom Ausbleichen des Fluoreszenzfarbstoffes wahrend der Messung. Die Genauigkeit der 
Objektpositionierung kann durch Messungen mit fluoreszenz-markierten Quarzpartickeln 
bestimmt werden [Bo98]. Der groBte anzunehmende Fehler betragt in axialer Richtung 85 nm, 
im besten Fall ca. 50 nm. Fur Objekte, die lateral zueinander liegen, ergeben sich im 
schlechtesten Fall 31nm und im besten Fall 15 nm. Fiir die Genauigkeit dreidimensionaler 
Distanzmessungen fiir Objekte mit gleicher spektraler Signatur erhalt man einen Wert unter 
60nm. Fiir Objekte verschiedener spektraler Signatur ergibt sich eine Genauigkeit die kleiner 
ist als 70 nm. Im ersten Fall ist eine prinzipielle Grenze durch die Trennung der Objekte 
gegeben, die im zweiten Fall nicht existiert. 


Um das Verfahren zu beschleumgen, muB man dem Segmentierungsalgorithmus im Programm 
von Bomfleth zunachst einen von Hand gesetzten SchweUwert vorgeben, mit dem man die 
interessierenden Objekte zur Segmentierung auswahlt. Ein weiterer Algorithmus bestimmt 
dann automatisch die Distanzen der nachsten Nachbam. 


IV.2 Observablen der Expert mentanalyse 


IV^.1 SchweUwert 

Um ein Maximum als interessierendes Maximum in dem in Abschnitt IV. 1.4 beschriebenen 
Bildverarbeitungsprogramm zu klassifizieren, wird ein SchweUwert gesetzt. Da bei der 
Setzung des SchweUwertes meist nur wenige Punkte auBer den interessierenden ubrigbleiben, 
steUt er eine GroBenabschatzung fiir den Hintergrund dar. Dies ermoghcht eine Aussage iiber 
die Qualitat der Praparation. 

IV.2.2 Brightness 

Die Bnghtness steUt ein MaB fur die Summe aUer emittierten Photonen liber aUe Voxel eines 
Signals dar. Sie gibt also Auskunft liber die absolute Emission eines Signals. 
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IV ^3 Signal volumen 

Das Signaivolumen ist ein MaB fur die Ausdehnung eines Signals. Dieses ist aber nur relaUv 
zum Bildverarbeitungsalgorithmus bestimmt. Das Volumen hSngt von der Ausdehnung des 
Klones, der Bindungsstarke Target-Probe, der fur die Detektion verwendeten Anzahl von 
Andkdrpem und AntikOrperlagen ab. Diese Effekte konnen auch glcichzcitig auftreten und 
eine Trennung nur schwer ermoglichen, ohne andere Observable zu berucksichtigen. 

IV .2.4 Integrierte Fluoreszenz Intensity (IF1) 

Die IFI ist das Produkt aus Brightness und Volumen eines Signals und macht sotnit eine 
Aussage iiber die emittierte Photonenzahl pro Volumeneinheit. Sie enthalt die Information 
iiber die Effizienz der Fluorochrome in Abhangigkeit der Preparation und ihrer lokalen 
Umgebung. 


IVii Distanzen zwischen genetischen Markern 

Um Abstande innerhalb eines Chromosoms zu messen, markiert man mindestens zwei Stellen 
auf einem Chromosom. Bestenfalls erhalt man vier Signale mit zwei rauralichen Abstanden 
auf je einem der beiden homologen Chromosomen. Lassen sich zwei Signale nicht mehr 
trennen, ist im allgemeinen das Volumen groBer, da die Signale dann 100% deckungsgleich 
sind. 

In einem Kern kommt immer ein kleinerer und ein groBerer Abstand vor. Man kann fur die 
beiden Abstande getrennte Verteilungen erstellen. Sofem die Mittelwerte signifikant 
verschieden sind, laBt dies Ruckschliisse auf die Struktur der homologen Chromosomen zu. 
Der signifikante Unterschied ist dabei nicht auf Sortiereffekte zuruckzufiihren [Ku96, Arif Esa 
-unveroffentlicht]: 

Wahlt man aus einer HaufigkeitsverteUung von Abstanden zufallig Abstandspaare aus, ist 
immer einer kleiner als der andere. Man konnte nun annehmen, daB bei der Sortierung der 
fn.^ gr °^ ere a Dlslanzen zw ei unterschiedliche HaufigkeitsverteUungen 

« PaarC iSt j6d0ch statistisch unabhangig, so daB es gleich 

zu wahlen. WcTein^olrh p inan ) derlie 8 endes Paar Oder ein weit auseinanderliegendes Paar 
Aufsnaltune einer H-i r ^ -t 8 *’ e ^e n ^* s gleichverteilt. Dies entspricht einer 
S7 Ciner HauP, g keitsv erteilung in zwei HaufigkeitsverteUungen mit demselben 


,V - 2 ' 6 “ S "‘ DZ zwi schen homologen Chromosomen 

hom2g^^omoS m a e n f ^ man den AbStand ^ 

den Markerschwerpunkten verwenden epai .. Man kann hierbei den Abstand zwischen 
Hieraus kann man feststellen ob Chr 8 W ° dlC Marker auf dem Chromosom liegen. 
zueinander besitzen Oder nicht. bfelSSST 1 flxierte raumliche Orientierung 
die GroBenverteilung des Zellkems ° 1 Un ® dcr d* 0, nosoraen ist auch ein MaB fur 
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IV3 Experimentdurchfuhrung 


Za Begin dieser Arbeit waren sowohl die MeBergebnisse von van den Engh [vdE 92], als auch 
diejemgen von Yokota rYo95] bekannt. Dabei wurden jedoch nomfalc Epifluoreszenz 
mdcroskope benutzt, d. h. nur die zweidimensionale Projektion des dreidimensionalen 
A^tandw wurden gemessen. Da dreidimensionale Abstande durch diese Projektion urn im 
Nlittel 2/3 reduziert werden, sinkt die Distanzauflosung vor allem bei einfarbiger Detektion. 
Zwex Markierungen wurden nur dann als einzelne Punkte wahrgenommen, wenn sie einen 
Mindestabstand besaBen, der von der SpotgroBe der Markierung abhangig war Die 
Praparationsmethode hatte dariiberhinaus den Nachteil, daB sie selbst stnikturzerstorend wirkte 
(Abs.1.6 und [Yo 95]). Die Aufnahmen der Zellkeme wurden auf handelsublichen photo- 
graphischen Film gespeichert und die Distanzen dadurch bestiramt, daB die Negative auf ein 
Digitalisierungsboard mit ca. 10.000 facher VergroBerung projiziert wurden. Von Hand 
wurden dann die Signalmittelpunkte bestimmt, deren Abstande von einera Computerprogramm 
berechnet wurden. Deshalb blieb unklar, wie sich raumliche Abstande fiir genomische 
Distanzen unter 500kbp verhalten wurden und welches Modell fiir ein Chromosomen- 
territorium nun eigentlich anwendbar sei. 


Urn zwischen dem RWGL- und dem MLS-Modell zu unterscheiden, und um realistische 
experimentelle Ergebnisse zu erhalten, muBten bei den Experimenten verschiedene Parameter 
variiert werden. 


■ Es wurden sowohl isotonische, als auch hypertonische Praparationen durch- 
gefuhrL 

■ Es wurde gezeigt, daB die verwendeten Klone nicht nur innerhalb des 
Chroraosomenterritoriums liegen. Die Klone selbst wurden dabei so klein wie 
moglich gewahlt 

“ Um Abstande messen zu konnen, die nahe oder unter der Auflosungsgrenze des 
Mikroskops lagen, wurden die Proben mit unterschiedlichen Fluoreszenz- 
farbstoffen detektiert. Bei Detektion desselben Targets mit zwei verschiedenen 
Fluoreszenzfarbstoffen wurde die chromatische Verschiebung und die 
Lokalisationsgenauigkeit der Proben gemessen. 

“ Zur Verstarkung des MeBsignals wurden ein- und zweilagige Antikorper- 
schichten zur Detektion benutzt, und die Abhangigkeit von Signalvolumen sowie 
die Integrierte Fluoreszenz Intensitat (Ih'l) wurden untersucht. 

• Es wurden geniigend Zellkerne gemessen, um statistisch aussagekraftige 
Abstandsverteilungen zu erhalten, die adaquat mit den Simulationen verglichen 
werden konnen. 


IV Experimente 


68 


die 


IV .3.1 Experimentansatz und Objekttragerpraparation 

Expenmentansata und [^!ff^ ti Mun C hen konnten unter Leitung von Irina Solovei 
te^Sloblktr^erLgestellt werden. Als Genregion fur die Abstandsmessungen wurde 

b,w. Angelroano-Region auf Chromosom ,5 gew M , Die 
Heiatellung der X-Klone und ihB Charakttrtaerang wurden freu™Uicherwe 1S e von Dr. Karin 
Biitig und Prof Dr Horsthonke air Verfiigung gestellL Fur die Auswahl der Genregion war 
die GroBe der Klone (Abb.IV.7) und der beschriebene medizuusche Hintergrund maBgebend. 
to fdgenden sollen die einzelnen Objektrager beschrieben werden (Tab.IV.2): 


Hypertonische (2D) Preparation: 

Hypertonische Praparationen wurden durchgefuhrt, indem auf 5 Objekttragern 7 
Hybridisierungen angesetzt wurden. Drei Objekttrager (2D279, 2D280 und 2D281) wurden 
mit Y48 und Y60 hybridisierL Damit sollten groBere raumliche Distanzen gemessen und ein 
AnschluB an globales Chromosomenverhalten erzeugt werden. Objekttrager 2D280 wurde 
weiterhin bei der Detektion mit einer doppelten Antikorperlage versehen, um eine 
Amplifizierung des MeBsignals zu untersuchen. Auf Objekttrager 2D280 wurde einer der 
Klone difluoreszent markiert, um den chromatischen Shift zu bestimmen und die 
Klonlokalisation eindeutig festzustellen. 2D278A+B wurden gemaB 2D280 erstellt, allerdings 
unter Verwendung von X-Klonen und einem Wechsel der Lokalisationsdetektion. 
Entsprechend wurde auch 2D277A+B, allerdings mit monofluoreszenter Detektion und dem 
Versuch der DAPI-Anfarbung zum besseren Auffinden eines Kerns angesetzt. Beide zuletzt 
genannten Objektrager wurden hauptsachlich auf das im Aufbau befindliche Wellenfeld- 
mikroskop abgestimmt (GroBe der A,-Klone). 


Isotonische (3D) Preparation: 

Zwei Objekttrager wurden mit Y48 und Y60 hybridisiert. Durch ihre GroBe sollten sie auf 
jeden Fall binden und detektierbar sein. Wahrend 3D288 monofluoreszente Detektion 
aufweist, wurde in 3D289 nicht nur eine Detektions-Amplifizierung vorgenommen, sondem 

dCr q_Ann VOn Chromosom 15 angefarbt. Als KontroUe wurden auf 3D290 
und 3D291 die Chromosomenarme von Chromosom 1 markiert. SchlieBlich wurden drei Paare 
von A,-Klonen prapanert, wobei die Markerpaare mit nur einem Farbstoff markiert wurden. 

16 ^ r ® e ^ n ^ sse von der Praparationsmethode, den jeweiligen Klonen, der 
Detektion und vom verwendeten Fluoreszenzfarbstoff ab. 

Hybridisierun^un^DefpIrt^ Prap ^ lion wur( * e die Fixierung aus Abb.IV.2 und die 
Fibrob^^^ew^l^diV^ 1 aUS Als illen wurden menschliche 

Dadurch wurde eneichL daR sM? 6 !? v wu chsen, bis sie ihn vollstandig bedeckten. 

wurde erreicht, daB srch alle ZeUen to gleichen Zellstadium (G 0 -Phase) befinden. 
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Zentromer 



Armjf^aJL 

23 



Prader-Willi Region (Paternales Chromosom) YAC 48 
Angelmann Region (Maternales Chromosom) YAC 60 


Klone und ihre Genomische GroGe 


1 IMbp 




X 48.17 
17.8 


X 48.14 A. 48.7 




13.0 


$ 9,2 


—79- 

J2 5J . 




19.0 


[kbp] 


1 9 4, 9 


144.7 


213.9 


-J 


Genomische Distanzen zwischen Probenpaaren 


Abb.IV.7: 

Die verwendeten Klone, ihre GroBe und Lokalisation auf Chromosom 15. 


IV.3.2 Durchfuhrung der Mikroskopaufnahmen 

Um Alterungseffekte und Beschadigungen der Objekttrager zu vermeiden, wurde nach der 
Durchfuhrung der Preparation sofort mit der Aufnahme von Zellkemen begonnen. Zunachst 
soil auf die gewahlten Einstellungen des Mikroskops eingegangen werden, bevor die 
Ergebnisse der Durchfuhrung genauer beschrieben werden: 

Die Aufnahmen wurden im XY-Modus (Scan-Ebene in der Horizontalen) vorgenommen, da 
hierfiir die entsprechenden Auflosungen optimal geeignet sind (lateral ca. 250 nm, Axial ca. 
500 nm). Achtmal wurde iiber dieselbe Bildebene gemittelt, um Rauscheffekte zu 

unterdriicken. . 

Als Bildformat stellte sich bei vierfachem Zoom eine 512*512 Pixel-Matrix als sinnvo 
heraus. Ein Bildstapel enthielt zwischen 8 und 25 Bildebenen. In Anbetracht der achtmaligen 
Mittelung ergaben sich hiermit keine allzu starken Ausbleicheffekte der Fluorochrome. Einmal 
aufgenommene Kerne hatten oft ein zweites Mai aufgenommen werden konnen, ohne mit 
einem nicht aufgenommenen Kem verwechselt zu werden. Die VoxelgroBe ergab sic in 
lateraler Richtung zu 78nm und in axialer Richtung zu 200nm. 

Die Spannungsversorgung der Photomultiplier wurde durch Probieren ermittelt. Sie betrug 
V im griinen Kanal und 732 V im roten Kanal. Die konfokale Blende wurde auf den Wert 51 
gesetzt, als Objektiv wurde ein 63* Leitz-Objektiv mit NA=1,4 verwendet. Die Laser 
fill die einzelnen Farbkanale wurde fur den roten auf 14% und fur den griinen 
der Gesamtlaserleistung eingestellt. Alle Bilder wurden dabei im GlowOverUnder-Modus 
aul genommen, um etwaige Uberbelichtungen feststellen zu konnen. 
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Objekttrager 

Gelabelte Hone 

genomischer 
Ahstand fMbp] 

Detektion 

Ziel 

2D 277 A 

A48.1-Dig 

17-Die 

0.069 

SaDig-FITC 

Hypertonie 

Wellenfelri ^ 

2D 277 B 

A48. 1-Dig 
148.1 7-Die 

0.069 

SaDig-FITC 

+DAPI 

Hypertonie, DAPI, 

Wellenfelri 

2D 278 A 

148. 1 -Dig + A48. 1-Bio 
148.17-Die 

0,069 

Av-CY3 

SaDig-FITC 

Klonlokalisation, 

Chrom. Shift 

2D278 B 

148.1-Bio 

148. 17-Bio + A48. 17-Dig 

0,069 

Av-CY3 

SaDig-FITC 

Klonlokalisation, 
Chrom. Shift 

2D 279 

Y48-Dig 

Y60-Die 

1,0 

SaDig-FITC 

Monolage 

2D 280 

Y48-Dig 

Y60-Dig 

1,0 

MaDig + 
SaM-CY3 

Amplifizierung 

2D 281 

Y48-Bio 

Y60-Bio + Y60-Dig 

1.0 

Av-CY3 

SaDig-FITC 

Klonlokalisation 

3D 282 

148.1-Dig 

148.14-Die 

0,018 

SaDig-FITC 

Monolage 

3D 283 

148.1-Dig 

148.14-Dig 

0,018 

MaDig + 
SaM-CY3 

Amplifizierung 

3D 284 

148.1-Dig 

148.17-Die 

0,060 

SaDig-FITC 

Monolage 

3D 285 

148.1-Dig 

148.17-Dig 

0,060 

MaDig + 
SaM-CY3 

Amplifizierung 

3D 286 

148.7-Dig 

148.17-Die 

0,194 

SaDig-FITC 

Monolage 

3D 287 

148.7-Dig 

148.17-Dig 

0,194 

MaDig + 
SaM-CY3 

Amplifizierung 

3D 288 

Y48-Dig 

Y60-Dig 

1,0 

SaDig-FITC 

Hypertonie 

3D 289 

Y48-Dig 

Y60-Dig 

15q-Bio 

1,0 

MaDig + 
SaM-CY3 
Av-FTTC 

Klonlokalisation 

auf 

3D 290 

lp-Bio 

lq-Dig _ 

- 

Av-CY3 

SaDie-FTTC 

Armpainting 3D 
KYmtrolle 

3D 291 
Tab.IV.2: 

lp-Bio 

lq-Dig 

~ ' 

Ethanol Postfix 
Av-CY3 
SaDig-FITC 

Ethanol Postfix 
Kontrolle 


Ej^rim^Unsate and hergestelhe Objekttrager. 

MaDig: Maus Antikorper gegen Digoxigenin, Av-CYV a m- 

Schafantikorper gegen den NUusantikorn^t ' S** “ ,t CT3 ““**«<. SaM-Cy3: 

Digoxigenin, mit FITC markiert. * J markiert, SoDig-FlTC: Schafantikorper gegen 


Der GlowOverUnder-Modus ist pine • r. 

ein Farbwert auf einer Skala von 0^fc 8 o« arSteUun8 eines Bil dpunktes, die anzeigt, wan* 
uberschreitet. Die Aufnahmen warden • 256 Cme Helli S k eit von 0 unter- oder von 25 6 
nmpAwanrielt n ln euie Glowdarstellung fur die Bildverarbeitung 
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Abb.IV.8: 


J^hnitteeri^iirc^ine^Cerr^oi^Objektd*age^D287^enoinischei^bstand^ 18kbp 


^ X 

* 


\ 


\ 

V 

V 

^ ' 

jf 


4 




Abb.IV.9: 

Projektion der Schnittserie von 
Abb.111.8 in eine Ebene (Extended 
Focus) nach der Medianfilterung. 

Folgende Signale wurden 
registriert: 


1: nur ein Spot. 
2: zwei Spots, d = 


135nm. 


3: zwei Spots, d = 182nm. 
4: zwei Spots, d = 424nm. 

5: zwei Spots, d = 469nm. 
6: nur ein Spot. 






72 


IV Experimente 


« Hie obise Skala iiberschreitende Werte auf 0 oder 256 gesetzt 
Glowdarstellung heiBt, dajj aic ! ’ t esteUtt so wurde die Laserleistung reduziert und ein 

werden. Wurden Uberbehchtung 6 ions6 , wur den diejenigen der Firmen Neolab und 
anderer Kem aufgenommen. Ais m 

Zeiss verwendet. 7un achst auf ihre Taugbchkeit hin durch ein normales Okular 

Generell wurden die Ke Gesicht sfeldes positioniert. Kerne wurden als tauglich 

ausgewahlt nnd un Zentmn $ich vermuten lie6 , d aB man durch die Aufnahme 

ael s' 1 ”" 51 ' Dan,ch wurde Hilfe 

oder die Bildver g ^ das zgntrum des ScannerbUdes verschoben. Dabei 

eI S °BuSb M e im Scanbetrieb ohne Micelung und bd « _3% Use,l« sl „ ng 
aufgenommen (nicht im Dauerbemeb!). Der Kem zeigte kemeriei Schadigungen ode, 
friihzeitiges Ausbleichen der Fluorochrome. 


IV.4 Ergebnisse der Experimentanalyse 


Wie Tab.IV.4 zeigt, wurden insgesamt 1303 Bildstapel aufgenommen, die insgesamt 1661 
Zellkeme enthielten. 166 Aufnahmen wurden von J. R., 1137 Aufnahraen von T. A. K. 
gemacht. Insgesamt wurden 212 Stunden fur die gesamten Aufnahmen benotigt. Dies 
entsprach 6,2 Aufnahmen oder 7,8 Kemen pro Stunde. Rechnet man ca. 1,5 Stunden fur die 
jeweilige Inbetriebnahme und das Brennen der CD ab, so erhalt man mit 7,0 Aufnahmen oder 
9,0 Kerne pro Stunde als MaB fur die Leistungsfahigkeit des Mikroskopes. Die gesamte 
aufgenommene Informationsmenge betrug fur aUe Aufnahmen 7,7GByte. 



Protokol) der AufnahiDen der 


prapanerten Objekttrsiger. 
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FV.4.1 Durchfuhnmg der Experimentanalyse 

Nachdem die Keme wie oben beschrieben aufgenommen waren, wurden die Bilder in das fiir 
die Bildverar^itung benotigte Rohformat konvertiert, wobei auch die Median-Filterung 
erfo gte- a e Bildverart>eitungsalgorithmen von Bomfleth nur fiir kleine Spots geeignet 
sind, wurden von den Armproben bisher nur diejenigen ins Rohformat konvertiert, in denen 
sich noch einzelne Spots befinden. Danach wurde der Schweliwert fiir jede Aufnahrae und 
jeden Farbkanal erner Aufnahme einzeln von Hand bestimmt 

Bisher wurde dieses Verfahren auf die Objekttrager 2D278A+B, 2D279, 2D280, 2D281, 
3D283, 3D285, 3D287 und 3D289 angewendet. Dabei erfolgten ca. 1200 Schwell- 
wertfestlegungen.Die benotigten Rechungen wurden auf der Silicon Graphics Iris der 
Arbeitsgruppe Angewandte Optik und Informationsverarbeitung durchgefiihrt. 


IV.4.2 Analyse der Objekttrager: 

Bisher wurden die Objekttrager 2D279 und 2D280 vollstandig ausgewertet. Die Objekttrager 
2D278 A und B wurden bis auf Distanzmessungen ausgewertet Die Distanzen miissen hierbei 
noch mit der chromatischen Verschiebung korrigiert werden. Fiir die Objekttrager 3D283 und 
3D285 liegen zwar die Daten schon vor, es konnten bisher aber nur die Distanzen analysiert 
werden. Leider lassen sich daher keine Aussagen iiber systematische Effekte bei Brightness, 
Volumen und IFI ziehen. Auch bei den Distanzen konnen Aussagen, die mit diesen GroBen 
zusammenhangen, nur schwer gemacht werden. Bei Vorlage aller Daten des Screenings sollte 
dies jedoch moglich sein. Deshalb sollen in dieser Arbeit nur die Ergebnisse der Single- und 
Double-Spots (Nomenklatur siehe unten) besprochen werden. Die Ergebnisse von SP1 bis SP4 
sollen dagegen nur in den Tabellen aufgezeigt werden. 

Bevor die Ergebnisse besprochen werden, sollen einige Begrifflichkeiten erlautert werden: 

In einem Kern konnen maximal 4 Signale registriert werden, wenn die Hybridisierung eines 
Markerpaares gelungen ist Jeweils zwei Signale befinden sich dabei auf einem der homologen 
Chromosomen. Da dieser Optimalfall jedoch nicht immer eintritt, wurden die Keme in 5 
Klassen unterteilt: A = 2*2 Signale, B = 1*2 Signale und 1*1 Signal, C = 1* 2 Signale, D = 
2*1 Signale und E = 1*1 Signal. 

Wahrend im Fall A und C immer eindeutig ist, daB ein Signal von einem Target herriihrt, laBt 
sich dies im Falle von B nur fiir das eine Chromosom sagen. Fiir die einzelnen Signale aus 
Klasse B, D, E ist jedoch nicht klar, inwieweit ein Signal von einem oder von zwei nicht 
trennbaren Targets stammt. Diese Spots sollten also im Mittel groBer sein, da sie Anteile aus 
beiden Gruppen besitzen. 

Fiir das Analyseprotokoll wurden die Signale nach der Zugehorigkeit zu einem Chromosom 
und ihrer GroBe nach sortiert Der Signale mit der groBten Brightness wurde dabei immer als 
SP1 bezeichnet, das dazugehorende zweite Signal als SP2. Ebenso wurde mit dem anderen 
Chromosom verfahren. Grundsatzlich gehoren SP1 und SP2 zu einem und SP3 und SP4 zum 
anderen Chromosom. Aus der daraus folgenden Auflistung im Analyseprotokoll laBt sich im 
Zusammenhang mit der Klassifizierung immer erkennen, welche Eigenschaften ein Spot einer 
Aufnahme besitzt. 

Spot SP1 besteht aus alien Klassen; SP2 nur aus Klasse A, B und C; SP3 nur aus Klasse A, B 
und D und SP4 nur aus Klasse A. 
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. mm ,r in Distanz 1 und Distanz 2 unterteilt, wobei Distanz 1 immer d - 
Die Distanzen wurden immer in uie 

groBere ist. tv 5 bis IV. 10 sind sowohl die Mittelwerte der Observablen 

Its s5»SpSsP 4, als auch die Mittelwerte dereindeutig einzelnen (single), ,| s auch 

der nicht dndeutigen einzelnen (double) Spots aufgefnhrt. 


IV.4.3 Schwellwert 

Der Schwellwert liegt zwischen 56 und 1 16,5 (Tab.IV.4). Werte von 56 bis 67 scheinen dabei 
vorherrschend zu sein. Eine weitergehende generelle Aussage erscheint jedoch bisher nicht 
moglich. 



Schwellwert 

IF 

Sin 

I 

ele 

IFI 

Double 


0 

± 

a 

N 

0 

± 

a 

N 

0 

± 

a 

N 

2D278Ag 

86,8 

6,3 

20,8 

11 

170 

14,7 

55,0 

14 

269 

44,2 

132 

11 

2D278Ar 

66,3 

5,0 

16,6 

11 

101 

11,6 

46,3 

16 

- 

- 

- 

_ 

2D278Bg 

56,0 

1,3 

11 ? 3 

71 

137 

5,4 

59,8 

122 

- 

_ 

_ 

_ 

2D278Br 

67,7 

1,8 

15 ? 0 

71 

118 

5,5 

47,7 

73 

188 

9,1 

85,1 

88 

2D279 

116,5 

JA. 

81,2 

129 

529 

21,1 

340 

260 

698 

54,0 

458 

72 

2D280 

57,4 

_LL 

18,1 

143 

249 

8,5 

154 

320 

318 

17,9 

152 

72 


Tab.IV.4: 

Ergebnisse der Schwellwert- und IFI-Anlyse verschiedener Objekttrager und Signale. 



2D280 
Tab.IV.5: 

Egebuisse der Brightness und des Volumens 


verschiedener Objekttrager und Signale 
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Brigh 

SI 

tnes 

>1 

S 

Brightness 

SP2 

Brightness 

SP3 

Brightness 

SPd 


0 

± 

a 

N 

0 

± 

a 

N 

0 

± 

a 

N 

0 

± 

(j 

N 

2D278Ag 

19,9 

4,6 

15,5 

11 

27 8 

6,2 

15,2 

6 

16,8 

3,3 

7,3 

5 





2D278Ar 

14,1 

1,9 

6,0 

10 

- 

- 

_ 

_ 

17,2 

5.3 

10,5 

4 





2D278Bg 

19,7 

1.6 

13,3 

67 

- 

- 

. 

_ 

14,6 

1,6 

11,3 

53 





2D278Br 

23,2 

1,7 

14,2 

69 

23,5 

2,9 

23,5 

33 

22,5 

1,9 

14,6 

58 

25,3 

3,6 

10,9 

9 

2D279 

14,5 

1,1 

13,1 

129 

7,2 

0,8 

7,8 

105 

10,2 

0,7 

6,9 

97 

11,7 

2,7 

1,7 

43 

2D280 

ToK TV 

6,2 

0,7 

8,5 

143 

2,6 

0,78 

3,1 

116 

5,0 

0,6 

5,9 

94 

2,5 

2,4 

0,34 

48 


Ergebnisse der Brightnessanalyse verschiedcner Objektragcr und Signale. 




Volui 

SP 

men 

1 


Volumen 

SP2 

Volumen 

SP3 

Volumen 

SP4 


0 

± 

a 

N 

0 

± 

a 

N 

0 

± 

a 

N 

0 

± 

a 

N 

2D278Ag 

22,5 

9,9 

32,9 

11 

15,5 

5,7 

14.1 

6 

20.0 

8,9 

19,8 

5 





2D278Ar 

8,00 

1,2 

3.65 

10 

- 

. 

. 

_ 

5,75 

1.1 

2,2 

4 





2D278Bg 

11.8 

2,0 

16,2 

67 

- 

. 

_ 

. 

14,3 

1,6 

11,8 

53 





2D278Br 

9,2 

1,7 

13,8 

69 

8,7 

1,2 

7,1 

34 

16,4 

3,2 

24,8 

58 

7,9 

2,0 

6.1 

9 

2D279 

70,1 

6,3 

71,6 

129 

101,2 

6,1 

63.3 

105 

89,1 

7.9 

76,6 

97 

81,8 

8,3 

52,7 

43 

2D280 

97 

6,2 

75 

143 

127 

7,0 

75 

116 

117 

7,8 

76 

94 

123 

7,5 

52 

7.56 


Tab.IV.7: 

Ergebnisse der Volumenanalyse verschiedener Objektrager und Signale. 


IV.4.4 Brightness 

Die Brightness fur ein Single-Signal betragt zwischen 1 1,4 und 21,5 und fur ein Double-Signal 
zwischen 4,3 und 17,9 (Tab.IV.5). Sie ist in fast alien Fallen kleiner als fur ein Single-Signal. 
Auffallig ist, daB die Brightness dabei fur alle Objekttrager keine definitiven Unterschiede 
zwischen den YAC und Lambda-Klonen zeigt. Eine Erklarung kann hier nur die gleiche 
Mikroskopeinstellung und die gleiche Laserleistung, mit der die Aufnahmen gemacht wurden, 
liefem. 

rv.5.5 Volumen 

Beim Volumen zeigen sich groBe Unterschiede zwischen den Objekttragem. 2D278 A und B 
haben ein Volumen zwischen 7,8 und 15 fur Single- und zwischen 14,9 und 19,3 fiir Double- 
Signale. 2D279 zeigt Werte von 80,6 und 67,2 und 2D280 zeigt Werte von 108 und 137. 
Hieraus ergibt sich, daB die YAC-Klone ein wesentlich groBeres Volumen als die Lambda- 
Klone besitzen. Dies korreliert mit der GroBe dieser Klone. Der Unterschied zwischen 2D279 
und 2D280 ist auf die Amplifizierung des MeBsignals mit zwei Antikorperlagen 
zuriickzufiihren. Fiir ein Single-Signal betragt sie 35% und fiir ein Double-Signal fast 100%. 
Der Zuwachs von 100% diirfte darauf zuriickzufiihren sein, daB hier der Anted von Spots, die 
aus zwei nicht mehr trennbaren Klonen bestehen, wesentlich hoher ist. Dies wird durch 
Zahl der gemesssenen kleineren Distanzen 2 gestiitzt: Bei 2D279 sind dies 89 und bei - 
48. also praktisch nur noch die Halfte. Generell gilt fiir alle Objekttrager, daB das Volumen r 
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5,8 
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- 
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12,6 

4 
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137 

7,7 

63,3 

67 

- 

- 

- 

- 

140 

7,6 

55,4 

53 

- 

- 
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2D278Br 

147 

9,5 

79,5 

69 

137 

9,3 

53,5 

33 

179 

11,2 

86,0 

58 

155 

26.0 

78,2 

9 

2D279 

454 

26,2 

297 

129 

503 

31,5 

322 

105 

581 

43,4 

421 

97 

444 

202 

32 

43 

2D280 

269 

12,8 

153 

143 

230 

9,9 

106 

116 

302 

21,4 

207 

94 

241 

130 

18 

48 


Tab.IV.8: 


Ergebnisse der IFIsanalyse verschiedener Objektrager und Signalc. 
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2D279 

694 

33,0 

395 

143 
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40,7 

415 
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43,6 
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89 

8558 

489 

3175 

42 

2D280 

848 

30,5 

389 

162 

922 

39,1 

418 

114 

673 

33,5 

232 

48 

10226 

517 

4602 

79 

3D283 

358 

33,6 

349 

99 

410 

41.9 

368 

77 

177 

37,1 

174 

22 





3D285 

Tnk n; n 

561 

41.2 

387 

95 

586 

43,8 

382 

76 

461 

90,4 

394 

19 

. 
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Ergebnisse der Abstandsanalyse verschiedener Objektrager und Signale. 

iSt ’ S ° daB diC ° ben gemachte v ermutung aus den Signalklassen 


IV.4.6 Integrierte Fluoreszenz-Intensitat (IFI) 

(Tab.in.8). Fur 2D278 A und B P r i n zipiell dasselbe wie fur das Volun 

bis 269 fur ein Double-Signal. Fur 20970 1 V ° n 101 bis 137 fiir ein Sin £ ,e ' und von 1 

und 318. Die niedrigeren Werte fiir die Werte 529 und 698 und fur 2D280 2 

werden, wobei eine verstandlichere a ° nnen nur niit der niedrigeren Brightness erkl 
Auch fur das IFI sind die Single- Werfe^i S ^ stemat * sc he Erklarung bisher nicht moglich i 
groBer gegenuber den Lambda-Klonen ^ S ^ ou ^ e- Werte bzw. fiir die YAC-Klc 

1VA7 Dislanzen ***« g. ntUsche|| 

Uer genomische Abstand fur 2D97Q . 

™“! e £ “ SU " 2 '»« 694±33„ m , 2Mbp F “ r 20219 wmde eiK ,0 ° 

Disianz ?92 2 'uo ' F “' 20280 »«n ZTL S™ 758±41 “»« Disia" 2 2 * 

Disianz betraet <!, ?" d fUr D «stanz 2 673+ 34 die ‘ otale mittlere Distanz 848±31 nm, ft 

° m,t nm . obwohl dieselbe ^ er ^ nter *>chied in der totalen mittlere 

ben Kl °ne verwendet warden. Der Abstand vo 
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Jr totalen Distanz zu Distanzl und Distanz2 ist fur beide Objektrager ungefahr gleich Alle 

Werte sind sigmfikam unterschiedlich und wurden alle auf mindestens 4 Arten mil IMP 
getestet. 

^ rkJ . arU " g h . egt c in d ® n unter schiedlichen Volumina der beiden Objektrager, die sich um 
3 ~ c furdl ® S ^ g e ‘ S ' gnale “ nterscheiden Ein Signal wachst dadurch in seinem Durchmesser 
um ca. 18%. Durch die VolumenvergroBerung lassen sich die Signale jedoch auch weniger 
leicht trennen. Da man beim Auswahlen der Kerne aber darauf bedacht ist, immer vier Signale 
zu trennen, wahlt man bei 2d280 automatisch Signale mit groBerem Abstand aus. Im CLSM 
konnte der Durchmesser eines Signals mit Hilfe der Funktion „Intensitats-ProfiT auf ca 450 
nm bestimmt werden. 

Eine VergroBerung entspricht also einem Anwachsen des Durchmessers auf ca 530 nm Fur 
den Abstand zweier Signale bedeutet dies, daB er ca. 80 nm groBer ist was in der 
GrbBenordnung des beobachteten Effektes liegt. Eine genauere Analyse kann’aufgrund der 
derzeitigen Datenlage leider noch nicht erfolgen. 

Fiir 3D283 konnte eine totale mittlere Distanz von 358±34 nm, eine Distanz 1 von 410±42 nm 
und eine Distanz von 177±37 nm fur einen genomischen Abstand von 194kbp bestimmt 
werden. 


Fur 3D285 und eine genomische Distanz von 69kbp betrug der Wert fur die totale mittlere 
Distanz 561±41 nm, fur Distanz 1 586±44 nm und fur Distanz 2 461±90 nm. 

Diese Werte mUssen jedoch aus zwei Grunden mit Vorsicht behandelt werden: Erstens waren 
durch die iiberbelichteten Aufnahmen viele Kerne nach der Segmentierung fragmentiert, so 
daB Fehler bei der Distanzbestimmung gemacht werden konnten. In welcher GroBenordnung 
dieser liegt, kann bisher allerdings nicht abgeschatzt werden. Zweitens waren vor allem bei der 


Totale Distanverteilung 



_ Test auf Unterschiedlich 

Verteilung sortierter D.stanzen keit der Mitt elwerte 



Haufigkeit 

# Distanzen: 305 

Mittelwert: (776±23)nm 

Std. Abweichung: 399nm 

Shapiro Wilk Test auf Normalitat: W=0.87 


# Distanzen: 218 # Distanzen: 87 

Mittel: (844±26)nm Mittel: (606±41)nm 


Faktor fOr Signifikanten 
Unterschied > 0 


Abb.IV. 8: 

Distanzverteiiung und Signifikanztest: 

Link^>: Uaufigkeits verteilung alter Distanzen von Objekttrager 2D279 und 2D280 

Rechts: Lm linken Bild sind die Sortierergebnisse in Distanz 1 und 2 dargestellt. Im rechten Bild zeigen die 
beiden sich nicht beriihrenden Rreise, dab der Unterschied der Mittelwerte signifikant verschieden ist. Es 
wurdc der HSU-MCB Test verwendet [SAS 96]. 
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„ . wpniee Distanzen auszuwerten bzw. genereU viele Signale nicht mehr „ 

Distanz 2 nur noc g die ser Distanzen nur die groBeren Distanzen erfaBt in 

SSjSSr £*£ daB fur eine groBere genomtahe Dis,anz eta Heine, e 

Standaniabweichung ft alle totaen Distazen zwischen 350 
und 400 nm liegt. SoUte es sich heraussteUen, daB die Standardabweichung auch bei weiteren 
Messungen in diesem Bereich liegt, so konnte dies auf einen EinfluB der Fluoreszenz-in-situ- 
Hybridisierung hindeuten, da eine solche Einformigkeit nicht zu erwarten ware. Ein 
signifikanter Unterschied zwischen den Distanzen bei verschiedenen Experimentatoren konnte 


nicht festgestellt werden. 

Eine Analyse ergab, daB die von Rauch aufgenommenen Zellkeme generell etwas kleinere 
Mittelwerte in den Distanzen zur Folge hatten. Dies ist darauf zuriickzufiihren, daB T. A. K 
hauptsachlich als erster Aufnahmen auf einem Objektrager vornahra. Es ist deshalb 
wahrscheinlich, daB der Autor besonders augenfallige Kerne mit gut sichtbaren 4 Signalen 
aufnahm, die Rauch nicht mehr zur Verfugung standen. 

Alle hier gezeigten Abstande sind signifikant verschieden, so daB damit geraachte Aussagen 
statistisch gesichert sind. Stellvertretend soli eine Distanzverteilung gezeigt werden, die durch 
Mittelung von Objekttrager 2D279 und 2D280 erhalten wurde (Abb.IV. 10). 


IV .4.8 Distanzen zwischen homologen Chromosomen 

Die Distanzen liegen zwischen 8,5 pm und 10,2+0,5 pm, was mit den Erwartungen und einer 
ZellkemgroBe von 20pm*12pm*6pm ubereinstimmt. 
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V Vergleich zwischen Simulation und Experiment 


Die Simulationsdaten entsprechen Abb. HI. 24, d.h. alle Simulationskurven staramen von 
Simulationen mit niedriger Excluded- Volume-Wechselwirkung. In Abb.IV.l wird zunachst ein 
Vergleich zwischen den Simulationen und den Daten von van den Engh [vdE92], Yokota [Yo 
92], Monier [unpublished], und den un experimentellen Teil dieser Arbeit gewonnen Daten 
dargestellt 

Die Daten von van den Engh zeigen ein auBerst zwiespaltiges Verhalten. Wahrend sie 
unterhalb von 500kbp zu einem MLS-Modell passen konnten, sind sie oberhalb von 500kbp 
mit verschiedenen RWGL-Modellen vereinbar. Oberhalb von IMbp entsprechen die Daten 
nicht einmal einem RWGL-Modell mit 5Mbp groBen Schleifen. Dies ist wahrscheinlich auf 
die von van den Engh benutzte und sehr strukturzerstorende Praparation zunickzufiihren. 

Die Daten von Yokota, die mit einer istotonischen Praparation gewonnen wurden, passen 
dagegen sehr gut zu einem MLS-Modell mit 126kbp groBen Schleifen mit einem Linker von 
1200nm. Hier ist zu beriicksichtigen, das die Daten von Yokota 2D-Aufnahmen sind, also um 
ca. 2/3 zu klein sind und um diesen Faktor korrigiert werden mussen. 

Die Daten von Monier passen sehr gut zu einem MLS-Modell. Aufgrund der noch zu geringen 
Auflosung laBt sich jedoch noch nicht sagen, welchen Parametern die MeBwerte genau 
entsprechen. Die Daten von Monier stellen bisher jedoch die besten und in ihrer Gesamtheit 
eindeutigsten Werte dar. 

Die Werte des experimentellen Teils dieser Arbeit zeigen dagegen ein uneinheitliches 
Verhalten: bei 2D279 spricht die totale mittlere Distanz fiir ein MLS-Modell von 126kbp- 
Schleifen und eine Linkerlange von 1200nm. Die Distanzaufteilung in groBere und kleinere 
Distanz ergibt fiir die groBere Distanz 1 ein MLS-Modell mit einer Linkerlange von 1800 nm 
und fiir die kleinere Distanz 2 ein MLS-Modell mit einer Linkerlange von 600 nm. Fiir 
Objekttrager 2D280 ergibt sich ein MLS-Modell von 1800 nm fiir die totale Distanz, etwas 
mehr als 1800 nm fiir Distanz 1 und 1200 nm fiir Distanz 2. Fiir die totale Distanz und Distanz 
1 wiirde jedoch auch ein RWGL-Modell mit SchleifengroBen zwischen 126kbp und 503kbp 
die Daten erklaren. Aufgrund der in Abschnitt III.3.C gemachten Aussagen ist jedoch 
anzunehmen, daB 2D279 die korrekteren Werte darstellt. 3D283 deutet zwar auf ein MLS- 
Modell, allerdings kann wegen der zu geringen Auflosung der Daten nicht gesagt werden fiir 
welches. Fiir 3D285 liegen die Werte noch generell zu hoch um irgend eine Aussage iiber ein 
Modell treffen zu konnen. 

In Abb.IV.2 werden dieselben Simulationen mit neueren Daten von Yokota [Yo 97] 
verglichen. Yokota wendete dabei zwei Praparationsmethoden an, die sich von seinen friiheren 
[Yo 95] unterscheiden. Mit MAA bezeichnet er eine Fixierung, bei der die DNA bei 70°C mit 
70% Formamid denaturiert und danach in 70% Ethanol dehydriert wird. Mit PFA fixierte 
Zellkeme werden bei 90°C unter denselben Bedingungen denaturiert und dehydriert. Diese 
Fixierungen wurden auf Fibroblastenzellen und HeLa-Zellen angewendet 
Folgende Genregionen wurden in Fribroblastenzellen mit der MAA-Fixierung prapariert: 
4pl6,3, 6p21,3, 21q22.2, Xq28, und Xp21.3. In Fibroblastenzellen mit PFA-Fixierung wurden 
Xq21.3 und Xq28, und in HeLa-Zellen mit MAA-Fixierung wurden Xq28 und Xp21.3 
prapanert. Von der MAA-Fixierung in Fibroblastenzellen ergibt der Vergleich mit den 
Simulationen fiir 4pl6.3 ein RWGL-Modell mit 5Mbp- oder 4Mbp-Schleifen. Fur 6p21.3 
ergibt sich ein RWGL-Modell mit 3Mbp- oder 2Mbp- Schleifen. Fiir 21q.2 und Xp21.3 ergibt 



, simulation und Experiment 




RWGW-Modell: 
♦Mbp 
3Mbp >Mbp 
IMbp 

SSlkbp LKkbp 

dunneLinien 


MLS-MbddI: 


Linkerlang^ 

1200 nm 
ISOOnm 
2400nm 
did* Linien 


NewLiiwar: 
dictoe schwarae Liiue 


O vaiui*jiEngh[va^2]: 
■ Ycteota (Yo75j! 

• Karin Mcni er: 

Q PWS: trial* Difttanz 
A WS: Diftaal tuu12 


Abb.V.l: £ 3 D 280 : Diftanzl uwi J 

v ergleich der Simulationen mit Dat 

dieser Arbeit gemessenen Daten. V ° n Va “ de " Engh [vdE 92 L Yokota [Yo95], Monier und den in 


fVni • ullLer ^alb von lMbn Pn aIa e * n RWGL-Modell mit einer Schleifengi 
der F^y 1 JCWeilS diesel ^ en Ergebnisse V)' a '^ xierter HaLa-Zellkerne ergeben sich fur 
V° r ali^^ Und dem hin 1CS eutet au ^ e ^ ne Unabhangigkeit der Ergebnisse vc 

« jedoch die Auflosung der Simulade* 


^ cm unemheitliches BilH h* 

entsprieht. Fiir die PFA fixierten FK^/^^'^odell mit Schleifen von ungefahr 1Mb] 
Modell mit 1800nm Linkerlan He un ^ 0bl ^ stenzellk erne ergibt sich fiir Xp21.3 ein MLS 
on eiwas unterhalb 8 d fur Xq28 ein RWGL-Modell mit einer SchleifengroB. 

Genloci jeweils Hiao<»ii „ 1Dp ‘ Fur MAA-fixiert^r Her „ sich fur dit 

; vor 



raiimliche Distanz [jim] 


V Vergleich zwischen Simulation und Experiment 


81 



RWGW-M odell: 
§Nbp 4Mbp 
3Mbp 2Mbp 
IMbp 50 35? bp 
2Slkbp l^kbp 
dimne Litiien 


MLS- Mod 41: 
Linkerlangc 

1200 nm 
18 00 run 
2400 nm 
dlckeLmien 


Yckota [Yc07]: 

4pl63 ♦ MAA-Xp21.3 

--4-- 6p<213 W MAA-Xq28 

2Iq22.2 y PFA-Xp2 L3 

Xp2L3 B PFArXq28 

Xq28 


NewLirueai: 
didae schwaize Lime 

Abb.V.2: 

Vergieich der Simulationen mit Daten von Yokota [Yo 95]. 
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VI 

Entspricht nun die dreidimensionale Organisation von Chromosomen in der Interphase 
RWGL- oder dem MLS-Modell ? 

Eine abschlieGende Antwort auf diese Frage kann in dieser Arbeit leider noch nicht geg e b en 
werden. Experimentelle Daten konnen zwar durch die in dieser Arbeit durchgefiihrten 
Simulationen interpretiert werden, ergeben jedoch kein eindeutiges Bild. In Teilaspekten 
konnten jedoch klare Aussagen und Fortschritte gemacht werden: 

Ausgehend von der Tatsache, daB Chromosomen Territorien bilden, wurden verschiedene 
Modelle der dreidimensionalen Organisation simuliert. Da nicht bekannt war, ob und unter 
welchen Bedingungen sich die Modelle unterscheiden lassen, wurde ein Feld von 
Parameterabhangigkeiten aufgestellt. Dieses lieferte zahlreiche Abhangigkeiten, deren 
wichtigste hier kurz erwahnt werden sollen: 

- Bei einer VergroBerung der Linkerldnge und kleiner EV-WW wuchsen dreidimensionale 
Abstande, so daB das Chromosom weniger dicht wurde, wie aus der Polymerphysik 
plausibel ist. 

Bei VergroBerung der Schleifen bei kleiner EV-WW wuchsen ebenfalls die 
dreidimensionalen Abstande, und das Chromosom wurde weniger dicht. Im RWGL- 
Modell ist dies durch die groGen Schleifen bedingt. Im MLS-Modell kann dies nur 
durch eine Wechselwirkung der Rosetten erklart werden. 

Bei hoher EV-WW waren die Chromosomen durchweg schlecht relaxiert und die 
me( * r '^ e u 3 * S ^ S er inger EV-WW, da Umlagerungen von Teilen des 

- SS. 5 , 

definierter Oh fi - I, et f, n . S * C ^ ^ urcb den hohen Potentialtopf kompakte Gebilde mit 
Oberflache de-fri^ 6 C n ‘ ec * r *S er EV-WW kann man jedoch nicht von einer 
Innem, sprechen r ° mosoms ’ au fgrund der hohen Dynamik und Zuganglichkeit zum 

zeigte, in einem RanHi^ r0 T OSOrnS s * cb > w ’ e die Aufenthaltswahrscheinlichkeit 

m R “ db "'^ von ca, 1 „ m auf „„d s.arke Dynamik. 

groGen Schleifen im RWrf'vff 11 * m MLS-Modell distinkte Einheiten, wahrend die 
m RWGL -Modell sich durchdringen konnen. 

RWri iSChen Vorher s agelfl on S^chl’iih 36 diC Simulati on eines RWGL-Modells mit den 

RWGL- und MLS-Modell ze u ub eremstimmt. 

hangen. Dies bezieht sich vor allem^anV!!- 6 ^* 1 ' 6416 ’ die jedoch stark von der Parameter^ 1 
Pajjt man d le Best >mmung raumlicher Abstandsverteilunge* 1 - 

sind und^ ^ nnabme ’ daG die Ergebnisse de ^ u ^igenden Einschatzungen der Modelle: 
konnen be 11111 Aussa S en nur von Sim / Simulati °nen bei hoher EV-WW schlecht relax* 
Aulkn^reich Zen Chr ° m0S0Inen keiLd^ 10 " 6 " mit nied nger EV-WW gewonnen werd«» 
Experimenten c Dynamik Un d Zuganelichk " ° berflachen und zeigen dariiberhmaus 
' dlC Gene "dativ zur Chromo e “ Dies ^iderspricht dem ICD-ModeU““ 

nromosomenoberflache genau lokalisieren. Das 9*° l 
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Model 1 widerspricht durch seine unscharfe Organisation Vorstellungen der gezielten 

Genomregulation, wahrend im MLS-Modell dies durch die distinkten Untereinheiten klar 
gegeben ist. 

Ob die Annahme schlechter Relaxation gerechtfertigt ist, kann bisher nicht mit letzter 
Sicherheit gesagt werden. Da jedoch in Experimente bisher keine solch scharf begrenzten 
Gebilde, wie sie bei hoher EV-WW auftreten, zeigen konnten, konnte dies einerseits ein 
Hinweis auf streBrelaxierende Mechanismen (Topoisomerase II) im Zellkern sein oder 
schlichtweg die Tatsache widerspiegeln, daB die gut relaxierten Simulationen der Natur am 
nachsten kommen. 

Im experimentellen Teil dieser Arbeit wurden dreidimensionale Abstandsverteilungen 
zwischen genomischen Markem gemessen, urn eigene Daten mit den Simulationen vergleichen 
zu konnen. Auch hier wurde ein Feld von Parametern verandert. Die bisher wichtigsten 
gefundenen Abhangigkeiten und Ergebnisse lauten: 

Es war moglich, mit CLSM und schonenden FISH-Methoden kleine genomische 
Sequenzen zu markieren und Abstandsverteilungen mit geringem Fehler zu messen. 

- Zwischen mutterlichem und vaterlichem Chromosom 15 in der PLWS/A-Region gibt es 
statistisch signifikante Unterschiede in den Abstandsverteilungen, die nicht auf einen 
statistischen Sortiereffekt zuruckzufiihren sind. 

" Fur Praparate, bei denen das MeBsignal verstarkt worden war, ergeben sich groBere 
Abstandsverteilungen, aufgrund einer Auswahl beim Aufnehmen der Zellkeme. 

■ Flir kleine genomische Abstande unterhalb von 200kbp kann keine exakte 
Abstandsverteilung bei einfarbiger Markierung der Markerpaare mehr bestimmt 
werden, da die Signale oft nicht trennbar sind. Dies fiihrte zu zu groBen Mittelwerten in 
der Abstandsverteilung. Allerdings ist auch hier noch ein Unterschied zwischen 
mutterlichem und vaterlichem Chromosom erkennbar. 


Ein Vergleich zwischen Experimente und Simulationen ergibt folgendes Bild: 

Die bisher durchgefiihrten Experimente verschiedener Arbeitsgruppen lassen sich durch 
verschiedene Modelle interpretieren. Eine Systematik ist hierbei nicht erkennbar. Die Daten 
von Yokota lassen sich jedoch haufig durch ein RWGL-Modell mit sehr groBen Schleifen 
beschreiben. Hier besteht die Vermutung, daB dies auf Praparationsartefakte zuriickzufiihren 
ist. Die in dieser Arbeit gemachten Experimente favorisieren eher das MLS-Modell. Fur kleine 
genomische Abstande zeigen diese Daten jedoch kein klares Bild, da sie nur einfarbig markiert 
waren. Die Abstande liegen dabei zu hoch, da die Verteilungen zu groBen Abstanden 
verschoben sind. 

Mit einer schonenden Preparation maB Yokota jedoch auch Daten, die einem MLS-Model 
entsprechen. Die entscheidenden Faktoren sind unbekant. 

Diese Interpretation ging davon aus, daB die Simulationen korrekt sind. Umgekehrt konnte es 
jedoch sein, daB die Experimente klares Verhalten zeigen, die Simulationen dagegen 
verbesserungswiirdig sind. Da aber das RWGL-Modell durch die Simulationen korrekt 
beschrieben wurde, ist dies eher unwahrscheinlich. 
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Urn zukunftige Fragestellungen der Zellkernorgamsation zu untersuchen warden fraktale 
AnalyLveSren in den Simulations eingesetzt Hterbe. konnte erstmals fur e.ne zellu,a re 
Orgs sadonsstruktur ihr multifraktaler Aufbau in den Simulationen gezeigt werden. 

Im RWGL- und MLS-Modell konnte geze.gt werden, daB die Bere.che, ,n der e.ne besting 
fraktale Dimension gultig ist, den Organisationsmveaus der entsprechenden Modelle 
entspricht. Bisher resultierten alle Parameteranderungen in den Modellen in emer spezifischen 
Anderung der fraktalen Dimension. Jede Simulation konnte damit gegen andere Simulationen 
abgegrenzt werden. Abhangigkeiten, die sich nicht in der fraktalen Dimension niederschlugen 
miissen mittlerweile mit Vorsicht betrachtet werden. 

Da es fur viele komplexe Systeme - wie zum Beispiel Polymere, Gele Oder Mineralien - gut 
dokumentierte fraktale Analysen gibt, konnten somit zur Interpretation herangezogen werden 
Zugrundeliegende Prinzipien der Organisation des Zellkerns konnten, dabei mit in 
GesetzmaBigkeiten aus anderen naturwissenschaftlichen Bereichen in Zusammenhang 
gebracht werden. 

Ausblick: 

Ausgehend von den Monte-Carlo-Simulationen konnen die Brownschen-Dynamik- 

Simulationen ausgewertet werden. Die Ergebnisse werden wegen der kurzeren Segmentlange 

nicht nur eine bessere Unterscheidung zwischen den Modellen aufgrund hoherer Auflosung 

zulassen sondern auch lokale Abstande zuganglich machen. Damit wird theoretisch eine 

aS u T, d | e , Position v °" Schleifenbasispunkten aufgrund der Form de, 

fKtTrauiwtn 8 f'". M '“ elWen mS * lich ' Sollte " ^ch die Schleifenbasispunkle ak 

H £pe2ifiSChe Se< > ue " z VOTha " de " sein. Eine solche 

Sequenz in den 97 % nicht ”l^d" IC * ® efunden werden - Es ist anzunehmeu, daB eine solche 

durchaefiihtt werden um „ • a en Pe 8 l0nen liegt. Weitere Simulationen miissen 
rchgefuhn werden, urn etn.ge der btsher unentsehiedenen Effekte genauer zu kldten. 

^ Wstanzvertei, ungen getnessen werden, um 

mutterlichem und vaterlichem Chmmr. ; . 1 ^ ei mu ® vor allem der Untschied zwischen 

Zellinie benutzt werden, bei der ein^T™ ^ lrekt nac hgewiesen werden. Hierzu kann eine 

mutterlichen Chromosoms anzeiet Anrh okat,on auf Chromosom 15 die Herkunft des 

um schnellere Kontrolle iiberdie Aufnahme'n zu l ^ I J i a *J riletec * ln ** c ^ ann noch verbessert werden, 




vielleicht auch eine 
erfolgen. 


Fur eine Fortfuhrung der fraktalen a , 

aufwendierBo^v ^? lches sich zwischen^en^h ^ fraktale Dimension des Volumens 
in den sfmni . ■ Verfah ren verwendet werden n hromatinfa sem befindet. Hierzu wird das 
Vorhersaeen ^ UCh eine Mul bfraktalitat Multifra ktalitat der Chromosomen 

r das Diffusionsverhalten im Zellk ** Zw,schenra umes folgt, konnen damit 

" ,,kem gemacht werden. 
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Anhang 


Konfi- 

gurations- 

name 

Schleifen 

GroBe 

Htbpl 

Schleifen 

Position 

Schleifen 

Form 

Linker 

Lange 

[nml 

Linker 

Form 

Excluded 

Volume 

fkTl 

MC 

Schritte 

[*10001 

acl5MLS96 

126 

trans 

loop 

600 

loop 

1.0 

395 

alcl5MLS96 

126 

trans 

lQ OP 

600 

loop 

0.1 

756 

bcl5MLS96 

126 

cis 

loop 

600 

loop 

1.0 

395 

blcl5MLS96 

126 

cis 

loop 

600 

loop 

0.1 

561 

ccl5MLS96 

126 

trans 

loop 

1200 

loop 

1.0 

411 

clcl5MLS96 

126 

trans 

loop 

1200 

loop 

0.1 

530 

dcl5MLS96 

126 

cis 

loop 

1200 

loop 

1.0 

411 

dlcl5MLS96 

126 

cis 

loop 

1200 

loop 

0.1 

426 

ecl5MLS96 

126 

trans 

loop 

1800 

loop 

1.0 

408 

elcl5MLS96 

126 

trans 

loop 

1800 

loop 

0.1 

348 

fcl5MLS96 

126 

cis 

loop 

1800 

loop 

1.0 

408 

flcl5MLS96 

126 

cis 

loop 

1800 

loop 

0.1 

331 

gcl5MLS96 

126 

trans 

loop 

2400 

loop 

1.0 

750 

glcl5MLS96 

126 

trans 

loop 

2400 

loop 

0.1 

335 

hcl5MLS96 

126 

cis 

loop 

2400 

loop 

1.0 

695 

hlcl5MLS96 

126 

cis 

loop 

2400 

loop 

0.1 

270 


Tab.A.l: 

Simula tionsansatze des MLS-Modelis, Feld 1 


Konflg 

uration 

Name 

Schlei 

fen 

GoBe 

fkbp] 

Schleif 

en 

Posi 

tion 

Schlei 

fen 

Form 

Linker 

Lange 

[nm] 

Linker 

Form 

Seg 

ment 

Lange 

[nm] 

Excluded 

Volume 

[kT] 

MC 

Schritte 

[*1000] 

clcl5SmaUMLS96 

80 

trans 

loop 

1200 

loop 

200 

0.1 

321 

dlcl5SmaUMLS96 

100 

trans 

loop 

1200 

loop 

200 

0.1 

443 

elcl5SmallMLS96 

158 

trans 

loop 

1200 

loop 

300 

0.1 

326 

flcl5SmalLMLS96 

251 

trans 

loop 

1200 

loop 

300 

0.1 

420 


Tab.A.2: 

von Simulationszyklus IIL 
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Konfiguration 

Name 

Schleifen 

GoBe 

[Mbp] 

Schleifen 

Form 

Linker 

Lange 

Tnni] 

Linker 

Form 

Excluded 

Volume 

[kTl 

Hvic 

Schritte 

r*iooft]_ 

acl5GIANT2D20 

5.0 

2D-random 

3600 

loop 

1.0 

312 

alcl5GIANT2D20 

5.0 

2D-random 

3600 

loop 

0.1 

574 

bcl5GIANT2D20 

5.0 

2D-random 

3600 

random 

1.0 

291 

blcl5GIANT2D20 

5.0 

2D-random 

3600 

loop 

0.1 

381 

clcl5GIANT3D20 

5.0 

3D-random 

3600 

loop 

0.1 

328 

dlcl5GIANT3D20 

5.0 

3D-random 

3600 

random 

0.1 

285 

ecl5GIANT2D25 

4.0 

2D-random 

3000 

random 

1.0 

374 

elcl5GIANT3D25 

4.0 

3D-random 

3000 

random 

0.1 

466 

fcl5GIANT2D33 

3.0 

2D-random 

2400 

random 

1.0 

800 

flcl5GIANT3D33 

3.0 

3D-random 

2400 

random 

0.1 

347 

gcl5GIANT2D48 

2.0 

2D-random 

1800 

random 

1.0 

543 

glcl5GIANT3D48 

2.0 

3D-random 

1800 

random 

0.1 

214 

hcl5GIANTlD94 

1.0 

loop 

1200 

loop 

1.0 

544 

hlcl5GIANTlD94 

1.0 

loop 

1200 

loop 

0.1 

420 

h2cl5GIANT3D94 

1.0 

3D-random 

1200 

loop 

0.1 

495 

icl5GIANTlD187 

0.503 

loop 

600 

loop 

1.0 

323 

ilcl5GIANTlD187 

0.503 

loop 

600 

loop 

0.1 

281 

kcl5GIANTlD338 

0.251 

loop 

600 

loop 

1.0 

730 

klcl5GIANTlD338 

0.251 

loop 

600 

loop 

0 1 

449 

1c15GIANT1D561 

0.126 

loop 

600 

loop 

V/. 1 
1 0 

175 

Ucl5GIANT1 H461 
Tab.A.3: 

0.126 

loop 

600 

loop 

1 .u 
0.1 

414 


Simula tionsansatze fur das RWGL-Modeli. 


Konfiguration 
Name 


_alcii>(i|25 S0TAK^nQ^ 
_blcl5GI2S SQTAK?no^ 
.LinearlSMI 
J Linear 15 Ml 
JNewLine^rigK/n S96 

Tab.A.4: 


Schleif 

en 

Grofte 

Jkbp] 

_ 126 

Schleifen 

Form 

3D-randnm 

Linker 

Lange 

[nm] 

7z:aa 

Linker 

Form 

Excluded 

Volume 

[kT] 

MC 

Schritte 

[*1000] 

126 ~ 
126 

_3D-random 

J2P-ranrinm 
loop 
_ loop 

_ loop 

jo(J(J 

3600 

random 

random 

0.1 

0.1 

213 _ 

218 _ 

126 

3600 

loop 

0.1 

329 _ 

, 126 

1200 

loop 

1.0 

2675 

126 

1200 

loop 

0.1 

2080 

1200 

loop 

0.1 

486 -J 
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Euklid 

Fraktal 

Fraktal 

Fraktal 

Fraktal 

Max 

Konfig 

Cut-Off 

Zone I 

Zonell 

Zonelll 

ZonelV 

Cut-Off 

u ration 

D 

R 

D 

R 

D 

R 

D 

R 

D 

R 

D 

R 

acl5MLS96 

1,07 

0,999 

2,84 

0,998 

3,68 

0,999 

1.70 

0,997 



4,32 

0,989 

alcl5MLS96 

1,08 

0,999 

2,67 

0,999 

3,13 

0,999 

2,07 

0,999 

2,43 

0,995 



bcl5MLS96 

1,07 

0,999 

2,71 

0.999 

3.64 

0.999 

1,601 

0.998 



6,33 

0,995 

blcl5MLS96 

1,08 

0,999 

2.63 

0.999 

3.07 

0.999 

2,07 

0.999 

2.39 

0,999 



ccl5MLS96 

1,07 

0,999 

2,85 

0.999 

3,56 

0.999 

1.60 

0,998 



4,95 

0.980 

clcl5MLS96 

1,08 

0.999 

2,50 

0.999 

2,89 

0,998 

2,06 

0.999 

2,25 

0.999 



dcl5MLS96 

1,07 

0.999 

2,98 

0.998 

3,66 

0,999 

1,86 

0,997 



5,34 

0,992 

dlcl5MLS96 

1,08 

0,999 

2,50 

0,999 

2,96 

0,999 

2,03 

0,999 

2,15 

0,999 



ecl5MLS96 

1,07 

0,999 

2,7 

0.998 

3,40 

0,999 

1,97 

0,997 



5,37 

0,997 

elcl5MLS96 

1,07 

0.999 

2,38 

0.998 

3,11 

0.999 

2,083 

0,999 





fcl5MLS96 

1,07 

0,999 

2,77 

0,998 

3,19 

0,999 

1,86 

0,998 



4,72 

0,997 

flcl5MLS96 

1,07 

0,999 

2,15 

0.999 

2,95 

0,999 

2,11 

0,999 





2C15MLS96 

1,07 

0,999 

2,51 

0,999 

2,70 

0,998 

2,17 

0,999 





glcl5MLS96 

1,07 

0,999 

2,23 

0,999 

2,49 

0,998 

1,98 

1,0 





hcl5MLS96 

1,07 

0,999 

2,62 

0,999 

2,83 

0,998 

2,12 

0,998 





hlc!5MLS96 

1,07 

0.999 

2,20 

0,999 

2,42 

0,999 

1,98 

0.999 

2,28 

0.986 




Tab.A.5: 


Fraktale Dimensionen mit in ihren Zonen mit Korrelationskoeffizient des Fits fur MLS-Feld I. 


Konfig 

uration 

Euklid 
Cut-Off 
D R 

Fraktal 
Zone I 
D R 

Fraktal 
Zonell 
D R 

Fraktal 
Zonelll 
D R 

c 1 c 15S mallM LS 96 

1,05 

0,999 

2,31 

0,999 

2,96 

0,998 

2,01 

0,999 

dlcl5SmallMLS96 

1,07 

0,999 

2,26 

0,998 

2,92 

0,997 

1,94 

0,999 

elcl5SmallMLS96 

1,07 

0,999 

2,47 

0,998 

3,035 

0,998 

2,08 

1,0 

flcl5SmalLMLS96 

1,07 

0,999 

2,20 

0,999 

2,95 

0,999 

2,09 

1,0 


F raktale Dimensionen mit in ihren Zonen mit Korrelationskoeffizient des Fitts fiir MLS-Feld 2. 
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Konfig 

uration 

Euklid 
Cut-Off 
D R 

Fn 

Zo 

D 

iktal 
»ne I 
R 

Fraktal 
Zonell 
D R 

Fraktal 
Zone III 
D R 

acl5GIANT2D20 

1,05 

0,999 

1,81 

0,999 

2,60 

0,998 



alclSGIANT2D20 

1,06 

0.999 

1,82 

0,999 

2,67 

0,999 



be 15GIANT2D20 

1,05 

0,999 

1,81 

0,999 

2,32 

0,999 



blc!5GIANT2D20 

1,06 

0,999 

1,76 

0,999 

2,45 

0,998 



clcl5GIANT3D20 

1,06 

0,999 

1,83 

0,999 

2,70 

0,999 



dlcl5GIANT3D20 

1,06 

0,999 

1,83 

0,999 

2,78 

0,999 



ecl5GIANT2D25 

1,05 

0,999 

1,76 

0,999 

2,40 

1,0 



elcl5GIANT3D25 

1,05 

0,999 

1,82 

0,999 

2,81 

0,998 



fcl5GIANT2D33 

1,06 

0,999 

1,83 

0,999 

3,05 

0,999 



flcl5GIANT3D33 

1,06 

0,999 

1,89 

0,999 

3,04 

0,999 



gcl5GIANT2D48 

1,05 

0,999 

1,79 

0,999 

2,50 

0,999 



glcl5GIANT3D48 

1,06 

0,999 

1,86 

0,999 

2,99 

0,999 



hcl5GIANTlD94 

1,05 

0,999 

1,78 

0,999 

2,77 

0,999 



hlcl5GIANTlD94 

1,06 

0,999 

1,92 

0,999 

2,63 

0,999 



h2cl5GIANT3D94 

1,06 

0,999 

1,84 

0,999 

2,74 

0,999 



icl5GIANTlD187 

1,06 

0,999 

1,91 

0,998 

2,69 

0,999 

2,36 

0.999 

ilcl5GIANTlD187 

1,07 

0,999 

1,87 

0,998 

2,75 

0,999 

2,47 

0,999 

kcl5GIANTlD338 

1,07 

0,999 

2,19 

0,998 

2,54 

0,999 

1,71 

0,999 

klcl5GIANTlD338 

1,06 

0,999 

2,16 

0.999 

2,43 

0,999 

2,04 

1,0 

1c15GIANT1D561 

1,07 

0,999 

2,76 

0,999 

1,26 

0,999 

1,46 

0,999 

11C15GIANT1D561 

1,07 

0,999 

2,23 

0,999 

1,93 

1,0 




Tab.II.7: 


Fraktale Dimensionen mit in ihren Zonen mit Korrelationskoeffizient des Fitts fiir das RWGL-Model. 


Konfig 

uration 

Euklid 
Cut-Off 
D R 

Fraktal 
Zone I 
D R 

Fraktal 
Zonell 
D R 

alcl5GI850TAK3D32 

1,06 

0,999 

1,80 

0,999 

2,73 

0,998 

alcl5G12550TAK3D96 

1,06 

0,999 

2,04 

0,998 

2,86 

0,999 

_ blcl5GI2550TAK2D96 

1,05 

0,999 

1,96 

0,998 

2,75 

0,999 

!NewLinearl5MLS96 

1,06 

0,999 

1,99 

0,998 

2,69 

0,999 


Tab.II.8: 


ktale Dimensionen mit in ihren Zonen mit Korrelationskoeffizient des Fits fiir andere Modelle. 
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